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In  the  next  several  sections  we  shall  be  considering  a  vari¬ 
ety  of  problem®,  each  associated  with  some  special  geometrical 
configuration. 

In  treating  a  particular  boundary  configuration  one  must 
first  consider  whether  it  is  tractable  at  all  by  the  theory  of 
infinitesimal  waves,  i.e.  whether  it  is  possible  to  select  a  per¬ 
turbation  parameter  t  satisfying  the  requirements  mentioned  in 
section  10.  On  this  basis,  for  example,  it  would  appear  unrea¬ 
sonable  to  try  to  apply  infinitesimal-wave  theory  to  the  waves 
generated  by  a  vertical  circular  cylinder  moving  with  constant 
velocity,  for  the  slope  of  free  surface  may  be  expected  to  be- 
c  >me  very  large  near  the  front  of  the  cylinr’ ^r.  On  the  other 
hand,  in  certain  similar  situations,  notably  the  theory  of  plan¬ 
ing  surfaces,  it  is  possible  to  strain  the  theory  to  accommodate 
such  a  situation.  The  choice  of  par -meter  will  be  discussed  In 
each  individual  case.  We  call  attention  to  the  fact  that  In 
many  cases  it  is  a  consequence  of  the  linearization  procedure 
that  the  boundary  condition  on  a  solid  bounda^  *  is  no  longer  to 
be  satisfied  oa  the  physical  boundary,  but  ins  d  on  some  neigh¬ 
boring  surface.  The  same  situation  occurred  earlier  in  lineariz¬ 
ing  the  free-surface  condition.  This  should  not  be  considered 
as  a  further  approximation,  but  rather  as  one  consistent  with 
the  inf inltesimal-wave  approximation. 

The  methods  for  finding  a  solution  to  a  boundary-value 
problem,  once  it  has  been  properly  formulated,  seem  to  fall  into 
two  or  possibly  three  groups.  One  method  is  a  combination  of 
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separation  of  variables  and  expansion  of  the  factors  in  Fourier- 
type  series  or  integrals.  This  requires,  of  course,  a  geometric 
configuration  related  in  a  suitable  way  to  the  coordinate  sur¬ 
faces  of  a  set  of  variables  which  allows  separation  and  a  complete 
set  of  associated  elementary  solutions  to  be  used  in  the  expansion. 
If  a  Fourier-series  expansion  is  to  be  used,  orthogonality  of  the 
elementary  solutions  is  desirable. 

If  the  motion  is  harmonic  in  time  with  frequency  er  and  if 
the  fluid  Ip  of  finite  depth  h,  then  the  functions 


^cosk  (  j 


(  "  0  ^  (16.1) 
occurring  as  factors  in  (13.2)  and  (13,4),  in  (13.6),  .  nd  in 
(13 JO  may  be  shown  easily  by  direct  computation  to  be  orthogonal 
on  the  interval  0>y>- It .  Weinstein  [1927]  has  shown  that  they 
form  u  complete  set  on  this  interval.  The  result  may  be  used  in 
the  following  way,  for  example.  Suppose  fluid  occupies  the  region 


and  that  the  boundary  conditions  on  the  wallr  and  bottom  are 


(fjo.'H,*)*  (rJvW'rt. 

♦j.  ( j'i ’j.’rt  - 

4v(x,y,-k,t-)=o. 

Then,  by  expressing  Ff^  ,  J.'j  as  a  double  series 


(16.2) 
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(16.3) 


(with  appropriate  restrictions  upon  F),  one  may  construct  a 
solution  from  the  elementary  solutions  in  (13.6).  Further  con¬ 
ditions  relating  to  boundedness  and  behavior  as  /  H»  °Q  are 
necessary  in  order  to  ensure  a  unique  solution,  but  will  not  be 
discussed  here.  The  ilementary  solutions  (13.8)  can  be  used  in 
a  similar  way  for  the  region  exterior  to  a  vertical  cylindrical 
boundary.  Still  other  configurations  are  possible  corresponding 
to  the  various  coordinate  systems  allowing  separation  of 

If  the  fluid  is  infinitely  deep,  it  is  possible  to  construct 
a  Fourier-integral  expansion  using  the  functions 

I 

|e.  ^  A  +  v /<*]  j  v-G j  o  <•  A  <  °o ,  16  4^ 


In  fact,  Havelock  [1929b]  has  remarked  that  the  usual  Fourier- 
integral  representation  of  a  function  may  be  altered  to  give 

t  W  “  7 r  f  '  Vv  ' 


Jr  1V  £V^  _[  *f  (*))  e  ^ 


(16.5) 


If  the  problem  is  such  that  rectangular  coordinates  may  be  used 
conveniently,  then  (16.5)  may  be  combined  with  a  Fourier-series 
or  Fourier-integral  expansion  in  y  and  the  elementary  solutions 
(13,5)  used  to  construct  a  solution  analogous  to  (16.3).  The 
necessary  expressions  in  both  rectangular  and  cylindrical  coordi- 


nates  can  be  found  "In  the  cited  paper  of  Havelock. 

If  the  fluid  is  of  bounded  horizontal  extent,  and  la  bounded 
by  vertical  surfaces  which  are  constant-coordinate  surfaces  in 
on©  of  the  coordinate  systems  allowing  separation  of  d^Cp  **  <M 
the  various  possible  modes  of  motion  of  the  fluid  may  be  obtained 
as  the  solution  of  an  eigenvalue  problem  of  a  classical  type.  If 
the  container  is  of  more  general  shape,  it  is  more  difficult  to 
obtain  explicit  solutions.  The  problem  will  be  discussed  in  sec¬ 
tion  23. 

The  orthogonal  functions  (16.1)  were  associated  with  a  single 

value  of  the  frequency  O'.  It  possible  to  derive  another  re- 

t 

suit  concerning  orthogonality  of  solutions  associated  with  differ¬ 
ent  values  of  '0*  Let  d*,  t  and  <•£][  (*«  'j  ®tt*j 

be  regular  ve loci ty  potentials  of  harmonic  oscillations  of  differ¬ 
ent  frequencies.  Furthermore,  let  any  solid  boundaries  be  fixed 
and, if  the  fluid  is  not  bounded  in  extent, we  suppose  that 
|  grad  cp|  =  QfR‘  "‘J  a®  7?"=  xV^-aoo  .  Consider  the  fluid  contained 
within  a  1.  ‘ge  cylinder  of  radius  "R  and  above  the  plane 

Hr-'R  .  The  fluid  will  be  bounded  partly  by  free  surface 
p^  »  partly  by  solid  boundaries  ,  partly  by  the  horizon¬ 
tal  plane  J3^  and  partly  by  the  cylinder  -H-r  .  Applying  Green's 
theorem  to  the  two  potential  functions,  one  obtains 

O  =  (<P,  <4*  = 

^ 1  (?'  ^  f  j  t  H  (f.  V-  <0  J  r. 

Fr  Sr+Jlft 


(16.6) 
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As  I?"**0  ,  the  integral  over  A  R  "f  ~3J>0,  and  on©  has 

f  \  ^  C°* 


(16.7) 


Pron  th©  tree-surface  condition 


,16.8) 


and  (10.7)  becomes 


ff.V-ir’-V  J"j  f,  fJC.o,  j'J  0,  j)  i<r  -  O, 

3-  f 


( 16,9) 


or  simply 


S5  <p.f«*d)9UJ'<M) '«-=o. 


(16. 10) 


Hence  and  are  orthogonal  over  the  free  surface  of  the 

fluid.  This  theorem  can  be  used  for  certain  initial-value  prob¬ 
lems  in  a  manner  analogous  to  that  in  which  the  orthogonality  of 
(16.1)  can  be  used  for  boundary-value  problems.  This  will  be 
done  in  section  22. 

A  second  method  for  solving  special  problems  is  the  method 
of  Green's  functions  or  source  functions  (cf.  Volterra  [1934)). 

In  this  method  one  constructs  first  a  potential  function  of  the 
form 

?•?,!)*  r  ■» 6. (*.?■}. if.?,.?), 

(16.11) 


rv*  (f-rjv  *■  (3-7)'’+  (j- 
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such  that  G0  is  regular  In  ^  <0  and  such  that  Q  satis  flea 
the  free-surface  condition,  condition®'  at  Infinity  appropriate 
to  che  problem  at  hand,  and,  If  the  fluid  Is  of  finite  depth, 
the  boundary  condition  on  the  bottom.  Such  solution®  are,  of 
course,  just  the  singular  solution®  derived  in  section  13^  . 
lfext,  if  there  are  surface®  S  in  (or  on)  the  fluid  upon  which 
certain  further  boundary  conditions  must  be  satisfied,  we  attempt 
to  satisfy  them  by  a  distribution  of  the  modified  sources  (16.11) 
over  the  surface (a)  S  * 


s 


(16.12) 


Hen  \f  i*  an  unknown  function  which  is  to  be  determined  from 
the  boundary  condition  on  S*  1®  most  problems  this  boundary 
condition  consists  in  specifying  (p^  on  5.  fell  known  proper¬ 
ties  of  surface  distributions  of  sources  then  allow  one  to  formu- 
A  ©A©  mi  An  m  1  ©(juntt  A  oil  for  i  « 


&(/,),$# » (i6. 13) 

(y0'i)  *" 

where  Y]  Is  the  exterior  normal  to  the  surface  S  (taken  here  as 
a  closed  surface).  When  it  is  convenient,  one  may  also  use  dis¬ 
tributions  of  dipoles. 

It  is  also  possible,  and  sometimes  advantageous,  to  construct 
solutions  satisfying  given  boundary  conditions  on  a  closed  surface 
5  distributing  the  singular  solutions  on  surfaces,  lines  or 


points  completely  inside  $•  Examples  will  occur  later. 

A  third  method  of  approach  is  to  seek  first,,  instead  of 
or  ,  the  functions 

These  functions  satisfy  a  simpler  condition  on  the  plane  Jj-O: 

0  ^  'R'e  ^  O- o,  ^  -  0 . 

If  the  other  boundary  conditions  are  such  that  they  can  be 
formulated  simply  in  terms  of  ^  or  ^  ,  the  new  problem  may  be 
simpler  to  solve.  After  finding  ^  or  ^  ,  one  must  then  solve 
a  differential  equation  in  order  to  obtain  the  desired  solution 
(p  or  .  This  procedure  is  called  the  "reduction"  method  by 
Weinstein  [1949],  It  was  apparently  firsx  introduced  by  Keldysh 
[19351  and  has  since  been  much  exploited  by  him,  Kochln,  Sedov, 
Haskind,  Lewy,  Stoker  and  others.  It  has  already  been  used  in 
the  derivation  of  (13.28)  and  will  be  applied  in  several  other 
problems.  (The  method  is  used  also  by  Uuskhelishvili  [Singular 
Integral  equations,  Noordhoff,  Groningen,  1953,  §74]  to  reduce 
a  mixed  boundary  condition  of  more  complicated  type  to  a  simple 
one.)  The  solution  of  the  reduced  problem  may,  of  course,  be 
carried  out  by  one  of  the  two  methods  already  described  above, 
or  any  other  one  which  is  convenient. 

The  methods  outlined  above  do  not  exhaust  the  possible  ones 
for  finding  analytic  solutions.  However,  tnoy  will  occur  fre¬ 
quently  in  the  next  several  sections.  Seve  al  of  the  special 
problems  treated  in  the  following  sections  can  be  solved  by  each 
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of  the  three  approaches.  The  choice  of  a  particular  one  has  been 
made  either  to  illustrate  a  method  or  because  it  happens  to  be 
convenient.  Techniques  for  finding  numerical  solutions  will  not 


be  discussed. 


IT •  Two-dimensional  progressive  and  standing  wave a  In  unbounded 
ggg.1 1 -boundir  1  e  a  ■ 

In  this  and  the  following  section  we  shall  consider  situations 
in  which  the  region  occupied  by  fluid  extends  to  infinity  horizon¬ 
tally,  the  solid  boundaries  are  fixed,  but  of  more  complicated 
shape  than  the  simple  flat  bottom  considered  up  to  now,  and  the 
motion  of  the  fluid  at  infinity  is  prescribed,  or  at  least  partly 
so*  We  shall  assume  that  the  velocity  is  bounded  at  all  interior 
points  of  the  fluid  and  also  at  the  Infinite  1 imi te  of  the  fluid. 

The  motion  lu  taken  to  be  periodic  everywhere  with  period  &  .  Hence 
we  shall  assume  (cf .  §  11)  that 

${*, J, t)=  +  *  1?e  <pe 

The  restriction  to  standing  or  progressive  waves  can  be  properly 
applied  only  at  x  ■  +  ».  Thus,  we  shall  look  for  solutions  which 
at  x  «  ■  behave  like 


or 


(A  cos  mx  ♦  B  sin  mx)  cos  t 


A  cos  (mx  +  6  t)  +  B  ccs  (mx  -  d  t), 
and  similarly  at  x  ■  -«  if  the-  fluid  extends  In  that  direction.  As 
we  shall  see  below,  the  coefficients  cannot  be  chosen  independently 
If  f  remains  bounded  everywhere. 

The  parameter  of  linearization  may  be  chosen  as 


£  »  max  (Am,  Bm) . 

If  the  solution  *p  is  bounded  everywhere,  then  as  £  ■»  0,  If  •*  0 
uniformly.  However,  if  a  singularity  Is  allowed,  then  <f  +  Q 
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uniformly  only  in  a  region  excluding  a  neighborhood  about  the 
singularity .  One  may  presume  that  the  solution  to  the  linearised 
problem  loses  physical  significance  within  such  a  neighborhood. 

(It  is  assumed  tr;  Stoker  ,[l947j  p.  5j  that  singularities  at  the 
surface  are  associated  with  breaking  of  the  waves.) 

We  shall  discuss  below  two  types  of  problems:  obstacles  In 
an  infinite  ocean  and  sloping  beaches.  For  each  type  a  special  case 
will  be  discussed  in  some  detail. 

17  <x  .  Obstructions  in  an  infinitely  long  canal. 


Consider  first  the  following  situation.  The  fluid  extends 
from  x  »  ••  to  x  ■  +«j  the  bottom  is  given  by  y  ■  -h(x),  where 


h'x)  «  h,  >  0  for  x  >_  x., ,  h(x) 


h2  >  0  for  x  <_  x2;  fixed  obstacles 


/  W  4.  W*  A  /  Aj 

may  be  present  in  the  fluid  or  on  the  surface  (see  Figure  13).  The 
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Figure  13. 

surface  at  x  ■  +•  will  be  assumed  to  behave  like 


^  »#  |  CQ&  (mtx>  "f*1  f  )  *f*  t  p,) 


and  it  x  m  -•*  like 

?  •  |  c os(n\JC  ts?t  #«wt)  +  BgCOs  (tV%.x  -  di  t  fiz) 


h  proof  of  the  existence  of  a  solution  to  this  problem  does  not 
seem  to  exist  for  the  general  case.  One  would  not  expect  a  unique* 
at'*#  theorem  since  no  statement  has  been  made  about  singularities 
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or  circulation.  For  Infinite  depth  and  a  submerged  body  Koehlnj[X959] 
has  proved  the  existence  for  sufficiently  large  values  of  m  (the 
situation  1®  slightly  different,,  but  the  proof  carries  over),  Krelsel 
(1949J  has  established  the  existence  of  a  solution  and  Its  uniqueness 
In  two  e».ses,  In  the  first  case  *  hg,  only  obstacles  on  the  bottom 
are  allowed,  <fi  is  assumed  bounded,  and  a  certain  constant,  defined 
in  terns  of  the  conformal  mapping  of  the  fluid  region  onto  the  strip 
0  <  y  <  ^1*  b®  l®8®  than  1*  Included  are  theorems  comparing 

the  values  of  this  constant  for  differ* nt  types  of  obstructions. 

The  second  result  allows  a  shallow  obstruction  in  the  surface,  but 
requires  a  flat  bottom  and  sufficiently  long  waves  and  again  bounded! 

4  *  .Roseau  [ 3 952}  has  proved  existence  and  uniqueness  for  no 
obstructions  within  the  fluid,  but  for  4  hgj  the  curve  Joining 
the  two  ends  is  of  a  special  sort.  John  [i960,  p.  78  ff , ]  has  proved 
uniqueness  lor  a  flat  bottom  and  for  a  body  in  the  free  surface  with 
the  property  that  every  vertical  line  Intersects  either  the  free 
surface  or  the  body  just  once;  certain  regularity  properties  of 
piust  also  be  assumed.  If  the  body  la  convex  and  Intersects  the  free 
surface  perpendicularly,  he  is  able  to  prove'  also  existence  of  a 
solution . 

Existence  and  uniqueness  theorems  have  also  been  proved  for 
several  special  configurations,  in  most  of  these  cases  explicit 
solutions  are  given.  A  vertical-line  barrier  extending  from  the 
free  surface  to  a  depth  /  in  an  infinite  fluid  has  been  considered 
by  Dean  [1948],  Ursell  [l<H7]  ,  and  Haskind  [1948],  Both  Dean  and 
Unwll  rise  consider  a  barrier  extending  from  to  a  distance/ 
below  the  surface.  John  [19483  has  generalized  both  these  problems 


to  the  case  of  a  slanting  barrier  of  slope  ir/ 2n ,  and  obtained  a 

more  general  solution  even  for  the  vertical  barrier.  Dean  £l948j 
and  Ursell  [i960]  have  also  considered  submerged  circular  cylinders 
in  an  infinitely  deep  fluid ,  and  Ursell  has  established  a  uniqueness 
theorem  for  this  case .  A  horizontal  obstruction  of  finite  width  on 
the  water  surface  (the  "finite-dock  problem)  has  been  treated  by 
Rubin  [1954J  ,  who  proved  existence  of  a  solution  by  a  variational 
method.  Other  references  concerning  the  dock  problem  will  be 
given  below. 

iteflectionandtran8misslon^  If  one  assumes  the 

existence  of  a  solution  to  the  general  problem  stated  above,  one 

may  establish  the  form  of  the  solution  for  x  >  and  x  <  xg,  by 
using  Weinstein's  theorem  [ 1927J  asserting  the  completeness  oi  the 
functions  [cf.  (13.6),  reduced  to  two  dimensions^ 

[  cosh  b^(y+h) ,  cosmn(y+h)] 

in  the  interval  -h  <  y  <  0  [cf .  Kreisel,  1949 ,  PP  26-29*  John, 

1948,  p.  152].  it  is 

i,*)m  [Ai  f  t*i)+  £.  cosfrri^Jc  -dt  ccnh  h^ji- 

ao  (17.1) 

+n?/aw“4firf  +  kn^ei*)  tx-pt'^nlxDoDsmyCythi), 

where  1  -  1,2  and  £  2  -  gm^tanh  —gm^tan  . 

Let  us  now  apply  the  formula  for  dE/dt  in  equation  (0.2)  to 
the  region  of  fluid  bounded  by  the  planes  x  »  Cg  <  Xg  ,  x  •  c^  >  x-^, 
the  bottom  and  any  other  obstructions,  which  we  take  to  be  between 
these  two  planes.  Itien,  if  I  x  +  f y  la  boumle<a  111  fcl1®  region 
considered. 


At  -i,,  -\ 

since  on  the  "physical"  boundaries  [cf.  (8. 3)1  either  p  ■  0  or  $n~0. 
Anticipating  that  we  are  interested  only  in  the  asymptotic  values 
for  Cj  •*  •  and  Cg  ♦  ,  we  compute  the  above  expres  ion  using 

only  the  first  term  in  (17*1)  and  average  over  a  period  2  X/xr  i 


-Inn: 


Since  the  average  energy  in  the  region  is  constant, 


ACu  j [_ak (^-2) 


This  is,  of  course,  a  statement  of  the  conservation  of  energy.  If 
A1  is  given  flQ  and  Bg  ■  0,  then  Ag  «  are  uniquely  determined. 
For  suppose  two  solutions  ^  and  tjt*  re  possible,  both  with  the 
same  A^  and  Bg  -  0,  but  one  with  Ag,  B^ ,  the  other  with  Ag,  B^. 
Apply  (17.2)  to  the  difference  f*  -  t*'' 
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Each  side  must  be  zero  since  they  differ  in  sign  and  are  equal. 
Hence  Ag  -  A^  and  B1  -  Bn' 
uniqueness  of  ^  Itself. 


Hence  Ag  m  a^  and  B^  -  B^.  IMs  does  not,  of  course,  imply  the 
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ftp? 


If  -  hgj  then  (17 *2)  simplifies  In  an  obvious  ways 


4®  -  Bj  .  -  B? 


(17.3) 


Here  h  may  also  be  infinite. 

Setting  Bg  -  0  and  fixing  Ax  as  above  corresponds  physically 
Bl-vlng  the  amplitude  of  an  incoming  wave  far  to  the  right. 

Is  then  the  amplitude  of  the  reflected  wave  and  Ag  of  the  transmitted 
wave. 


Th  theorem  of  the  preceding  paragraph  states  that  Aj  fixes 


them  uniquely.  We  '>eflne  j  B^/A^ 1  as  the  reflection  coefficient  R 
and  /  Ag/i  /  as  the  transmission  coefficient  T.  They  are  uniquely 


determined  and  »  1 .  Properly  one  should  define  both  left 

and  right  coefficients  since  the  channel  is  not  symmetric.  However, 
the  uniqueness  theorem  implies  that  both  have  the  same  value.  One 
can  clearly  arrange  the  phases  bo  that  A1  and  have  the  same  sign. 
If  this  is  done,  2  "  **  j  will  oe  the  phase  shift  caused  by  the 
otofl  trlStC  0  # 

Kreisel  [19493  has  proved  several  general  theorems  concerning 
the  reflection  coefficient  if  h1  -  hg.  In  particular,  if  there 
are  no  obstacle®  within  the  fluid,  he  determines  upper  i nd  lower 
bounds  for  the  reflection  coefficient  in  terms  of  the  conformal 
i  mapping  z(  )  of  the  infinite  strip  0  >  -h  onto  the  region  occu¬ 

lt  pied  by  fluid,  with  Infinities  corresponding .  KIs  bounds  become 
"t  closer  as  the  wave  length  increases.  He  gives,  for  example,  asymptotic 


j 

i 


expressions  as  mQ  0  for  the  reflection  coefficient  from  a  horizon¬ 


tal  reef  of  width  a. and  height  €  and  from  a  flat  plate  in  the  sur¬ 
face  of  beam  b,  namely. 
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An  Interesting  special  result  of  Dean  [1947]  (see  also  Uroell 
£l95o])  Is  that  the  reflection  coefficient  from  a  submerged  circular 
cylinder  In  Infinitely  deep  water  vanishes.  The  proof  may  be  briefly 
sketched.  Let  a  be  the  radius  and  let  the  1  jnter  be  at  fO*-b),  b  >ft. 

Let  the  velocity  potential  be  written  as  a  sum  of  an  Inc,. . »ng  wave 

and  a  diverging  wave:  11,1 


f  »  A  vlVi2  00s  ( vX  4  <£  t )  +•  $0  j 


and  suppose  that  can  be  expressed  aa  a  sum  of  multipoles  (13.31), 

starting  with  dipoles: 


i* m' ^  ^h-  (•*•>]>*)  +  4*,  %>*)+  c+sfZ  (<*,*}  ,t  +  f4)  +  («*.  y>t+  f~t)  t 

lfs» 

where  fn  la  the  potential  for  the  symmetric  potential  of  order  n 
and  strength  Q  «  1,  and  that  for  the  antisymmetric  one.  The 

boundary  condition  on  the  cylinder  (using  the  notation  of  (13.31)), 

~j^T~j  ^  11  A  Vtf"1,  V  ^  -cqs(vq.cqs &)cc's  e»4  trt 

t  Si*,  (V  AS***  &)C4M&J 

gives  the  relation  &n  ».  -dn,  bn  *  cR.  The  reflected  wave  at  +• 
from  the  antisymmetric  functions  then  just  „anoela  that  from,  the 
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symmetric  functions.  The;  reinforce  each  other  at  x  ■  The 
phase  change  for  b/a  -  5/4  ,  <**&/&  m  4/3  was  computed  numerically 

by  both  Dean  and  Ursell  and  for  this  case  was  very  close  to  90*. 

Aa  mentioned  above,  straight-line  barriers  have  been  considered 
by  Dean  [1945]  ,  Ursell  [1947]  ,  Haskind  [lQ48'J ,  John  [1948J ,  the 
latter  having  treated  also  barriers  inclined  at  an  angle  JT /2n,  and 
Levine  [l958] .  The  last  three  authors  use  the  reduction  method, 
whereas  the  first  two  use  a  Fourier- integral  method  which  leads  to 
a  singular  integral  equation.  We  shall  treat  this  problem  by  the 
reduction  method. 

Vertical  barrier.  Let  the  barrier  extend  along  the  y-axls 
from  y  ■  0  to  y  ■  -  /  and  suppose  an  Incoming  wave  is  given  at 
X  m  +«  as 

*1*  A  on  (v  X  i-  4t  +  .  / 

We  shall  look  for  a  velocity  potential  f  having  the  form 

and  satisfying  the  following  boundary  conditions  on  the  free  surface 
and  the  barrier: 

fti  f  9  ^  (Wi*),  /-*/>*,*"  J  $  (Q 

As  x  •*+*»,  cp^cos  ert  +  sin  «rt  must  represent  outgoing  waves. 

In  the  neighborhood  of  (0 ,-Jl  )  it  will  be  assumed  that 

In  '#he  neighborhood  of  the  intersection  of  the  barrier  and  the 
surface  (0,0)  as  well  «a  in  the  region  of  fluid  bounded  away  from 
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cite  barrier,  we  a  ha  11  assume  ^  bounded.  It  should  be  noted, 

however,,  that  this  assumption  excludes  a  large  class  of  solutions 
of  possible  physical  Interest  [cf.  John,  1946). 

If  we  Introduce  the  stream  functions  'Jp  ,  IP,  ,  and  *2- 
corresponding  to  /  .  rf(  and  ^  and  the  corresponding  complex 

p Ot» €111. jLul®  F  |  Mid  t  'WWS  ]flAV4t 

F  «<•  (- -~5-i  t  ^  Oadt  *  ^  ^  +fz)  ^  m  Ft  c#a  &t  +  f 

and  the  boundary  conditions . 

Ref-vPd^Fl,!  -  o  ,  =  o  „  \  x  \  >  o  mu,  i,; 

f?t  F„/  ,  X-O,  0  >  *  >J  ^  n:  1,1. 

After  finding  and  F2  satisfying  these  conditions,  constants 
occurring  in  the  solutions  must  be  adjusted  so  that  and 
satisfy  th  radiation  condition: 

(f/±vK)  =o  ,  h 


X  -* 

Consider  the  function 


X  -%  -at  OO 


(h'  *  ^ 


3,  =.  f/4<vF(  -  sT"^  (VV*F, )  . 

Then  the  boundary  conditions  imply  that  satisfies 

Gt  =©  -for  l*l»o, 

U  G/  *-0  F*r  >-9,  ixUo  x-  o  j  o 

The  function  may  be  extended  into  the  upper  half*-plane  by 
defining  x  <S,  £  ">  0.  Since  we  have 

assumed  |F/  |  ~  3  I-})**  4  j  we  may  conclude 

that  1 6*  |  <J3”+Clj|  for  |  'J.  |  >  ^  and  expand  G^  in  a  Laurent 
series 
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where  all  coefficients  are  real  since  Ira  fi^x  4*  10)  »  0.  The  con- 
ditlon  |  Pi  |  *  0  as  y  -  -  Implies  J  0£ »  -  0  as  y  -  -«  and  hence 
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0. 
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Further,  we  nay  show  as  follows  that  »  0.  Consider  a  contour 
containing  the  obstruction  and  lying  in  the  region  |  s )'  >  b.  Then 
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Let  this  contour  b©  contracted  onto  the  barrier.  Then,  from  the 
assumed  behavior  of  F,  on  the  barrier,  the  Integral  vanishes)  hence 
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a,  «■  0.  Thus 
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Let  us  now  exploit  the  boundary  condition  for  Ql  by  mapping  the 

■«ys < 

z -plane  Into  a  j -plane  by  the  mapping 
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where  the  branch  of  the  square  root  *.s  chosen  which  makes  J 
for  large  2“  .  This  maps  the  z -plane  cut  from  ~i£  to  it  9 
along  the  barrier,  onto  the  z -plane  cut  from  -  l  to  +  t  ,  with  In 
flnltles  and  upper  and  lower  half-planes  corresponding.  Ihen 
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and 


Ira  Hj(  j?  +  10)  -  0.  Since  ( jf  )  mi 
HAP  must  have  the  form 
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The  condition  on  &  near  the  edge  of  the  barrier,  Implies  that 
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or 
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and  hence  that  b,  «  0,  n  >  4.  Thus 
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Integrating,  and  writing  D1  *  U 
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where  the  constant  of  Integration  has  been  chosen  so  as  to  make 
G^(z)  behave  like  z"^  for  large  z.  Then 

F  (V-  *  e 1 ' 'l.  A £' ‘  ^ e 

where  the  path  of  Integration  will  be  taken  around  the  right-hand 
side  of  the  barrier.  The  boundary  condition  Re  p^O  +  iy)  -  0, 

0  <  |y  I  <  i  ,  relates  Ex  and  a3  follows.  Prom  P1(z): 

F/(V e'^t iv£’ +  D>  -  ‘ d}  l 


Take  the  path  of  integration  along  the  y-axls,  so  that  the  integral 
becomes 
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We  now  compute  the  asymptotic  expressions  for  F-.fzJ.for  x  +  •  • 

JL 

If  the  path  of  integration  Is  taken  on  a  large  arc  of  radius  R  in 
the  first  quadrant  and  then  to  z,  and  if  R  is  allowed  to  become 
infinite,  it  follows  from  Jordan’s  lemma  that  the  Integral  may  also 


be  written 
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Clearly, 
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whelm  the  path  of  integration  passes  below  the  barrier.  By  couple 


ting  this  path  toy  a  large  semicircle  in  the  upper  half-plane,  which 


gives  a  mt'Q  contribution  in  tb*  limit,,  and  than  contracting  the 
contour  aoout  the  barrier,  one  sees  that 


Hence 

Ff(^)  ^  c  I  [c,  1 1*  tHi  l<t (i>  t)  -  2.W  D,I,  ( v-t  J  * 


Similar  expressions  hol<*  for  F0(z)  with  constants  e0  and  D0* 
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For  f1  and  f0  ’#©  ha\#©  the  asymptotic  expressions 5 
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The  radiation  condition  gives  simultaneous  equations  for  the  deter¬ 
mination  of  ©2#  J>1  and  Dp.  The  solution  may  be  written 
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Substitution  In  the  expressions  for  f,  and  fg,  and  computation 
of  P,  cos  <r t.  +  p_  sln<  o' t  gives,  afttr  a  somewhat  tedious  cal- 

X  Iff 

culatlon,  the  following  asymptotic  expressions  for  ci  : 
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Lj  'being  a  Struve  function  of  Imaginary  argument  [illation,,  p»  329 S 
Lj  Is  tabulated  in  J.  Math.  Phya.  £$  (1946),  252- 269 ]  •  Let  the 
force  and  moment  about  the  origin,  per  unit  length  of  barrier,  be 
denoted  by  X  and  1,  the  former  being  positive  If  directed  along  OX, 
the  latter  la  eountereloekwls  .  Then 

X  *  iZf^At  X®  0&  (let  +o(  4-  fit#,)  } 

fi  »  t  2  ?  j  A  l  v  Ha  on  ( 4.  t  foe  ) ,  (17 .7  j 


where 
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Hasklnd  also  computes  the  average  force  and  moment  per  unit  length 
of  the  barrier.  The  result!  ares 
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where 


I*)  0 


Figure  IS  displays  all  four  functions  in  dimensionless  form. 


The  method  of  Integral  equations .  This  method  for  finding 
solutions  has  been  frequently  used,  especially  by  Kochin  [1937,  1939, 
1940J  and  his  colleagues.  One  of  its  advantages  is  that  approximate 
solutions  to  the  integral  equacion  can  frequently  be  found  even  when 
an  explicit  solution  cannot  be  easily  obtained.  The  following  expo¬ 
sition  follows  approximately  Kochin  [ 193?]  and  Keldysh  and  Lavren- 
t'ev  [1937] . 

Consider  a  submerged  obstacle  whose  contour  C  is  given  para¬ 
metrically  by  z  m  z(s)  and  is  oriented  counterclockwise.  Let  /S  (») 
be  the  angle  between  the  tangent  vector  and  the  positive  x-direction 
(see  Figure  16) .  We  shall  assume  that  as  x  *  «  the  motion 
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The  other  boundary  conditions  In  terns .Of  the  complex  potential 
f(ss)  »  cp(x,yj  +■  1  Y  (x,y)  are 

•*/  (17.19) 

u.  1/  7  -<?* 

4  ♦  —  «o 

Write  f(«)  In  the  form* 


f  (i  )  ■  fj  f:5 )  t 
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Then  f.(*)  oust  satisfy 

/.  (5)  *  0 

X  -*■  (17-12) 
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as  well  as  the  free  surface  condition  and  the  condition  as  y  •+  -■*. 

We  shall  try  to  express  f 1 (* )  as  a  distribution  of  vortices 
over  the  contour  c.  However,  the  vortices  are  chosen  so  that  the 
conditions  on  the  free  surface,  at  x  «■  *»  and  at  y  •  -•*  are  satisfied. 
Aa  is  apparent  from  the  derivation  of  (13.28),  the  complex  velocity 
potential  for  aweh  vortices  Is  given  by 

•o  j 


(17.13) 


We  set  P  «  1  and  try  to  exp. esa  f1(i)  as  follows: 


i(}?-frvf,().j(*)d*  ,  (1T.U) 

where  f  (a)  must  be  chosen  so  that  the  boundary  condition  on  the 
body  is  satisfied. 

In  order  to  derive  an  Integral  equation  for  f  (a),  consider 
the  following  expression  for  fj(z),  a  direct  consequence  of  Cauchy's 
Integral: 
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The  function  g-^z)  Is  regular  everywhere  outside  C1  and  g^z)  la 
regular  everywhere  Inside  C2«  One  may  contract  C^  onto  C  and  extend 
Kgf2)  analytically  Into  the  whole  lower  half-plane  (or  fluid  strip 
If  the  depth  is  finite). 

Consider  now  (for  Infinite  depth j  the  finite  depth  case  Is 
analogous)  the  following  function: 

$V  ‘z7i  [ +~zJT[f'(l)( -^-nv6tV,|}j£L_ 

Ci  i  C  i  ’  *© 


The  first  summand  is  identical  with  g^z)  and  the  second  is  also 
regular  in  the  whole  half-plane.  g(z)  satisfies  the  same  boundary 
conditions  as  fj(z) .  Hence  f-[(z)  -  g(z)  is  regular  in  the  whole 


lower  half-plan®,  satisfies  the  free-surface  condition  and  vanishes 
as  y  ■*  -•  and  x  *  +»•  Ttte  uniqueness  argument  used  In  the  cerlva- 
tlon  of  (13.28)  shows  that  f '  (z)  «  g(z).  Thus  we  have 

f  (j)  (T) {tts- 2t»e  v*  f~fc. ry 

OX  'c,  '  *  -  (17.15) 


Mow  contrAot'  tjo  C  •  Then 
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If  one  substitutes  above,  one  flnda  that  the  contribution  from  the 
second  summand  in  fj(£  )  vanishes  and  that,  since  dT/da  -  e1  P  (8), 


-  [V*M  fa' (y,:w)d*. 


(17. IT) 


This  identifies  J-  (s)  as  the  tangential  velocity  vfc(s)  at  a  point 
of  the  contour. 


Let  us  now  consider  the  effect  of  letting  z  ■*>  z(s'),  a  point 
of  the  contour  C.  Then,  according  to  the  Theorem  of  PleraelJ- 
Sokhotskll, 
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Hanoe  we  have  the  integral  equation  for  {#)s 

-4*W  *  p/rw /  V*  jwV^< (17.19) 

This  is  really  two  integral  equations.  The  imaginary  part  gives  a 
singular  integral  equation  of  the  first  kind: 

p¥(  [(*)K(a\ a)M  -- 2jt^  6  *  cm  [■vx(a‘) -ft(A') -t  ©cj  (17.20) 


rrn-.j  x-jal  part  gives  a  Fredholm  equation  of  the  second  kind  with 
continuous  kernel: 

-i  rM+  ^  fr <*>L(*l*)di .zi iAi  tT'} *.*%>*(*') -w)  <-°<] .  (17 >21) 


Here 
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The  kernel  K(s'#s)  is  of  the  form 


K(s',s)  - - i -  +  C ( B ' , a ) 
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where  u(s',®)  la  continuous;  the  first  tern  cones  from 

el  A (a* )/[z(u> )  .  z(a)].  if  the  curve  C  ia  sufficiently  smooth, 
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where  f  (a)  la  the  radius  of  curvature  of  C  at  z(a) . 

If  the  obstacle  constat#  of  a  smooth  arc,  an  analogous  argument 
leads  to  only  the  singular  Integral  equation  above,  but  with  (a) 
now  Identified  with  the  Jump  in  ( s )  as  one  goes  from  the  left  to 

the  right  side  of  the  arc. 

There  does  not  seem  to  be  a  published  proof  that  a  solution 
to  either  lntegr  -1  equation  exists  for  all  i>  .  However,  Kochln 
[1936,  pp.  119-126]  shows  the  existence  of  a  solution  for  both 
sufficiently  large  and  sufficiently  small  values  of  v  for  the  equa¬ 
tion  of  the  second  kind . 

By  adjusting  the  phases  in  (17. 9)  one  may  obtain  two  's  which 
may  be  added  to  give  an  outgoing  progressive  wave.  The  behavior 
as  x  -►  -»  will  then  be  a  superposition  of  an  incoming  and  an  out¬ 
going  wave.  However,  one  may  also  modify  the  preceding  arguments 
in  order  to  treat  the  progressive  wave  problem  directly.  One 
specifies,  say,  an  incoming  wave  from  the  right,  writes 

,  f  )  *  C  *1  oc  /"vx  +  t ^  -4-  (<£)  Ayt  1 

(17.25) 

where  must  satisfy  the  radiation  condition,  and  tries  to  express 
the  corresponding  complex  potential  as  a  distribution  of  the  vortices 
(13.28)  since  they  already  satisfy  the  radiation  condition.  We  shall 
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not  dwell  on  the  detail#  except  to  remark  that  the  problem  lead# 

to  a  pair  of  coupled  Integral  equations  since  one  need#  a  dlatrl- 

* 

but Ion  not  only  of  (13.20)  as  It  stands,  but  also  of  the  vortices 
obtained  by  replacing  t  by  t  -  If  /2 4  .  This  method  could  have 

beer  applied,  for  example,  to  the  problem  of  the  vertical  barrier 
considered,  above . 

problems .  This  term  Is  generally  applied  to  water-wave' 
problems  In  which  the  obstruction  Is  a  horizontal  plane  of  finite 
or  semi-infinite  extent,  either  submerged  or  lying  on  the  surface. 
The  solution  lor  the  semi-infinite  dock  in  infinitely  deep  water  was 
given  by  Friedrichs  and  Lewy  [1948] ,  and  at  about  the  same  time 
the  same  problem  in  water  of  finite  depth  was  treated  by  A.  Heins 
|]l94'9]  who  also  allot.ed  a  restricted  type  of  three-dimensional  motion 
The  methods  were  quite  different.  Subsequently  Heins  [195^  and 
Greene  and  Heins  [1953]  extended  the  treatment  to  submerged  docks 
In  water  of  finite  and  Infinite  depth.  As  was  remarked  earlier, 
Rubin  [1054]  has  shown  the  existence  of  a  solution  for  the  finite 
dock  In  inflhiicely  deep  water.  Sparenberg  [  1957 J  has  deduced  an 
integral  equation  of  the  second  kind  for  this  problem. 

As  an  example,  consider  a  submerged  dock:  at  depth  b  and  ex¬ 
tending  from  x  »  -a  to  x  »  a.  The  integral  equation  (17 .20)  then 
becomes 

'Gy 
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where  K(x,  J )  .  K(x-  J  )  with 


(17 .26) 
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Without  actually  establishing  the  existence  of  a  solution  to  (17.26), 
Keldysh  and  Lavrent'ev  [  937  J  in  treating  the  flow  about  thin  hydro¬ 
foils  (see  section  20/5  )  propose  an  approximate  method  of  solution 
by  expanding  (x)  and  K(x)  in  a  series  in  X  «  a/2b: 

Hx,)  m  %0W  t  i]  l *)T  +-  •  •  > 

i  s  f  x  s 

^  t  ^  r 

and  determining  the  n(x)  recursively.  In  the  problem  treated  by 
them  the  total  vortielty  was  fixed  by  the  Kutta-Joukowski  condition, 
in  the  present  problem  the  corresponding  condition  is  still  to  be 
determined. 

If  the  dock  extends  from  -»  to  0,  one  may  modify  the  earlier 
arguments  so  as  to  apply  to  an  unbounded  body  and  derive  the  Integral 
0  qu£i  t*  1  on 

r(f)IC(x-  t)<Lj  =  -  l*A6  i  .  (17.28) 

An  integral  equation  of  this  form  is  known  as  a  Wiener-Hopf  integral 
equations  and  in  many  cases  can  be  solved  by  use  of  Fourier  trans¬ 
forms.  It  does  not  seem  possible  to  expound  the  method  briefl", 
so  we  refer  to  the  paper  of  Greene  a.id  Heins  [1953]  where  this 
problem  is  treated,  but  with  the  kei-nel  expressed  differently. 
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When  the  -«ml -infinite  dock  is  on  the  surface,  the  dock  may  be 
considered  as  a  limiting  case  of  a  beach  in  which  the  angle  between 
the  bottom  and  the  free  surface  ia  l8o°.  Although  waves  on  beaches 
are  discussed  in  the  next  section,  the  methods  which  allow  exten¬ 
sion  of  the  angle  to  180°  are  also  difficult  and  will  not  be  con¬ 
sidered  there.  They  may  be  found  in  Stoker's  Water  waves  [l957,$£4], 


17  /3 


•  Waves  on  l)1®  a  o  h1 o  a  * 

Let  the  fluid  at  i*est  be  contained  in  the  wedge  defined  by 
tan  <  l  ,  x  >  0,  f  >  0  , 


l.e.,  the  bottom  is  the  plane  x  sin  <x  +  y  cosoc  «»  0.  For  such  a 
body  of  fluid  one  may  look  for  periodic  waves  which  are  either 
standing  or  progressive  .  The  appropriate  mathematical  problem  for 
standing  wavep  is  to  find  a  velocity  potential 

(  J' |  ^  t  )  *  T  ( Ji,  C&3  +  T)  ggj 

satisfying 


AV  =  o 


^  0)  *  O  f 

^  03  f  z  o  l**  J(,  <MA~f  f  y  c*.  -  o 

//  ♦  f r  JO  ^ 

Jt’-fg1-*  oo 


This  p  oblem,  in  both  this  form  and  the  three-dimensional  form 
to  be  considered  in  £  18,  has  received  intensive  study  In  recent 
years  £e.g.,  Miohe,  1944;  Lewy,  1945;  Stoker,  1947;  Friedrichs,  1948; 
Isaacson,  1948,  1950;  Weinstein,  1949;  Peters,  1950,  1952;  Roseau, 


195?;  Lehmann,  1954 %  Brillaudt,  19&7]  ■  In  particular,  the  cited 
work  of  BrlllouHt  and  chapter  5  of  Stoker's  Water  waves  [ 1957J 
contain  a  general,  exposition  of  the  mathematical  theory*  We  ahull 
restrict  the  present  tieatment  to  simple  cases. 

Kirchhof f  [l879]  was  apparently  the  first  one  to  treat  the 
two-dimensional  case .  The  problem  was  taken  up  again  by  Hanson 
(1926]  who  considered  both  the  two  and  three-dimensional  cases.  Both 
these  authors  restricted  the  solution  to  be  bounded  everywhere.  Itils 
has  the  effect  of  excluding  a  physically  Important  class  of  solutions 
with  singularities  at  the  origin,,  One  ma>  see  this  easily  if  ■  90°, 
l.e,  when  there  la  a  vertical  cliff.  A  bounded  solution  Is  obviously 
<P(x,y)  ■  Ae  cos  1)  x,  u  -  i£^/g.  This  generates  a  standing 
wave  behaving  like  cos  v  x  at  x  «  ».  However,  if  we  wish  to  con¬ 
struct  a  solution  behaving,  say,  like  an  Incoming  wave  at  Infinity 
we  need  also  a  standing-wave  solution  behaving  like  sin  v  x  at  In¬ 
finity.  No  ou<~.i  solution  exists  which  is  bounded  everywhere.  How¬ 
ever,  as  we  shall  see,  it  Is  possible  to  construct  such  a  solution 
by  allowing  a  singularity  at  the  origin.  If  the  two  standing-wave 
solutions  are  used  to  construct  an  incoming  progressive  wave,  the 
consequent  loss  of  energy  associated  with  the  singularity  Is  some¬ 
times  Interpreted  physically  as  representing  loss  of  energy  in 
breaking  of  the  waves,  at  least  when  od  is  sufficiently  small  for 
this  to  happen.  There  is,  of  course,  no  a  priori  method  of  selecting 
the  mathematical  solution  best  representing  the  physical  phenomena. 

The  comparison  between  physical  waves  and  mathematical  solutions 
is  discussed  briefly  In  Stoker  [  1957#  pp.  69-77J » 
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Kirchhoff's  approach  to  the  solution  1®  Interesting  historically 
because  of  Its  similarity  to  the  method  used  later  by  Peters  [1950J 
and  Roseau  \l9Sl]  In  solving  the  problem.  HI a  reasoning  runs  as 
follows,  with  a  slight  change  In  notation.  Let  f(z)  ■  <p  +  1  T 
be  the  complex  potential.  Then 

2  ( x ,  ^ )  •»  +  +  jf"  (x  -  ir  )j 

y* cx,*j)  ®  f (*+  ‘  y)-  f(&-  i  y . 

The  free-aurface  condition  becomes 

'  {'*)]  -*[{(*)  +  {!*)] j 

But  then  ali.o 

{'(})].  (17-30) 

The  bottom  must  be  a  streamline.  Hence 

{(rtif)~  4(  re  Lf)  "  const.; 
we  mty  take  this  constant  as  0.  Prom  this 

{ (})-{(y'l*r)  (17.31) 

Hence 

4-  [  Hi  )“  f  0  *Ul r)  ]■ " -  «  [f  (}  )  7-  f  <}  c  Uf)  ].  (17.32) 

This  differential  difference  equation  muat  hold  for  all  z  for  which 
_  l  9  y~ 

f(z)  and  f  (ze  "  )  are  both  defined,  namely  for 

-T  *  }  <  f- 

Kirchhoff's  formal  arguments  need  to  be  supported  in  terms  of 
analytic  continuation  by  the  reflection  principle,  but  the  essential 
idea  is  the  same  as  that  used  more  recently  [_cf.,  e.g.,  Lehmann, 

1954,  f  3,  or  Peters,  1950,  §  3 j . 
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Klrchhoff  proceeds  to  aolve  this  equation  In  the  special  case 
»  m  T/n,  m  and  ri  relatively  prime  integers,  by  assuming 
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Substitution  In  (17.12)  gives 

4  it, 

l+t)  •  Ak„,(fr  A-  f)  ,  U*ot:,n-l,  (17.34) 

with  A_1  —  An-1.  Multiplying  all  equations  together  and  remembering 
.  ITi  •“X 

that  l»p,  .  ..,p  are  all  nth  roots  of  unity,  one  finds 

a r-/. 


which  can  hold  only  if  n  is  even,  say  n  -  2q  (hence  m  is  odd).  With 
JL 

A  ...  -1  »  fS  *  ,  the  cbove  equations  determine  Successively 

A.,  ,  A  ,  in  terms  of  A  .  and  A„  -  . . .  -  A  ,  ■«  0: 
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Aq  is  still  an  arbitrary  complex  constant.  The  differential-differ¬ 
ence  equation  is  a  necessary  condition  for  f(z),  but  not  sufficient 
to  ensure  that  all  boundary  conaitiona  are  satisfied,  if  one  sub¬ 
stitutes  the  above  expression  for  f(z)  in  (17.31),  one  finds  after 
some  computation  that  one  must  take 

A.s.e^HA  ( 


(17.37) 
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where  BQ  is  pure  imaginary  (say  IB^}  if  both  ijj-(m+l)  and  q  are  even 


tnd  otherwise  is  real,  with  this  choice  of  A  one  has 
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(17.38) 

As  Klrchhoff  points  out,  the  solution  is  physically  acceptable 
for  the  problem  at  hand  only  if  m  *  l;  otherwise,  does  not 
remain  bounded  as  x  *►  +<»,  If  m  -  1,  then  for  y  »  0,  the  dominant 
term  as  x  •*  «  is  given  by 

f(x)  ~  £>•  (**  iv  x  -  i  W  ) 

or 

„  ,  o,  (17.39) 

fU ,o)  ~  bcw  (~oxt  *) . 

Here  are  several  easily  computable  special  cases  of  (17.36): 
<f“  90° (  m  -  1,  q  -  1,  b  -  -1): 

f  (} )  -  Bv  eVl  -  60  (cot  vx)  , 

f-  45%  -  1,  q  -  2,  f>  -  -1); 
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numerical  computations  for  ^  (x»y)  for  f  «*  6°  (q  ■  15)  as  we1!  a® 
for  the  above  cases  were  carried  out  by  Stoker  [l9*7]  and  are 
presented  graphically  In  his  paper, 

Kirchhoff's  aolutlon  lo  limited  to  the  special  choice  of  angle 
noted  above  and  furthermore  presents  only  solutions  which  are  bounded 
at  the  origin.  The  solution  of  the  differential-difference  equation 
(17.32)  for  arbitrary  (  ,  0  <  <f  <  I ,  has  been  given  by  both 
Peters  [1950] ,  Isaacson  [1950],  and  Roseau  [^1952,  chap.  vj.  All 
use  Laplace  transforms.  However,  the  method  cannot  be  expounded 
briefly  and  ks  refer  to  either  the  original  papers  or  Stoker's 
Water  waves  for  the  details. 

T  >e  special  cane  f  •  3T  /2q  can  be  treated  fairly  simply  by 
the  reduction  method  used  In  the  problem  of  the  vertical  barrier. 


Prom  ( 17  •  32 )  we  have 
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The  free  surface  condlt1  on  |]cf .  (11.7)]  implies 
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Thin  last  equation  can  alto  be  written 
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If  the  numbers  &k  and  are  real,  (17.43)  and  (17.44)  imply 
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Me  wish  to  find  numbers  f j  and  ^a£  j  such  that 
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Comparing  coefficients  of  lerivatives  of  the  aame  order,  one  finda 
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These  relations  will  be  satisfied  if  one  takes  s  ■  q  -  1  (for  ft  ^—1) 
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(17.48) 


We  note  that  a  ^  *  ~v  With  this  choice  of  the 


coefficients  (\J  >  define 
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where  the  last  equation  follows  from  (17.42)  and  where 


|r'  4 

p(A)  «Z  \  t 

(5_  -  O 


(17.50) 


Prom  the  assumptions  originally  made  concerning  f(z)  and  from 
the  method  of  selecting  the  { >  it  follows  that  g(z)  Is  regular 


everywhere  In  the  wedge 


-  2  i  £  <9-0 


except  possibly  at  the  origin,  that 


Im  J  =  O  -J-o-r  ^  *  jo  >  o  <jW  -  r  e 


-2  i/- 


finally,  from  the  last  of  equations  (17. 24 )J  that 


g(/4z)  -  -g(z). 

Since  f(z)  is  assumed  bounded  as  x  +  *,  this  is  true  also  of  g(-). 
These  various  conditions  imply  that  g(z)  must  have  the  form 


oC  ^ 


i  <vj  rt  ax 


(17.51) 


We  have  thus  shown  that  f(a)  satisfies  the  differential  equation 
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(17.52) 


Prom  the  definition  of  P(  )  It  follows  that 

PfA^A  +  tv)  s 

Since  the  coefficient®  In  P{  are  real,  A  1*  '•  roo*-  of  P(  A  )  -  0 

If  A  is  a  root.  Furthermore,  from  the  Identity  al>  ■«  also  A 

Is  a  root  providing  /*>  X  •/  li>  .  since  A  “  1  ^  ls  cbvio\s  root 

of  the  left  hand  member,  -i/5v  la  alao  a  root  an  ..ence  -1/3>2u>  , 

3  a 

“i/5  ^  .  Since  ^j>H  -  -1 ,  no  new  roota  are  added  by  going 

further  than  -1/3  q"1v  ,  and  since  i/£"kV  — lyfcq”k-w  ,  a  complete 
set  of  roots  of  p(  A  )(  A  +  iv)  Is 

i  S'1 

-  iv,  -  i, fix >,-<./!>  v,  -  t,  v 


Thus  the  solution  of  the  homogeneous  equation  can  be  expressed  In 

the  form 


SoA<c 


(17.53) 


Thl3  la,  of  course,  exactly  the  form  of  Kirchnoff'3  solution  cf 
(17.36).  Since  we  have  already  determined  the  necessary  form  of 
the  In  order  to  satisfy  the  boundary  condition  on  the  bottom, 
we  need  not  pursue  further  the  solution  of  the  homogeneous  equation. 

The  solution  of  the  nonhomogene oua  equation  Is  straightforward. 
However,  just  as  fox'  the  hotnogeneouf.  equation,  one  must  take  care 


1 


ii»^  ,  m  im<-  mi  t 


. mm . . . 


to  satisfy  the  boundary  condition  on  the  bottom,  l.e. 

,  j  .  o.  ».  detailed  con.ld.ratlon.  »ay  be 
found  in  the  several  cited  papers,*  Brilloult  [1957]  treats  the  mat¬ 
ter  thoroughly.  If  one  considers  (17.52)  with  the  right-hand  side 
replaced  by  only  one  of  its  summands,  say  b  *”(2n+l)q^ 

Jil 

complete  solution  can  be  put  In  the  following  form,  as  shown  by 
BrillouStt 


(17.54) 


where  cn  is  an  arbitrary  real  constant,  Bq  oi'  (17.37)  has  been  set 

equal  to  1  ,  and  where  r  k  indicates  that  the  integral  is  to  be 

carried  out  over  each  of  the  paths  T  *  and  ['  ~  shown  in  Figure  17 


* 
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Figure  17. 


However,  one  my  obtain  a  variety  of  other  forms  for  the  solution. 
An  asymptotic  expression  as  x  •+•  »  and  for  y  -  0  is  given  by 


fU)~  few* i  c 


(17.55) 


*9~t 


~  cn  cm(vx  +  r  -~r-  )  t  ^('iy  f  3  *Z~). 


In  the  neighborhood  of  *  »  0,„  f(x)  behaves  like  log  z  for  n  *  0 
and  like  for  n  >  0. 

It  ia  not  clear  physically  what  type  of  singularity  at  z  «  0 

most  nearly  describes  the  behavior  of  real  waves.  However,  most 

writer®  have  restricted  their  treatment  to  the  weakest  possible 

singularity,  :!*•»,,  the  logarithmic  one. 

Prom,  the  asymptotic  expansion  as  r.  *►  «®  one  sees  that  it  ia 

now  possible  to  construct  an  incoming  progressive  wave  by  proper 

choice  of  the  constant®  c  and  b  .  Thu®,  if  we  select 

n  n 

■n-2  '  ,  _  . 

cv  *  61  'Dec!  (4  t  f  T  )  t  4*.  -  -  (—  I  )  JT  (l  *  g  - 1)  !  fp  &  \+  0 

then  the  resulting  solution  will  behave  like 
(X  Ctn  (  v  jt  f.  4 1  i  X) 


an  x  *  for  y  -  0.  In  connection  with  (17.55)  and  the  selection 
of  ton  just  made,  it  is  apparent  that  the  formulas  (17.54)  and  (17.55) 
will  be  more  directly  connected  with  parameters  with  a  simple 
physical  interpretation  if  we  replace  bn  by 


For  n  -  0  companion  singular  solutions  to  the  regular  solutions 
(17 .40)  and  (17*41)  are  not  difficult  to  write  out: 
m  90° (q  m  l,  n  *  0): 


IMMMMNMI  'M* 


45° (q  «  2j  n  -  0)  j 


V'iJK'  ■*! 

f(*,y)  -  —  e"1’  [({+  S  Ci(vx)cm(v Jt +  *)+£& £ 


(1757) 


Further  formulas  lor  $  —  ho  1» Jl  I  -1  <5  nwiy  be  found 

In  Brllloult  (1957,  p.  93ff , ] . 

1 8 .  _Th ree - dimej Lionel  progressive  and  a«~ ending  wave s_ i n 
unbounded  regions  with,  fixed  boundaries , 

The  general  remarks  at  the  beginning  of  section  17 
apply  here  also.  Although  most  of  the  solvable  problems 
In  the  present  category  are  such  that  they  can  be  reduced 
to  two-dimensional  ones  (however,  see  the  end  of  section  19ft) , 
the  methods  of  complex- function  theory  are  no  longer  applicable 
to  the  same  extent.  The  division  of  topics  Is  the  sane  as 
In  the  last  section,,  namely,  diffraction  of  waves  by  obstacles 


and  waves  on  beaches., 
18oc,  Diffraction  of 


_w«  ter. _  waves., 


In  a  horizontally  unbounded  ocean,  of  uniform  depth  k 
assume  that  an  incoming  wave  is  s peel fled  by 


lPl[  •*'  >  J,  0  "  --~M  mf *  k)  c&i  ( mx  +  +  **) 


(18.1) 


and  that  it  is  scattered  by  one  or  more  obstacles  in  the 
water.  Me  wish  to  find  the  velocity  potential  for  the 
motion  of  the  water  In  the  form 


]•*)-  4>  [  +  fs  . 


(18.2) 


condition  xf  the  body  lit  of  bounded  extent. 

A®  usual,  we  may  write  (b  in  the  form 

where  Cf  must  be  a  potential,  function  satisfying 


(jtfO,  JJ )  ~  v  «f  (x,o,  j  ) 


- 


VJv  *  fJn  +  'fSri  *0  on  tht  obstacles t 


(18.4) 


(aav,  \  R  ( I; 
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^  )  -  0  ,  Jr  ^  •  d(t)  as  R  -  oo 


Cfe^ra.1  obstructions.  Consider  a  single  submerged 
obstacle  bounded  by  the  surface  5  .  We  shall  try  to 

express  the  scattered  wave  <p$  *  &t  %  €  by  a  distri 

button  of  sources  over  5  .  However,  in  order  to  satisfy 

the  various  boundary  conditions,  we  take  sources  in  the 
complex  form  (13,18)  (or,  la  the  case  of  infinite  depth, 
in  the  form  (13.17'n: 


(18.5) 


wher*.  we  have  written  Q  •  C,  4  L  for  the  complex  for., 

of  (13.18),  The  boundary  condition  on  the  body  now  becomes 


0  _  2ft  +  M.„£fL 
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2 f  -  J)  •  (Id. 6) 


Since  J  dn  la  a  known  function,  this  Is  a  Fredholm. 

iiwft:irbrainiltthai*rL^i  "Wh,- 

.  ill! 


. . . . . . . . . . mln,iliilillin*iiut,l'>!iini||||(| 


ip-  M  --  T^IW-  -nn-^r-  ,,,.1-,,..  '  ?  } 

WlflWiPliPIIII  ir»iHa  1  tiJP*itW1!I!!FP^  |p 


'  ■» . ■■**■■ 

Wwplnwlmill 


1  1  '■' 


Integral  equation  of  the  second  kind  for 

if1** 

(We  note  in  passing  that  If  the  motion  of  the  surface  o  had 
been  prescribed  to  be  »  then  the  #ame  Integral 

Mw» 

equatl  m  for  0  would  ha 'e  been  obtained.) 

This  equation  has  been  considered  by  Kochln  [1940]  in 
the  case  of  Infinite  depth,  and  he  proves  that  a  solution 
exists  If  v*  <Sl/y  is  large  enough.  Iterative  pro¬ 

cedures  for  computing  if  follow  from  the  theory. 

Haskind  [1946]  hat  extended  the  argument  to  finite  depth. 

John  [1950]  has  treated  both  the  uniqueness  and  exist¬ 
ence  problem  in  great  detail  and  has  shown  that  a  unique 
solution  exists  for  a  body  whose  surface  Intersects  the  free 
surface  perpendicularly  and  which  can  be  represented  as  a 
single-valued  function  over  the  area  enclosed  In  the  inter¬ 
section.  His  result  holds  for  all  values  of  m  (or  V  if 
the  depth  Is  infinite) .  He  also  reduces  the  existence 
problem  to  solution  of  an  integral  equatior . 

Vertical  cylinders .  When  the  obstacle  or  obstacles 
are  vertical  cylinders  extending  from  above  the  free  surface 
to  the  bottom,  it  Is  possible  to  reduce  the  problem  to  one 
in  diffraction  of  sound  waves  for  wf ich  many  special 
solutions  are  known  (see,  e.g.,  Havelock  [1940]).  In  this 
case  we  may  separate  the  ^  variable  in  the  manner  shown  in 
section  13 o( : 

(f(x‘V  a«.7) 


where 


and 

f  rnV/  -0.  (18.8) 

Here  kn  must  be  the  came  as  in  (18.1)  since  the  frequency 
is  fixed  by  the  incoming  wave.  (x,^)  must  now 

satisfy  (18.8)  and  the  second  two  conditions  of  (18.4).  This 
is  exactly  the  same  mathematical  problem  encountered  in  the 
diffraction  of  sound  waves  by  a  cylindrical  body  (in  that 
case  the  air  pri.isure  replaces  ).  Thus,  any  solutions 

known  for  sound  diffraction  by  cylinders  may  be  taken  over 
immediately  for  water-wave  diffraction.  For  example,  if  the 
obstacle  is  a  vertical  circular  post  of  radius  ft,  the  velocity 
potential  of  the  scattered  wave  is  given  by  [see  Morse, 
Vibration  and  sound,  2d.  ed, ,  New  York,  1948.  pp.  347 f£,,  449] 

%(*,  9,  eM.  W,r  c«Q  ^(mR), 

1  (18.9) 


where  (JTk.  "  {**'&*) j  Vk  (m  an8 

£<,  ■  i  ,  t  K  =  2  f*r  n  £  I  . 

Various  approximations  for  large  and  small  values  of  meu 
are  known.  The  maximum  wave  amplitude  at  any  point  is 
given  by  ||f| 

The  diffraction  of  water  waves  by  a  vertical  half-plane 
may  also  be  treated  by  transferring  known  solutions  due  to 
Sommerfeld  for  sound  and  electromagnetic  waves  to  the  present 
context.  This  has  been  done  by  Hasklnd  [1948]  for  normal 
incidence  and  by  Penney  and  Price  [ 1952a]  for  both  normal 
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and  obliqua  inci dance .  Peter*  end  Stoker  {1954)  (tee  also 
Stoker  [1956]  and  [1957,  pp.  109-133] )  have  also  «olved  this 
problem  by  a  new  and  rather  easy  method,  following  an 
investigation  of  boundary  conditions  which  will  ensure 
uniqueness.  The  solution  has  an  obvious  application  in 
predicting  the  effect  of  breakwaters.  Let  the  breakwater 
be  the  half-plane  j  -  0  ,  X  >0  and  the  incoming 
wave  be  given  by 

*1  m  A  CM  (/y^  JC  on  o<  +  J  6  Of 

“  A  C6*  (aw  K  cel  (0  -  f  <o  , 

where  o<  is  the  angle  between  —OjC  and  the  direction  of 
propagation,  measured  clockwise.  Then  the  solution  given 
by  Peters  and  Stoker  is 

£2  cal  m(vL)[J.<K)*rf i  ]■  (18-10) 

The  result  can  also  be  expressed  by  means  of  Integrals.  In 
the  case  of  normal  incidence  these  reduce  to  Fresnel  integrals 
for  which  tables  exist.  Graphical  representations  of  the 
behavior  of  the  wave  amplitudes  may  be  found  in  Penney  and 
Price  [1952a]. 

Penney  and  Price  also  apply  this  analysis  to  an  approxi¬ 
mate  treatment  of  diffraction  by  a  breakwater  of  finite 
length  and  through  a  gap.  The  results  are  presumably 
applicable  if  the  wave  length  is  small  compared  to  the 
length  of  the  breakwater  or  the  gap. 


In  two  physical  situations  tha  dependence  upon  J  may  be 
precipitated  out,  leaving  a  two-dimensional  problem  which  in 
many  cases  can  be  solved  by  methods  analogous  to  those  used 
for  the  two-dimensional  problems  of  section  17. 

Let  the  obstruction  be  an  infinitely  long  horizontal 
cylinder  parallel  to  Qj  .  This  might  be,  for  example,  a 
semi-infinite  dock  or  submerged  plane  barrier,  say  ^  *  *  if,  M < ot 
a  finite  horizontal  barrier,  say  -  -.p  ,  |x|<  a  ,  a  vertical 
barrier,  xmO,  C  Y  d  o  ,  e  beach,  ^  »—  jc  h ,  etc. 

Let  en  Incoming  plane  wave  at  infinity  propagate  at  an 
angle  d  to  the  X  axis: 


7  5 1 0  *  ^  <*•)  +*^3.  (i8.il) 


Although  one  will  not  expect  the  velocity  potential  <f> 
to  he  periodic  in  X  ,  it  teems  reasonable  to  assume 
that  it  will  be  periodic  in  J  .  In  feet,  we  shall  assume 
that 


^,<7,3,0  - 
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(18.12) 


where  y)  must  now  satisfy,  with  4z 

fxx+ 


/vk. 


/ 


(18.13) 


and  the  usual  conditions  on  the  free  surface  and  rigid 
boundaries . 
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He  should  have  come  to  the  seme  conclusion  If  we  had 
a  e  e  uaie  d  an  incoming  wave  at  Infinity  of  the  form 

A  — CfllXcd  ®h03T t C 3T© l!J di  VP49IVC  (  liOt^P  fchttt  WC  -fll 0' JJUttmC  me#  )  • 

That  la,  we  shall  now  look  for  a  solution  in  the  form 

s  f  CUt f  (18.15) 

satisfying  equation  (18.13)  and  the  conditions  on  the  free 
surface  and  rigid  boundaries.  Thus,  a  solution  for  one  of 
these  cases  carries  over  easily  to  the  other. 

The  problem  is  thus  reduced  to  one  almost  identical 
with  that  of  section  17,  with  the  exception  that  the  two- 
dimensional  Laplaclan  is  replaced  by  (18.13).  Many  of  the 
same  methods  may  be  carried  over,  e.g„,  the  reduction  method 
and  the  integral-equation  method.  Haskind  [1953]  has  con¬ 
sidered  some  general  aspects  of  the  problem  which  will  be 
outlined  below,  has  derived  the  source  solution  of 
(18. 13),  and  has  treated  the  diffraction  about  a  vertical 
barrier  (an  analogue  of  the  problem  treated  in  17  ol  )  and  a 
finite  dock,  all  is  infinitely  deep  water.  MacCamy  [ 1957 1  has 
derived  a  source  solution  of  (18.13)  and  treated  the  finite 
dock  problem  in  water  of  finite  depth.  Heins  (1948,  1950, 

1953]  has  given  source  solutions  of  (18.13)  for  finite 
depth  and  formulated  and  solved Wiener-Hopf  Integral  equa¬ 
tions  for  semi -infinite  docks  and  submerged  horleontal 
barriers.  Greene  and  Heins  [1953]  treat  the  submerged 
barrier  In  water  of  infinite  depth.  The  literature  for 

i  '  Jill  i  , _  u  ,  (  ■  •  ,  |  i  i  (  i'|i«  r 
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beaches  will  be  given  In  section  !8£  . 

Suppose  the  fluid  infinitely  deep  and  let  a  cross- 
section  of  the  obstacle  have  contour  C,  We  wish  then  to 
find  a  solution  ?'7jr(4j)  *  V*,  +  t  ^  of  (18.13)  such 
that 


V  tftXfO)  -  o  on  ihc  free  Surface  f  (18.16) 

f~  £-2  evVe‘4-  e«  x~+-. 


/  t  i. 

where  ^  V  ,  Haskind  [1953]  applies  the  reduction 

method  In  the  following  manner  (we  follow  hla  presentation 
closely).  Introduce  the  function  •{( X ,  by 


3f  3^  * 


^  O  ^ 


Then  f  also  satisfies  (18.13)  and 

{^(XjQ)  -  0  on  fhf  /ret  surface^ 


(18.17) 


(18.18) 


Consequently,  f  may  be  extended  Into  the  upper  half-plane 
by  defining  *f  6* ,  and  f  now  satisfies 

(18.13)  in  the  whole  plant  outside  the  contour  C  and  its 
mirror  image  C  .  Moreover,  lfl~*0  as 
Assuming  that  f  is  known,  one  must  now  try  to  reconstruct 
f  from  f  tr,  such  .  «y  th.t  condition.  (18.16)  sre 


satisfied.  In  order  to  do  thl»,  Hasklrd  differentiates 
(18.17)  with  respect  to  ^  ,  subtracts  from 

■fjlX  ^  'f'lj  Aj  ~  ^  M  0 f 

and  after  sane  eat/  manipulation  obtains 

C (f +  »{).  { 


(18.19) 


Treating  this  as  a  differential  equation  for  y^-  -f  ,  he 
finds  the  following  solution  for  ^ : 

-i  ktx  f*  Z  /  ,  \  ,  i  7  49  v*f  -L^c,X  (18.20 

-e  _L e ’  U»  +  5  t)dil  +  7i*  ~e  / 

the  Integrals  being  taken  along  half-lines  parallel  to  the 
x-axis  and  below  C.  One  may  verify  without  great  difficulty 
that  satisfies  (18.13).  The  asymptotic  form  of  f  as 
jc  —*  ±  qno  may  be  written  down  immediately,  and  gives 


*±xi;itc<x 

C  ^  /  (18.21) 


the  path  of  integration  being  a  line  below  the  body. 

Consider  now  the  region  D  bounded  externally  by  this  line 
and  a  large  semicircle  containing  C  +  C  and  internally 
by  C  +  C  .  Application  of  Green's  Theorem  to  -f  and  \  = 
exp  (-v^  •+■  i  4;  iX  )  shows  that 
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a  definitely  related  Government  procurement  operation,  the  U.S. 
Government  thereby  incurs  no  responsibility,  nor  any  obligation 
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furnished,  or  in  tiny  way  supplied  the  said  drawings,  specification 
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in  any  manner  licensing  the  holder  or  any  other  person  or  corp 
tion,  or  conveying  any  rights  or  permission  to  manufacture,  use  c 
sell  any  patented  invention  that  may  in  any  way  be  related  ther 


Hence,  the  asymptotic  conditions  are  satisfied  and,  more ove r 


l[fx- 


(13.2 


By  a  similar  application  of  Green’s  Theorem  Haskind  shows 
that  one  may  also  write 

r-  i**pfcv\**  £p-  <18.. 

0<5  ' 

where  the  +  sign  is  used  for  points  to  the  right  of  C 
and  the  -  sign  for  points  to  the  left.  It  is  easy  to 
verify  directly  that  "atisfies  (18.13)  and  (18.17); 
however,  (18.20)  allows  one  to  investigate  the  asymptotic 
behavior  more  simply.  If  has  no  singularities,  then 
(17.3)  must  also  hold  here,  i .e . ,  ^  +££)*  +  2A3  mO 

This  result  may  be  used  to  find  the  source  solutions 
giving  outgoing  waves  at  ± -o  .  For  equation  (18.13)  the 
singular  solutions  for  the  whole  plane  are  known  to  be  the 
Bessel  functions  cr )  >  where  *  ,0c  -&-)1  +  (y  -  t 

To  find  the  solution  corresponding  to  (13.22),  one  assumes 
it  may  be  expressed  as 

£(*,*1  ;^,4)  -  <fo  +Ko(4c  r)-!  (4tr,), 

with  r/  *pc  cif  f  ( ^  i  4>V /  ,  where  i  has  no  singulari¬ 

ties  for  ^  <  c  ,  Then  foym  9*^  -  v  'fo  may  be  extended 
as  a  regular  solution  of  (18.13)  to  the  whole  plane.  Also, 

{<»/* ,o)  *  2  ^  :  , j  * 

One  may  then  show  that  this  relation  holds  for  all  <  0  ; 


Substitution,  in  (18*23)  with  A  «  o  #nd  direct  computation 

JL 

of  B  from,  (10,22)  by  talking  C  as  a  small  circle  about 
the  singularity  gives 

/«  V  1;  .  .  *  V4I  r*  .  I  J 


*"i N, 

C*  *K0(kr)  +  kt(kr,)+  2ve,4t  f  t  K0fkr,)du 

043*  <1 


(18.24) 


2Tir  e"V4)* 


For  Raskin**  a  application  of  this  method  to  the 
diffraction  about  a  vertical  and  a  horizontal  barrier  we 
refer  to,  the  original  paper.  Force  and  moment  are  obtained 
in  terms  of  Kathlau  functions.  Fox  the  horizontal  barrier 
in  water  of  finite  depth  we  refer  to  MacCamy’s  paper  11957) 
where  a  formula  analogous  to  (18.24)  is  derived. 

18 />  * 

Much  of  the  Immediately  preceding  discussion  of  diffrac¬ 
tion  of  plane  waves  approaching  at  an  angle  or  of  short-created 
waves  approaching  normally  applies  also  to  this  case.  One  is 
led  to  the  following  boundary-value  problem  for  +  if2 


15  “£*  +  't*n  -  C'f-O  ,  4l< 


2)  t,,)  -  V  </(#.,<>)- O, 

))  %  «S*r:+  f*t  mfO  J*r  1** *m.f 

k :  i  f.  ■■ }  il  j  Ji !  i 


(18.25) 
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Many  of  the  authors  cited  in  section  17  considered 
this  problem  along  with  the  two-dimensional  one.  In 
particular,  we  refer  to  Hansct  [1926],  Miche  [1944], 

Stoker  [1947],  Weinstein  [1949],  Roseau  [1952],  and 
Peters  [1952].  Both  Peters  and  Roseau  solve  the  problem 
for  arbitrary  angle  f  ,  7  (thus  including  the 

semi-infinite  dock  problem  treated  differently  by  Heins 
[ 1948] .  The  use  of  fcha  xr 48' (Jiac t f on  RM8t*hod  1. Icni cmtm®  hanra-ii 
as  in  the  two- diman •ional  case,  to  angles  fMP  i 
We  shall  illustrate  the  procedure  briefly  for  J" m  Tfa  a nd 
following  essentially  Weinstein's  [1949]  treatment  (see 
also  Bril  lout!  t  [1957,  chaps.  I,  II]}. 

Since  the  boundary  condition  on  the  free  surface  and 
bottom  is  the  some  in  the  two-  and  three-dimensional  cases, 
we  may  make  use  of  the  auxiliary  function  g  constructed  in 
(17.49)  by  using  only  the  real  part  of  the  complex  potential. 
Thus,  for  ^  one  finds  from  (17.48)  that 

3- 1  =  Q.b/V  •  Hence,  from  (17.50) 


p  (X)  -  a0{  l  +  ^h) , 


7 

Thus,  the  boundary  conditions  2)  and  3)  of  (18.25)  imply 
that 

V  r~  *0  On  yO,  x  >0  (18.26) 

^  Jt  *■  ^  ^  £  Q  • 

We  recall  that  the  definition  of  cflx,  y)  has  been 

extended  from  the  original  wedge  by  reflection  in  the  bottom. 
One  must  now  find  a  function  'b(Xi^)  satisfying  equation  1) 
of  (18.25)  and  the  boundary  conditions  (18.26)  and  which  is 
regular  everywhere  in  the  extended  wedge,  Op®  >/ ■%.  J, 
except  possibly  at  the  origin,  bounded  as  «*  +  y  *“*"  00  / 
and  symmetric  about  the  line  =  -  X  ,  it  is  I'nom  that 
the  general  solution  of  this  problem  is  given  by 

)  -  (^2  +  (^j  (*1  *])  «Z  Aw  ^ ifznH)(kr)s^2(Zn+p0X  18.27) 

A  similar  analysis  for  waves  approaching  a  vertical  cliff 
(  Jf  —  ^  IT  )  lead®  to 

^  (18.28) 

Let  us  take  the  weakest  possible  singularity  In  each 
case,  iL.e. ,  for  the  90°  cliff  and  K2  for  the  45°  beach. 
Consider  first  the  vertical  cliff.  Teking  AO  O  OUtl  t  0  te  tlA 
relation  k' (u,)—k,U)  ,  «  kwa 

'  v)  /(*•  V  -  -  jt  TZ &(*  o'. 


t#  .limiii 


2,  o  *2* 


Ws  may  than  Identify  -  j4tj  with  f  and  from 

MO  no\  r»  •$*  — 


(18*23),  with  £>  m  o  i  we  have 

10  A .  .  _  ✓  .  v  .  -U  fi  _  ,  , 


^ - -  £  K.(Lr)-A.2tV1  fl-'l 


whare  /!„  must  atlll  be  determined  ao  that  f^(o,  m  Of^  <0 
In  computing  fx  v.  £-*  0  ,  one  must  remember  that 

k0(t* )  ~  A*  **  U  -*  o  .  Hence  i  one  If  In  els 


'{*(0.*})  -  -  ^  v*  -  ^  f  -  l--?4'  tl 

®  fe  4  ■*  0  J(  "*  C  ,  7  ^  • 


-4?  vie)  -  t  e/T 

*■  6 


Setting  this  equal  to  zero,  one  find® 


A*  _  ^  ^  e  4 , 

fc  JTV 


Substituting  above  and  separating  the  real  and  Imaginary 
parts  of  cf„cf,  +  i(fz  ^  obtain  fln  everywhere 

regular  solution  ^  a-nd  eo  solution  fi  with  a 
singularity  at  the  origin  and  90°  out  of  phase  at  X  m  oo 


yf  (*jCi  ^  ) 3,1  e  cm  4r,  x , 


%  *■  “  jt<{0({cr)+  V  <°V^  j  ^  ^  (k  f^+f  j  ] 


(18.29) 


+  4*  ^s^4(x. 


The  cot responding  equation  for  (18.27)  can  be  written 


in  the  form 


I . 

1 . 

I 


r 

i 

I  a 

II 

|ii 

I 


II 


illHiHMHIHi  IHmMn 


hi ^  nmnn  i  i 


2  o  3 


On#  can  find  it*  Integration  discussed  in  Roseau  [1932, 
eh,  IV].  A  solution  for  the  next  simplest  case,  f  *  *0* 
does  not  seem  to  have  been  published.  For  f  *  the 

regular  solution  9*i  ,  and  singular  solution  fi  as  given 
by  Roseau  are: 


A 


^t’^kcmkfX  -  v  »*«■  4c,  Jt]  f  t  ’  X/4(re-*4,^,  t  ******  it,  y]  J  t 


*Z  *  A ije  ^VCesi, Jt--4l‘uw.4,jtJ+C  X^C#9^i^+ t 
+4?  j]*L  f ~  4,/4  UV^xy  r  u e  eZ*  (4t  I  jfVy2- ) 

-V>  4f  . , 

+  vt"*  ft  '"I n. (u  fcvy ) ^  - xV*'? j d  n l'"7 (iM fj.f 

Ja 


(18.31) 


A(J  ~'o*> 

In  order  to  satisfy  condition  4)  of  (18.25)  one  must:  take 

A  .  ~A±  ,  kdJ±  ,  a  m£±.  v> »  «,  4  ,  . 

/Ml*  y  ,  ...  1  hi  j  7-r 
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4*  ^  4  4V 

Edge  Waves »  In  the  Investigation  of  diffraction  of 

waves  on  hori cental  cylindrical  obstacles  and  of  waves  on 

1**  • 

beaches,  It  was  specifically  assumed  that  4  ./Vvv  m 

This  was  automatically  fulfilled  for  plane  waves  approaching 
at  an  angle,  but  needed  to  be  assumed  for  short-crested 
waves.  For  the  short-crested  waves  there  also  exist  standing- 
wave  (solutions  which  can  be  exhibited  in  certain  cases  for 


4.  ^  )  awZ 


Such  solutions  were  apparently  first 


noticed  by  Stokes  [1846,  p.  7-1880,  p.  167]  in  connection 
with  the  propagation  of  waves  in  a  canal  of  non-rectangular 
erftfcjfaJMfCloa.  Cnrtptn  peculiarities  of  these  solutions 
been  pointed  out  by  UrseU  [1951,  1952]. 


rnmmmm 


*2  <>4 


Consider  the  first  three  conditions  of  (18.25)  for 


V 

7  v  Then  one 

may  verify  directly  that 

*  <**fJ 

7  (x,  )  »  e  ^ 

is  a  solution.  This  gives  a  velocity  potential  for  stand¬ 


ing  waves: 


(18.32) 


where 


OA*vf  •  */q 


The  wave  amplitude  is  bounded  at  the  origin  and  drops 
off  very  quickly  as  x  increases.  Clearly,  one  must  have 

f  C  jTt  .  Urs.-ll  has  pointed  out  other  interesting 

* 

aspects.  For  a  given  f  and  6  there  is  only  one  allowable 
4  ,  i.e.,  it  Is  -  discrete  point  of  the  spectrum.  In  the 

case  discussed  earlier  with  all  values  of  4 

between  0  .nd  -v  were  allowable.  In  addition,  the 
total  energy  per  unit  length  in  the  J  direction  is 
finite  for  the  Stokes  edge  wave. 

From  (18.29)  one  may  construct  a  progressive  wave 
moving  in  the  direction  O ^  with  velocity. 


There  is  evidence  that  such  waves  have  beer,  observed  in 
nature  [cf.  Hunk,  Snodgrass  and  Carrier,  1956;  Dunn  and 
Ewing,  1956J . 


. . .  . . . . . . 


x  or 


Urscll  [1952]  has  shown  that  (18.32)  Is  only  the 
first  in  a  sequence  of  solutions  of  this  nature  for  a 
sloping  beach.  He  shows,  in  fact,  that  the  following 

velocity  potential  also  satisfies  the  condition: 

,  ,  f-Ux  c^r-  *u„r]  ?  a  r-kf  c*(2^-Qf4yV~(Z*»-0r 


‘  (18.33) 

^L[xor>  (ZMv+i)f-  y **w j |  J  «.  £)e^Ui  -hr) 


where 


r-  I 


ft*,  (i-v-r  f  1  jj 

tow  (k.  +  r)  f 


f  &l*  <^4«.  (?-*»■*■!)  f 


It  follows  from  the  last  condition  that  one  must  have 

(Z*v  +  \'Jf‘£  \  or  ^  ^  ^  +  2  ; 

where  n.-O  will  be  taken  to  indicate  the  Stokes  edge  wave. 
Thus,  for  fixed  wave  number  iz  ,  the  above  formula  gives 
one  frequency  &  if^T>f>gJT  ,  two  if 

1  if  >  f  >  ^  IT  ,  etc  An  experiment  carried  out 

by  Ursell  confirms  the  existence  of  these  other  inodes  of 
motion.  The  solution  (18.33)  for  ^  =  l>0  has  also 

been  given  by  Lamb  [ Hydro aynamics ,  p.  450,  eq.  (30)].  At 
the  critical  angles  }X/q  ,  */to  ,  etc.,  the  solution 
(18.33)  does  not  vanish  at  X  -*  ««=>. 

Keldysh  [1938]  has  stated  without  proof  that  for 

45°  the  Stokes  edge  wave  and  the  function  ‘fi  from 
(18.31)  constitute  a  complete  set  of  bounded  solutions  in 
the  sense  that  for  any  absolutely  integrable  function 
f  («*»  5)  J£  •  O  *  th®  follo**n8  Fourier -integral -like 

theorem  holds  lc£.  formula  (16.5)] : 


. . .  I! 


0  "  1  «/**  o 

■  {[k,e'x+  ^Jk.x-V^k,  a]  [1,^-1-  <c,c~ij  -V*~bj] 

+  24l  e.Aftf-'hUtD/nj  I  ■ 


It  is  possible  to  construct  other  types  of  edge  waves. 
First  we  rederive  the  Stokes  wave  from  the  third  formula 
in  (13.5)  with  iX m  O  .  A  surface  satisfying  ^  *  O  la 


de lined  by 


or 


} 


-  x  +  C  ,  v/V4l-yv  / 

-  n*t 

where  we  may  set  C  »  0  since  xc  '•Van /pro vide  essentially 

different  solutions  for  the  bottom.  This  is  just  Stokes' 
solution. 

One  may  expect  to  find  a  different  type  of  solution  by 
using  the  third  equation  of  (13.6)  with  &  •  O  .  Here 
the  corresponding  solution  is 


(18.34) 


where  again  ve  have  dropped  an  added  constant.  This 
describes  a  bottom  which  starts  as  a  sir  >ing  beach  and 
approaches,  as  JG  -*  oo  ,  a  flat  bottom  at  depth  C 
The  Initial  slope  of  the  beach  is  Nt-**©  ■ 

The  velocity  potential  describes  edge  waves  for  such  a 
c  on  t*  JLon « 


•nn{ifMnp 


1 


One  nay  proceed  to  the  seme  fashion  with  the  last 
formulas  of  (13,. 5)  and  (13.6).  They  turn  out  to  give 
Identical  hot tuns: 


-Coo  JdLJassiL^icii 

®  'tiw  ,Hn»i  (v 


(16.35) 


This  corresponds  to  edge  waves  along  an  overhanging  cliff 
in  water  of  finite  depth.  The  initial  backward  slop©  of 
the  cliff  is  4  Vj  +**3^  - 

A  particularly  interesting  sort  of  edge  wave,  although 
the  name  is  now  a  misnomer  since  there  is  no  edge,  her'  been 
discovered  by  Uraell  [1951].  He  has  shown  the  existence  of 
standing  waves  of  the  form 

(j£  |  M  )  00^.3  'Ll.  J  C#~.t  &  ?■ 

in  the  neighborhood  of  a  fixed  submerged  cylinder  of 
radius  <3.  If  i*dt  is  small  enough.  The  waves  are 
symmetric  about  the  vertical  plane  through  the  axis  of 
the  cylinder  and  decay  exponentially  as  |J6  I  increases. 

One  can,  of  course,  also  construct  waves  progressing  along 
the  cylinder. 

Ursell  calls  such  modes  of  motion  "trapping  modes" 
since,  if  they  occur  in  a  canal  with  sides  given  by 

5  ^  J  ‘  ■  n°  •,n'r»  l*  r*dl*”d  a“*y’ 

even  though  there  is  a  path  of  escape.  In  fact,  the  motion 
le  similar  In  this  respect  to  standing  waves  in  a  basin  of 
finite  extern:.  The  edge  waves  considered  above  also  tan  be 
used  to  construct  trapping  modes. 


YWMMMM  . . . 


properties  of  (18.36)  exists  only  for  ang.'es  Jr*45>#  end 
jf»  SO*  .  This  does  not  exclude  the  possibility  of  the 
existence  for  other  angles  of  a  periodic  progressive  wave 
with  a  curved  wave  front,  for  these  would  not  be  described 
by  the  assumed  form  O'f  « 

The  solutions  for  can  be  obtained  from  the 

fundamental  solutions  of  (13.6),  but  it  is  nearly  as  easy 
to  find  them  directly.  In  the  third  formula  of  (13.6)  let 
<%  *  Jr  m  ,  4c2  *  2  .  This  gives  the  velocity 

potential,  after  foming  a  progressive  wave, 

-j=r(y+4i)cett  c*>(U.ytt) .  (18.37) 


Let  the  sides  of  the  canal  be  given  by  JC-'-C  .  Then 

it  is  easy  to  verify  that 


L\v,-c - 4**t 

so  that  the  boundary  conditions  are  all  satisfied.  Since 
<2»  •»  /Vkv®  ^  ^ 


the  wave  velocity  la  given  by 

£  *  l^y^Csy  y=^“  1  (18.38) 

’'V 

If  *%v#  p  avw;  (in  the  notation  of  eq.(13.6)),  there 

will  be  t  further  symmetric  modes.  In  (13.6),  formula  4, 
set  1  A  and  add  this  to  formula  i  with  A  ,4rm& 


iMim4n 


-MM 


i  . . . . .  1  ” 
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This  give* 

f(&,*ji}lt)mA[C0a**i  th)ooUv  ifiifrWj? X  f  m0(4j  +4J a® {+£-($ "'jtj Cm(4t, .$ -<tf) 


One  may  again  verify  easily  that  o  on  the  two  aides 

of  the  canal  if  4t  m  4M>o  ~  ^  ■  .  Hence,  this  mode 

of  motion  will  exist  only  if  a^#v  y  *1*^  ,  For  given 

there  will  be  no  modes  of  this  sort  if  4,  is  small  enough, 
for  then  .  The  number.  gradually  increases  as 

'iv  increases.  If  4  and  it  ere  fixed  and  d  allowed 
to  Increase,  there  will  be  an  infin'tj  sequence  d,,dv 
for  which  4*.  **/'-**>  'will  be  satisfied;  C^l<$ 
as  4 ..  -*  oo  ,  The  situation  is  easily  visualized  by  plotting 
on  one  graph  41*  4  and  (d>u^jc^)  j  .  One 

may  write  the  potential  function  in  the  form 


f  rt)  m  A  [t*3  +-U,)c&s  M*,\  ij c+3(4tj- ii 


where 


AM,0  bcuUv  Mk„  c  V  *>vt  l&w  4wt-4/  *  ~V  f  it.  a4vt„ 


(18 . 39) 


The  velocity  is  given  by 


t-  ^A*u  /*&w(v 


-v-v. 


4*v; 


(18.40) 


Asymmetric  modes  of  motion  also  exist,  having  first  been 
noticed  by  GrcenhJll  [1886].  These  cannot  be  deduced  from 
(13.6)  but  must  be  found  directly.  The  velocity  potential 
corresponding  to  (18.37)  is 

uUy  Jfiy  c*(4tyXi)  (19.41) 
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The  wave  velocity  la 


cL-^2_.  c&tL 


(18.42) 


which  approaches  Infinity  as  >t*.  *L  -*■  o  .  In  addition  to 
this  node,  other  asymmetric  modes  may  exist  under  condition® 
similar  to  those  required  for  '18.39).  The  velocity  potential 
for  these  modes  is 

y,  §tt)M  4  L)  £aaM  •  (18.43) 

■  Cm  (i t  1-  <tf  } , 

wnere 

,  J  l  X.  \n"r 

K(y  'll/  *  ^  (  H'/  C4H  4*v  /  ^  i  *"  . 


The  velocity  of  propagation  la  given  by 

« 

t  /Wo  P^(»  <Kq/vo 


-W-o1 


‘o  Aa. 


(18.44) 


The  solution  for  the  angle  f *  will  not  be  discussed 

here.  It  can  be  found  in  Lamb's  Hydrodynamics  [1932,  p.  449] 
as  well  as  in  Macdonald's  paper  cited  above. 

One  may  construct  other  possible  contours  for  the  canal 
cross-section  by  starting  from  one  of  the  solutions  (13.5)  or 
(13.6)  and  finding  surfaces  for  which  ^  »<D  .  Thus, 

from  the  third  equation  of  (13.5)  for.. 

f  -  /4e.V?  oaaJU  jO  fL1-^  Csn  -£i)  . 


r 


IM 


Solution  of  the  differential  equation  c!y/Jx  *4^/4* 
leads  easily  to 


■ . -•■I——-  c*a JC  J  "V 

t  _  y *■  # 


as  an  equation  for  the  contour  of  a  possible  canal.  The 
contour  is  reasonably  shaped  but  varies  with  the  choice 
of  ■&*  ,  Also,  the  me  thou  is  unsatisfactory  In  that  it 

gives  no  information  about  other  possible  modes  of  motion. 
19.  Problems  with  steadily  oscillating  boundaries. 

Such  problems  include  waves  resulting  from  forced 
oscillation  of  a  submerged  body  and  the  waves  associated 
with  steady  oscillations  of  a  freely  floating  body  in 
oncoming  waves,  In  this  section  vc  shall  assume  the  fluid 
of  infinite  extent.  Waves  in  an  oscillating  bounded  basin 
will  be  discussed  later.  The  mathematical  treatment  has 
much  in  common  with  that  of  the  last  two  sections,  the 
scattered  wave  of  those  sections  becoming  the  forced  wave 
of  this  one. 

19  K  .  Forced  oscillations. 

Suppose  that  the  surface  of  the  oscillator  in  its 
equilibrium  position  is  represented  by  )  »  0 

Let  us  take  it,  for  example,  to  be  oscillating  vertically 

with  amplitude  £  .  Then  th*  equation  of  the  oscillating 
surface  $  may  be  written 

where  O-  is  some  length  dimension  oC  the  osciliatoi 
This  £  will  be  taken  as  th*  suspension  parameter  in  the 
perturbation  procedure.  U  the  perturbation  theory  of 


2 'i 


section  10,  we  were  concerned  only  with  the  function* 

4(<X  f,t)  and  *f,(x,  ,  t)  .  However,  we 

most  similarly  expand  f  before  substituting  it  into  the 
boundary  condition  satisfied  on  the  surface  of  the  oscillator, 
namely, 

(19.1) 


Ft  -O  F(x,yyt) 

The  expansion  for  this  case  is 


•O 


F(x,yt£^4tj))  •Ffay,  j)  t  FyyfX'f^yj  19.2) 


We  don't  wish  to  restrict  ourselves  to  this  one  mode  of 

motion  for  the  oscillator,  but  an  examination  of  the  form 
of  this  and  similar  expansions  indicates  that  we  may  assume 

in  general  that  the  surface  of  the  oscillator  can  be 
represented  by  the  series 

ip) 


\o }  |*  0/  ft ) 

ftJC''3')'*)  *  r  C*,y,  })  +  *[£  (19.3) 


+■  time-periodic  terns  in  higher  pawn  s  of  £  - 

Co) 

where  F  •*.  o  £#  t^e  e<qtall  brlum  position  of  the 


oscillator.  We  may  now  assume  either  that  ^  is  periodic, 
l.e. , 


or*  more  a  imply,  that  It  la  simple  harmonic, 


(19.4) 
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(19.5) 
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where  each  function  Tf**.  or  *f-  is  a  till  to  be  expanded 
In  ®  jjc jt t*umrb a t ioini  ®(srl.^?s  i*  TfcMS!  two  assumptions  are  not  quite 
equivalent,  even  for  the  first-order  theory,  but  si  nee 
under  certain  conditions  (18.4)  leads  to  the  same  first-order 
equations  ms  (18.5),  we  shall  *ssume  the  latter  form,  together 
with 

J  i  F')  ^  t(x,  1,(^1  j)  &*»*>  &F  Q9  gj 

Substitution  of  the  perturbation  series  into  the  exact 
equation'  and  boundary  conditions,  ms  in  section  10,  then 
leads  to  the  following  first-order  equation  and  boundary 

conditions: 

JO 

1)  (  K  -  i,l, 

2)  y^U.o,))-  ^  ft ’(*.<>, ))  -o  ,  k  • i , 

3)  (frraxb  F  •  <?!”+  i  F<?-0  Or,  ftosj'jyo,  (19.7) 

4)  F(0>*  <p<U  p('>  .  o  o*  F(jt,yj)-0' 

One  should  note  that  it  is  a  natural  consequence  of 
the*  method  that  the  boundary  condition  on  the  oscillator 
is  to  be  satisfied  at  its  equilibrium  position,,  lu  we  let 

'1.'Jr'1'5,'|yWP«|  '  **4»)*j (19.8) 


then  conditions  3)  and  4)  of  (18.7)  may  be  written 
^  =  A  (*,  J  )  oh,  F  ) «  O  f 

where  'f  -  'f,  -f  i  o><*cL  Am  A  >  +  ■ 


We  shall  henceforth  drop  the  superscripts  and  consider  only 
the  first-order  equations.  In  addition  to  equations  (19.7) 
the  functions  f/(‘  must  also  satisfy  the  usual  conditions  on 
fixed  solid  boundaries,  ,  and,  if  the  fluid  l« 

infinitely  deep t  j ^ f -*•  O  as  .  Finally,  one 

needs  a  boundary  condition  to  ensure  only  outgoing  waves  at 
Infinity.  As  has  been  pointed  out  by  Ursell  [1951],  the 
foregoing  conditions  are  not  always  sufficient  to  guarantee 
uniqueness  of  solution. 

Kochin  [1939,  1940]  has  considered  the  general 
mathematical  problem  In  water  of  infinite  depth  for  both 
two  and  three  dimensions.  Haskind  [1942b  ,  1944,  1946]  has 
extended  Kochin’ s  methods  to  water  of  finite  depth.  The 
frequently-cited  paper  by  John  [1950]  treats  the  theoretical 
aspects  of  the  problem  in  a  thorough  manner  and  Includes  many 
of  the  results  of  Kochin  and  Haskind.  Special  problems  have 
been  considered  by  numerous  authors.  Havelock  [1929k  ]  considers 
the  waves  generated  by  oscillation  of  a  vertical  plate  extend¬ 
ing  to  the  bottom  in  water  of  infinite  depth  for  both  two 
and  three  dimensions,  and  in  water  of  finite  depth  for  two 
dimensions;  he  also  considers  waves  generated  by  oscillations 
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of  a  vertical  cylinder.  MacCa^y  [1957]  has  treated  the 
three-dimensional  problem  in  water  of  finite  depth. 

Kenoard  [1949]  haa  treated  the  two-dimensional  problem  aa 
an  Initial -value  problem.  Uracil  [1948]  has  considered 
waves  generated  by  oscillation  of  a  vertical  strip  with 
finite  depth  of  immersion  in  water  of  infinite  depth; 
the  treatment  is  two  dimensional.  In  a  later  paper 
Ursell  [1949fe  ]  considered  the  waves  generated  by  the 
rolling  of  long  cylinders  of  ship-like  cross-section.  In 
addition,  Ursell  has  treated  the  waves  generated  by  a 
heaving  half-submerged  circular  cylinder  (1949a,  1953c, 

1954]  and  by  a  pulsing  submerged  cylinder  [1950]. 

Havelock  [1955]  has  treated  the  wave  motion  generated  by 
a  half-submerged  sphere.  Certain  mathematical  aspects  of 
this  last  problem  have  been  examined  in  more  detail  by 
MacCamy  [1954j.  Because  of  its  Interest  in  connection 
with  the  heaving  motion  of  a  ship  there  exists  many  papers 
attempting  to  compute  approximately  the  force  and  moment 
on  a  heaving  shiplike  body  resulting  from  wave  formation. 

We  mention  particularly  one  by  Grim  [1953].  In  the  cited 
papers  by  Kochin  and  Haskind  certain  special  problems  are 
solved  approximately.  In  addition,  Haskind  [1942,  1943b  ] 
has  considered  the  motion  resulting  from  forced  oscillation 
of  a  plate,  or  a  system  of  plates,  on  the  surface.  In  a 
more  recent  paper  Haskinr  [1953a]  has  developed  a  method  for 
finding  solutions,  and  in  particu  ar  the  force  and  moment 
on  the  body,  for  a  wide  class  of  two-dimensional  contours 
cf  ship-like  cross-section.  One  should  also  consult  a  recent 
expository  paper  by  Maruo  [1957]. 


This  brief  swill, .ary1  of  papers  on  forced  water  wave* 

Is  by  no  means  complete  but  lists  many  of  the  Important 
papers  and  l.uUcatea  the  richness  of  the  literature. 

As  was  tutted  in  the  introductory  remarks,  the  theory 
of  forced  water  waves  Is  mathematically  almost  identical 
with  the  diffraction  theory.  If  one  interpret®  the  value 
of  -  /3*v  on  the  body  as  the  function  describing  the 

motion  of  the  oscillator,  then  it  is  clear  that  the  problem® 
are-  the  same.  Hence,  the  general  remarks  about  existence  of 
solutions,  uniqueness,  and  special  methods  carry  over  directly 
and  will  not  be  repeated.  However,  we  wish  to  consider  here 
one  further  topic  in  the  general  theory. 

Kochln * a  H- function.  The  H- function  was  apparently 
first  introduced  by  Kochln  [1937]  in  connection  with  the 
theory  of  wave  resistance.  He  later  extended  it  [1939,  1940) 
to  waves  generated  by  oscillating  bodies,  and  it  has  become 
a  standard  device  among  other  Russian  workers  in  this  field, 
especially  Haskind,  who  has  extended  its  definition  to  other 
situations. 

Each  potential  function  satisfying  the  boundary 
conditions  has  associated  with  it  an  H- function  which  is 
related  to  it  much  in  the  same  way  that  the  Fourier  transform 
of  a  function  is  related  to  the  function.  One  of  its  chief 
virtues  is  that  it  allows  one  to  give  compact  formulas  for 
force  and  moment  on  an  oscillating  body  (in  the  present 
context)  as  well  as  certain  other  quantities.  It  is  also 
sometimes  helpful,  in  suggesting  approximate  solutions. 


II* 


Let  us  suppose  that  the  surface  ^  of  a  body  of 
bounded  extent  ia  oscillating  in  some  manner  in  fluid  of 
infinite  depth  and  let  <?.  f,  +L  be  the  solution  to 

the  potential-theory  problem  formulated  earlier.  Let  and 
$2  be  two  closed  surfaces  lying  below  with  •^closing 

S|  and  Sj  enclosing  S.  Let  us  denote  the  source  potential 
Introduced  in  (13.17")  by  ; £ ,  *[,  ))  , 

wliere  Q ,  ,  j)  are  the  coordinates  of  the  singularity, 

and  let  us  write  it  as  a  contour  integral: 


Af  /  1  I 


10) 


-I  0(C) 


where  the  path  L  passes  below  the  singularity  at 


4l  »  v-  6lfy  .  (The  residue  at  this  point  gives 

exactly  the  imaginary  part  of  (13.17").) 

Now  apply  Green's  Theorem  to  the  region  between  Sx  and 
S2  (the  following  argument  is  very  similar  to  a  two- 
dimensional  one  used  in  section  17 cv  in  discussing  the 
integral- equation  mef.nod) : 


(19.11) 


^U'V)  %&(?)]  it (f[r  |£  -  f$i(r)]  d6 

5,  Sz 

Then  iff  may  be  extended  to  &  function  harmonic  in  the 
whole  space  exterior  to  sl*  ft  is  harmon.'c  in  the  whole 
interior  of  S^,  but  since  S2  may  be  indefinitely  enlarged  as 
long  as  it  remains  below  /tj  -  O  ,  may  be  extended  to 
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bit  harmonic  in  tha  who  la  half- space ,  i  0  *  Consider  now 


the  functi  ,  on 
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satisfies  the  free -surface  condition  and  the  condition 


at;  infinity.  Moreover 


C Jf  :x  f- 


harmonic  in  the  lower  ha  1  f- space »  (f-  Y  is  harmonic  in  the 
lower  half-space  and  satisfies  the  other  boundary  conditions. 


But  then  Cf-Y3*0 


as  follows  from  a  uniqueness  theorem. 


proved  by  Kochin  [  1940 ,  5  1],  Hence,  we  may  write 
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Now  define 

Ifjri  7l 

H(4  (9)“J/«  (U(},'},)) 

s', 


4  <f{0*^(  'K  , '*7 )  f  l  C4K>  6^  6#<#  5  )  +■  ^  jstd 


gJtll 

k 


4,  (  1  +  L  ^  60»  o  +  l'  7  Vav 


'Hx 
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(19.14) 
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tw'O*  r  some  manipulation  with  (19.13),  one  can  show 
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W®  give  a  few  of  Kochin’a  derived  f  omul  as.  The 
••Wtotie  form  of  the  free  surface  in  a  direction  o( 
l«  given  by 

1  *  (19.16) 


The  rate  at  which  energy  is  being  carried  off  by  the 
wave#  (and  hence  also  the  power  Input)  la  given  by 

N-ra  -y-1/ /h(»,*;/W. 


(19.17) 


The  force  components  on  the  oscillating  body,  averaged 
a  period,  ere  given  by 

X„  -  jf  7 lH(v, 9)  I  1c*9jl9 

*  JT 

t  -nv <{-,'£ 

1  r  V 

Z„  •£*'/ IH<V,9)I  Ua~&  d& . 
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over 


(19.18) 


The  formulas  can  be  derived  from  (8.4),  (9.4),  and  asymptotic 
express  ions  for  /  . 


In  formulas  (19.14)  and  (19.15)  the  surface  Sj  over 
which  the  integrals  are  taken  may  be  contracted  to  S,  This 
some times  makes  It  possible  to  express  H  directly  in  terms 
of  known  boundary  values.  If  f  can  be  expressed  by  means 
of  a  source  distribution,  say 
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(19.19) 


then  on*  'ha* 
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In  order  to  find  approximate  answers,  Kochin  frequently 
uses  the  distribution  ,;jr  which  would  be  proper  In  an 
infinite  fluid  without  free  surface,  substitutes  this  in 
(19,20')  and  then  uses  the  resulting  approximation  to  M 
in  (19.17)  and  (19.18)  above.  The  procedure  may  be  looked 
upon  as  the  first  two  steps  In  an  alternating  type  of 
approximation  in  which  one  first  satisfies  the  boundary 
condition  on  the  body,  neglecting  the  free  surface,  next 
corrects  this  so  as  to  satisfy  the  free -surface  condition  , 
but  now  disturbing  the  condition  on  the  body,  then  corrects 
again  to  satisfy  the  condition  on  the  body,  etc.  This 
method  of  approximation  has  frequently  been  used  by  Have¬ 
lock  fe.g.,  1929a  1 . 

Kochin  [1939]  ha®  also  defined  the  H- function  for 
two-dimensional  wave  motion  excited  by  an  oscillating  body. 
We  simply  reproduce  the  formulas.  Let,  as  usual, 

+  b* the 
complex  potential  and  let  Cj  and  Cg  be  two  contours  in  the 

lower  half-plane  containing  C,  CL  Inside  C,.  Define 
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where  W^s  Ht  .  This  follows  Immediately  from  <s,  fc-mula 
similar  to  (17.15).  For  Che  asymptotic  form  of  the  waves 


one  get a 
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as  x-++oo  j  (19.23) 
t,  (*u  +  «££) 

as  X  -r  -  omO1  , 


The  rate  of  dissipation  of  energy  Is 

H-{  f*i[lH,l?>)l''+IHL(v)  r]  .  (19-24) 


The  mean  value  of  the  force  and  moment,  averaged  over  a 
period,  la 


^C'V  f  X  fyv\/  (  H,(v)  Hx(y)j , 

r«  %T  pf-£${lV*>l',+lHiVl''}  die  , 

i®»0 
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+  2  v  ^  |  W/6»)  J-Jttv)  ~  M  t(v)  |  . 

Roughly  the  same  remarks  apply  to  the  use  of  the  two¬ 


(19.25) 


dimensional  formulas  as  of  the  three-dimensional  ones. 


Hflg&flL-faSBLJn^gellletor  in  a  wall.  In  order  to 
Illustrate  the  use  of  the  H- function,  we  consider  the 
following  problem.  Let  the  («  ;  j) -plane  be  a  rigid  wall 


I25 


axi-ept  for  41  certain  bounded  area  S  in  which  there  Is  1 


membrane  oscillating  according  to  a  given  law 

X9  ?(*)•!)  ***■  vf  ,  fyj)  f"  S- 


(19.26) 


The  boundary  condition  which  has  to  be  satisfied  on  the 
plane  Ji  -  q  Is  then 


(19.27) 


where  we  still  have 


<A- 


This  boundary  condition,  as  well  as  the  ones  at 
Infinity,  can  be  satisfied  by  distributing  "modified"  sources 
(13.17")  or  (19.10)  over  S  with  density  -  4  F j)  /'2  JTf 


‘w  ffF(Tt )£(*']' }i  1 1 })  . 


(19.26) 


In  order  to  compute  the  H- function,  we  shall  Interpret 
the  source  distribution  as  representing  a  thin  body 
making  symmetric  pulsations  in  an  infinite  fluid.  Hence, 
we  may  as  time  that:  the  wall  is  removed  and  the  membrane 
replaced  by  a  doubled  one.  (That  the  requisite  motion  is 
physically  impossible  doesn't  Invalidate  the  considerations; 
a  more  realistic  model  can  easily  be  devised.)  In  (19 . 1*2) 
we  take  to  be  both  sides  of  the  thin  body.  Then, 
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remembering  that 
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one  finds  easily  that 

2.4,  j)C  d«j  dj  .  (19.29) 


From  (19.17)  one  then  finds  immediately  after  carrying 
out  the  &  integration,  that  the  rate  of  dissipation  of 
energy  to  one  side  is  given  by 

A/«  *  ffdjdjffdydfFfopFtyfit^^L&yi))  (19.30) 

s  s 

Expressions  for  and  Z.^  may  al»°  b®  written  down. 

The  result  Xcur  *  0  ls  not  reaUy  significant  because  the 
integral  is  over  both  sides  of  the  thin  pulsing  body. 

The  theory  for  genera'  on  of  two-dimensional  waves  in 
a  semi-infinite  channel  by  a  vertical  wave  maker  in  the 
end-wall  is  easily  derived  in  the  same  way.  If  the  motion 
of  the  wave-maker  is  described  by 


x  *  F(  <^)  At  f  O-  **  ^  ^ 

then  4 

u,(4t)  •jf'1  F(n)d*i  , 
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and,  for  cs xani pic ,  tlhic  rate  of  dijJS'*.4  acton  of  ©nair^y  ta 
given  by 

^  V  V 

W  *  ?  i  le  '  F(Vdj  j  (19.33) 

The  generation  of  short-crested  wav“s  is  subject  to 
the  limitations  described  In  section  14  f  .  Suppose,  for 
example,  that  the  water  is  of  depth  h,  ,  the  channel  of 
breadth  ,  and  that  the  motion  of  the  wave-maker  is  described 
by 

X -*  F(*j)  oso 4^  stu, to  t  (  4  *  4t.IT/4  ,  '  C  '  o  ,  (19.34) 

Then,  since  c<«4v., .  (<^  +  L)  ,  (^t4.)  form  a  complete 

set  of  functions  in  -  -4  6  <  o  ,  there  is  no  difficulty  in 

representing  F )  by  a  series  of  the  fundamental  solutions 
o 

(13.6),  but  if  &*>*%v*v  ,  no  progressive  waves  will  move 

down  the  tank  (within  the  limits  of  applicability  of  the 
linearized  theory,  of  course).  The  analysis  of  the  filtering 
effect  of  the  tank  on  more  complicated  wave-maker  motions 
can  easily  be  carried  through  by  Fourier  analysis. 

Waves  from  an  oscillator  not  in  a  wall.  Let  us  now 
suppose  that  we  have  a  two-dimensional  oscillator  m 
infinitely  deep  water  moving  according  to  the  law 

Ftys^il  ,  (19.35) 

but  with  no  wall  present.  This  small  change  complicates 
the  solution  of  the  problem  in  a  substantial  way,  the 
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complication  being  associated  with  the  now  possible  flow 
under  (and  over  if  l<  o  )  the  oscillator.  In  addition,  in 
order  to  ensure  a  unique  solution  some  further  condition 
analogous  to  the  Kutta-Joukowski  condition  in  airfoil  theory 
is  required;  such  a  condition  does  not  seem  to  have  been 
formulated.  Aside  from  this  lack,  the  boundary  condltl  -ns 
to  be  satisfied  on  the  oscillator  by  the  velocity  potential 
*  f,  cm  &  i  4  f2  4>t  are 


(19.36) 


7*(°i  t)  m 
jpO,  OL,b)  *  o, 

4>y(0<  t,t)  •  o  u/  <>  <0  , 

The  problem  is  clearly  closely  related  to  that  of 
diffraction  of  plane  waves  by  a  vertical  barrier  and  could 
be  treated  by  a  modification  of  the  method  used  in  section  17*f 
for  that  problem .  It  may  also  be  solved  by  the  integral- 
equation  method  discussed  in  section  17  c*  .  A  modification 
of  this  method  has  oeen  used  by  Ursell  [1948]. 


where 


Introduce  the  complex  potential 

$  +  1  f  -  Re  j  ft  (pc'***  J 


(19.37) 
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We  try  to  construct  a  solution  by  means  of  a  distribution  of 
vortices  of  the  form  (13.28) 
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wi  th  in  tens)  '  ty 
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along  the  oscillator: 


Q.  <  M  £  4 


jlP  '  j  A  l*])  /*  ( ]  >  u  7 )d 7 


(19.39) 


(19.40) 


An  analysis  almost  identical  with  that  in  17 x  leads 
quickly  Co  the  integral  equation 


l  i^}  /  /v  l  ~  '  j  i c  7  )^7 


W°  , 


(19.41) 


Separating  ^r(  and  and  noting  that  /y  (i  4  ,  i*j) 

is  real  with  respect  to  i  ,  one  finds 


J ( y/7)fvi  fvt  (i  7  * 1 V J  ^  7  •  ^  F(y)  t 

6 

if  iW(w  A>lV  + W/«  (:v  '\)]drl  - 5  • 


(19.42) 


The  equations  can  be  uncoupled  by  applying  the  operator 
I  9/3^  -■  v  ]  to  each  (so  that  the  reduction  method  enters 
In  after  all!).  Introducing 

y^'£>*  ^4,  “  v  ^  ,  G(*j)  •  h  i?  F t  0 

one  finally  obtains  the  pair  of  equations 

lAt 7>vi  -  M  -  JT<  ^ 


(19.43) 


(19.44) 


lAere  we  have  taken  advantage  of  the  fact  that 

0f  *  '»  (f  f-y)  end  hence  fta)  -  C>  ;  if 

4<*0  alio  f(4t)*0  ,  The  integral  equations  are 

easily  reduced  to  a  known  t>ne  occurring  ir.  airfoil,  theory 
[see,  e.g.,  W.  Sch raeidler,  uitegvalgleichungen 
Akademlache  Verlagsgesellachaft,  Leipzig,  1950,  pp.  55-56, 
or  Mikhlin,  Integral " nye  uravneniya  . . . }  Gostekhlzdat , 
Moscow,  1949,  pp,  149-154]  by  the  transformation 

r  -  +  f « ^ {(ax^  l) ' 

The  solution  may  be  written  in  the  form 


A-W *  i i* u 4  r1s,'i> jrtj 


/WV<? 


(19.45) 


It  Is  evident  that  the  solution  i  i  not  uniquely  determined 
without  some  statement  about  the  total  circulation.  Fixing 
the  total  circulation  Is  equivalent  to  fixing  l[  (JL)  ,  as 
follows  easily  from  the  form  of  JU  (\)  and  the  relation. 


j  t Vli  u(t) ds 


(19.46) 


It  is  possible  to  compute  the  H-function  directly  in 


terms  of 


First,  we  note  that 


<f&  ^  df  ($(*]){ v  (l )  ll)^1 

c,  c.  **  9 

i  t 

*  (u (1(1 

A  *  'L  *  A  f 
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It  then  follow#  from  (19.46)  that 
■A-(  * 

.HO)  -  r(i)  -  jrzfAWe  ‘ .  (19.47) 

One  may  now  apply  formula#  (19.23)  to  (19.25)  to  find  the 
quantities  Indicated  there  (note  that  W  *  H  ). 

One  notes  again  that  the  function  M  fX)  la  determined 
only  alter  |(4)  Is  fixed.  Taking  \(k)  f  O  Is  equivalent 
to  having  a  singularity  at  the  end.  If  the  oscillator  Is 
totally  submerged.  It  seems  reasonable  to  set  jf (4)**0  , 

a®  we  assumed  in  (19.36),  for  then  the  vertical  velocity  is 
continuous  at  the  end,  i.e.,  ^ (+o,  A)m(f^(- o,  4>)  ,  as  has 

already  been  assumed  for  the  lower  end  at  -y «  A .  It  is  not 
clear  what  is  the  proper  assumption  if  4*0  ,  i.e.,  if  the 

oscillator  extend®  through  the  surface.  In  the  similar 
problem  of  diffraction  by  a  vertical  plate,  treated  by 
the  reducl Lon  method  in  section  17 or,  the  assumption  of  no 
singularity  at  the  surface  is  ecm-'-alent  to  assuming 
)((0)  ss.Q  .  We  note  that  li  '(&)  *  O  ,  then  it 
follow®  from  (19.46)  and  the  form  of  in  (19.45)  that 

a o  ,  and  hence  that  fas  O  •  This  is  not  true,  of 
course ,  for  | |  . 

Waves  generated  by  a  heaving  hemisphere.  We  describe 
briefly  a  procedure  used  by  Havelock  [1955]  and  MacCany[1954] 
and  before  them  alao  by  Ur sell  [1949a  1  for  an  analogous 
twO'^ilsienslopal  problem.  let  a  hemisphere  of  radius  Ou 
have  Its  cancer  on  the  free  surface  in  its  undisturbed 
position  and  let  It  undergo  forced  vertical  oscillations 


described  by 


xV 


(y~  i,  Um^  6i)l+  - 


(19.48) 


Then  the  boundary  condition  to  be  satisfied  by 

-  f,  +  i  ^  on  Lu«t  hemisphere  is 

on  yto.  (19.49) 


*f  must,  of  course,  also  satisfy  the  free-surface  condition 
and  the  radiation  condition,  as  stated  in  (19.7). 

The  method  of  the  above-named  authors  is  to  represent 
Cf  as  a  series  In  which  the  first  term  is  a  source  at 
the  center,  say  (13.17),  and  the  remaining  terms  represent, 
only  local  disturbances  of  the  sort  shewn  in  (13.21),  with 
m  -  O  since  we  have  radial  symmetry.  The  source  term  is 
actually  taken  in  the  form  (13.17" ').  Since  the  source  is 
r*K  (O,  0,  o  )  t  r  *  r,  In  the  formulas  and  certain  terms 
cancel  and  others  double.  Let 

«  1-  uf[ 

Y0(vR)tiTv  I*!  Solvit.)  ,  (19.50) 

-v  P..^  (c+°  . 
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Then  the  assumed  form  for  ti?  is 
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Substitution  In  th«  boundary  condition  (19.49)  leads  to 
««  infinite  net  of  linear  equation#  for  the  coefficient® 

»  B  tu  *  Numerical  method®  may  then  he  ueed  to  find 
any  desired  number  of  term®.  MacCamy  thovi  that  the  solution 
(19.51)  is  convergent  if  va,  <  ** !j 

laving  found  ^  approximately,  one  may  proceed  to 
compute  the  vertical  hydrodynamic  force  on  the  sphere  by 
integrating  the  pressure  £> «  -  j$t  over  the  hemisphere. 
Havelock  carried  through  an  approximate  calculation,  expressing 
the  result  in  the  font 


Y-  £  w  i v  [A  Cck.  4  l  ~  2  4/  flee 

,  -  nit.Xk-n-  iLt  &  , 


(19.52) 


■ap-  •  ■  —  m 

where  M  is  the  mass  of  displaced  fluid  and  the 

coordinate  of  the  center.  The  parameter  is  usually  called 
the  added-mass  coefficient  or  the  virtual-inertia  coefficient; 
k>  i#  called  the  damping  parameter.  Figure  18  from  Havelock's 
paper  shows  A  and  2  L  as  functions  of  v  a.  . 


average  rate  at  which  work  la  being  done  by  the  sphere  Is 
k>  and  does  not  involve  ft  . 

It  is  of  interest  to  compare  the  same  parameters  as 
computed  by  Ursell  [ 1949a  ]  for  a  circular  cylinder 
(per  unit  length).  The  nre  shown  in  Figure  19. 


Figure  1  ^ 

The  asymptotic  behavior  of  4  is  given  by 

(IS. 53) 

as  w-40j 

fc.(va)  «  I  -  +  0 (yT£  )  * 5  va  ^  ■ 
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If  f£>  .  S toady  aaclllttloos  of  a  freely  floating  body  In  wave*. 

Lot  us  suppose  that  «  rigid  body  1«  floating  in  an  Infinite 
ocean  with  proscribed  plane  wave®  approaching  from  a  fixed 
■direction,  say  from  X  *  +  +*>  .  If  the  motion  has  persisted 

for  some  time,  we,  mey  suppose  that  the  body  is  moving  with  a 
staple  periodic  motion  of  the  same  frequency  os  the  waves. 

With  this  assumption  the  proper  formulation  of  the  linearized 
equation*  and  boundary  conditions  has  been  derived  by  John  [19491 

Suppose  the  body  is  at  re»t  In  still  water  and  let 
(  T  '  d  *  ^  th*  coor<*inate»  of  its  center  of  gravity. 

Let  6  x  Cj  J  be  a  coordinate  system  fixed  in  the  boCv  with 
Q  at  the  center  of  gravity  and  the  axes  parallel  to  the 
space  axes  0  X  y  j  ,  When  the  body  is  displaced,  one  may 
describe  its  position,  by  giving  the  new  position  of  the 
center  of  gravity  t  *j }  jj! )  -  (x,  +  CX,,  f  l  +  £  Ji) 

and  the  Eulerian  angles  £*,  if),  (we  change  notation 

from  the  customary  f,&,  Y  to  avoid  confusion  with  our 
other  use  of  these  letters) .  Thus  the  choice  of  £  implies 
that  the  amplitude  of  motion  is  small  compared  to  some  typical, 
body  length.  The  assumption  of  section  ^0  <x  that  £  4  +  ’" 
implies  that  the  amplitude  of  the  prescribed  incoming  waves 
i»  also  small  compared  to  this  length.  The  relationship 
between  the  two  sets  of  coordinates  may  be  easily  found  from 
the  usual  formulae  eOMh  tning  Eulerian  angles  to  be  of  the 
font 
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<£  (xi  ~  *  y*)  +  /b(^  -  y)J  f  £'[■  ■  ■  j  +  *  *  * /  (19 . 54) 
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Let  the  surface  of  the  body  be  described  by 

F(jt,y,j)*°  (19.55) 

in  body  coordinates.  To  find  the  position  In  space 
coordinates  one  must  substitute  from  (19.54)  in  (19.55). 

The  kinematic  boundary  condition  (see  (19.1))  then  becomes 


Aj  -  /j  -  _  £jf|(jC~X0) 


(19 . 56) 

i[r<UF(x-xe,rri  5  yy)k**((rr>Myifl' 

t  €V  f  ■  j  +  •  •  ■  =0. 


Letting  rjjt,  be  the  components  of  the  unit 

normal  vector  to  the  surface  at  rest,  i.e., 


F  y.  ,  *0  (19.57) 

(we  shall  call  this  surface  $©  ),  and  fy-t?  ~ T9)  *  * t  ‘'e. 

•» 

J*) **)n] /  “(* "X) *y“/ y - y») f  (J9.5g) 


we  may  express  the  first-order  term  in  (19.56),  after  dropping 
the  superscript,  in  the  form 

=  h*  +  ^  i  ‘•Vy  i  j,  tm  *  ¥  f]  /cc  £*(  y,  o-h,  • 


. . mm . . . . 


We  call  attention  to  the  fact  that  it  follows  as  a  natural, 
consequent®  of  the  1  Inter lxation  that  the  boundary  condition 
Is  to  b«  satisfied  on  the  surface  in  Its  undisturbed  position, 
In  order  to  state  the  dynamical  conditions  on  the  body 
we  Introduce  the  following  notations.  Let  M  be  its  mass 

*"d  t" l* '  h '  l“-  hn  lca  and  mt”'anta  and 

products  of  inertia  about  the  body  axes  selected  above. 

Let  V  be  the  volume  bounded  by  the  plane  y  *  o  and  the 
submerged  part  of  the  surface  in  its  rest  position,  and 

»V  -  V  r  V  *V  .  V  *  V 

let  *  (  1^  i  ,  I.  ,  lMX  If.ut  ••  •  be  the  volume,  the 


i\r 


.  V 

1  j  •  lxx  ,  1 /  •  •  • 
moments,  and  the  moments  and  products  of  Inertia  of  this 

volume  about  the  body  axes  in  their  rest  position.  Let  A 


be  the  Intersection  of  the  body  in  its  rest  position  with 
surface  y  >o  »  and  let  1*  /*,  ]*f 
denote  the  area,  the  moment!  end  the  moments  and  products 


the 


of  Inertia  of  A  with  respect  to  a^es  through  (a9l  o, 

and  parallel  to  the  body  axes;  e.g., 

4 

/ <  j  *  ff(*~  “**)(]“  j*)  dx  dj  ■ 


) 


The  exact  dynamical  equations  are 

M  f‘jjp  o*s(nrx)  < Idf 


(19.60) 
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f-  -  f<3r  fil  ^  H 


and  three  similar  equations  for  oc  ,  jb  ,  jf  .  Substitution 
of  the  perturbation  series  gives  for  the  zero-order  tarns 


M-  fl 


•  V 


i.r  - 1 


(19.61) 


i.e.,  Archimedes 1  law  and  the  statement  that  the  center  of 
buoyancy  and  center  of  gravity  are  on  the  same  vertical  line, 
Ihe  first-order  equations,  after  dropping  superscripts,  are 


(19.62) 
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We  note  that  the  boundary  conditions  have  been 


derived  for  general  motions  of  the  body  and  fluid,  not 
just  for  the  simply  periodic  ones  for  which  they  will  be 
used  below. 

John  [1949]  has  used  the  equations  to  investigate  the 
stability  of  a  floating  body.  Wa  shall  not  reproduce 
the  results  but  remark  that  he  shows  that  the  usual  condi¬ 
tion  for  stability,  namely  that  the  me  t«  center  lie  above 


Ml, 


-ml  inw»fWMHWH  "W*“  I"  ' lll—"  1 '  . . 

,  rr  WW wrB  '^.awgwwiiwniwr; WMW«ll*«i|'»in WtWH^WWWIMIIIIKlfWWIWtriWWW^W M»MWIW*WI»IW»IWIM«WIW^>>S«^  f  IM WtllW WMIWH»|I 

m,  m  mm 

&  *#  f 

the  center  of  gravity,  derived  from  purely  hydrostatic 
considerations,  is  in  fact  still  a  sufficient  condition  for 
stability  when  the  hydrodynamic  equations  are  considered 
(within  the  limitations  of  th#  linearised  theory). 

It  is  also  shown  by  John  that  the  above  equations  have 
a  unique  solution  for  an  initial-value  problem,  i.e. ,  if 
at  some  instant  the  position  and  velocity  of  body  and 
fluid  are  prescribed.  Hfwever,  for  the  problem  with  which 
we  are  concerned  in  this  section,  steady  simple  harmonic 
motion  with  a  prescribed  incoming  wave,  he  proves  uniqueness 
only  for  sufficiently  large  values  of  <4  and  for  bodies 
such  that  a  vertical  line  intersects  the  immersed  surface 
only  once  (e.g.,  a  floating  sphere  with  its  center  at  or 
above  the  f  ree  si  face) . 

Knowledge  of  the  motion  of  e  floating  body  in  surface 
waves  Is  obviously  of  great  importance  to  ship  designers, 
and,  as  might  be  expected,  there  is  a  large  amount  of 
specialized  literature.  However,  most  of  this  literature 
may  be  considered  irrelevant  to  this  article  for  it  is  based 
upon  the  assumption  that  one  may  neglect  the  kinematic 
boundary  condition  (19.59)  completely  and,  thft  ciy  0.21x11 i.  c 
boundary  condition  (19.62),  that  one  may  take  for  <£> 
simply  the  velocity  potential  for  the  oncoming  wave,  thus 
neglecting  the  effect  of  the  diffracted  waves  and  the  waves 
generated  by  the  shlp'e  own  motion.  This  assumption  is 
usually  called  the  Froude -Krylov  Hypotheala .  W.  Fraud*  [1861] 
introduced  it  in  connection  with  an  investigation  of  ahip 

a  "■  .  ,  !  .  " 
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rolling  In  wtvnt  and  A.  N.  Krylov  [1896,  1898] 
investigated  Its  Implications  rather  thoroughly  for  general 
motion®.  In  spite  of  its  apparent  crudmess  this  assumption 
has  been  useful  in  elucidating  many  aspects  of  ship  motions. 

In  recent  years  there  have  appeared  a  number  of  papers 
in  which  an  attempt  has  been  made  to  take  account  of  the 
proper  boundary  conditions,  but  no  attempt  will  be  made  to 
summarize  this  literature.  The  most  systematic  investigation 
of  the  matter  has  been  made  by  John  [1949,  1950],  Haskind 
[1946a],  and  Peters  and  Stoker  [1957].  The  papers  by  John 
consider  the  proper  formulation  of  the  linearized  problem 
for  a  bod'*  with  no  average  forward  speed  and  the  uniqueness 
and  existence  of  solutions.  Both  Haskind  and  Peters  and 
Stoker  are  primarily  concerned  with  ships  having  a  constant 
average  forward  speed.  Peters  and  Stoker  treat  carefully 
the  proper  formulation  of  the  linearized  problem  and  conclude 
that  Haskind' s  fundamental  equations  are  not  properly  formu¬ 
lated  In  that  some  of  his  varms  really  belong  with  the 
second-order  terms  and  should  have  been  discarded.  The 
objection  applies  also  to  part  of  his  results  for  a 
stationary  ship.  The  other  part  will  be  summarized  below. 

The  motion  of  a  ship  in  waves  when  it  has  a  nonzero 
translational  velocity  will  not  be  considered  in  this  article. 
For  this  theory  one  should  refer  to  the  cited  papers,  to 
Stoker's  Water  waves  [1957,  ch.9]f  or  to  a  recent  survey  by 
Maruo  [1957].  The  transient  oscillatory  motion  of  a  float¬ 
ing  body  in  calm  water  will  be  considered  later. 
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tMt  u®  return  to  the  problem  of  steady  oscillation  of 
a  floating  body  in  oncoming  waves .  Since  we  assume  steady 
oscillation,  we  shall  write 


,  ^  jjj.  - 

f  *  Rt  I  f  t  f-t '1^-/3') RtfctU/l*,  C.)e  ^  J-)  j?e  (fa  fa  ft) 

where  /-  /,  f  i  ft  »  <3.0»  <x»  f  i  ol  ,  etc.  The 

unknown  function  f  .nd  ch.  con.t.nt.  a,  .  ■  ■ ,  g.  «re 

to  be  determined  from  the  equations  and  boundary  conditions 
We  shall  assume  that  ^  can  be  expressed  as  the 
sum  of  the  velocity  potentials  of  the  incoming  wave,  say 


(19.63) 
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w 

If  the  fluid  is  infinitely  deep,  a  diffracted  wave 
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ijjj'  •  f*^Ut  and  e  forced  wave  <^y  *  ff  C 

resulting  from  the  body's  own  motion: 


is  t 


(19.65) 


We  shall  express  <P|  in  the  following  form  (following 
Hasklnd) : 


f>t  -fi,+  fi  +  fyf*  i  + 

Then  the  kinematic  boundary  condition  (19.59)  implies: 
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all  Co  b®  satisfied  on  £0  ,  the  rest  position  of  the  body. 
The  function®  .  .  .  t  <*  ,  are  to  satisfy  also  the 

radiation  condition  and  the  condition  at  •  -  «»®  (or  at 
^  M  -  U,  for  a  flat  bottom).  The  dynamical  condition 


(19.62)  will  be  used  to  determine  the  amplitudes  and  phases 
(!.«.»  the  complex  amplitudes),  k.:  first  we  Introduce  some 


notation.  Let 


(19.68) 


constants  jU  and  X  depend  only  upon  the 


geometry  of  the  body.  It  may  be  shown  by  an  application 


of  Green’s  Theorem  that 


7V-. 


and  X/  "  X ,  t 


Let  us  now  substitute  the  expanded  expression  for  ^ 
Into,  say,  the  first  of  equation* (19 . 62)  (the  others  may  be 
treated  similarly) ,  remembering  that  on  S6  : 


M*;.  -f  H(r*n**“  -  r‘t)  <  d*  •  <19.6  9) 

S*  S. 

Consider,  for  example,  the  second  term  of  the  second  integral; 


§ ff  d6  m(/*n  f  ^  X2l)j,  + Xzl  /  (19.70) 


where  we  have  used  the  special  form  of  aj  , - e 
Thus,  (19.69)  may  be  written 
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,r®  rf*m  tC*^  and  h»  •]<»*£  represent 

the  x-conponenttof  the  force*  re  sal  ting  from  the  Inching 
end  diffracted  wevee  and  are  to  be  computed  front  the  first 
integral  in  (19.69).  The  form  of  (19.71)  explains  the 
names  given  to  the  yAi^  and  A-  :  the  are  called 


added  masses,  the  A; 


damning,  forces .  If  one  now  writes 


Jt 


•  ! 


in  their  assumed  forme  in  (19.63)  and  sub¬ 


stitutes  in  (19.71),  one  obtains 


>4\ 


{ n  )  A«  ^  ^ *  “idA^dte' 


(19.72) 

€\fo{c*  +  /.at 


of  symmetry  rather  than  on  the 


2  *11 


and  five  similar  equations.  Since  the  amplitudes  CLmf  •  •  -,#• 
are  complex,  this  gives  twelve  equations  to  determine  the 
twelve  unknown  quantities.  It  is  thus  clear  that,  providing 
these  equations  can  be  uniquely  solved,  the  problem  of  finding 
the  steady  oscillatory  motion  of  a  freely  floating  body  can 
be  reduced  to  the  solution  of  several  problems  of  che  type 
studied  In  sections  18  and  19  o<  .  From  the  form  of  (19.72) 
and  the  similar  equations,  it  is  clear  that  the  complex 
amplitudes  are  all  proportional  to  the  amplitude  A  of  the 
Incoming  wave  as  would  be  expected. 

Hasklnd  has  applied  the  method  outlined  above  to  a 
body  symmetric  with  respect  to  the  (x,  y) -plane,  e.g.»  a 
ship  heading  Into  waves.  The  only  possible  motion®  are 
heaving,  pitching  and  surging.  In  carrying  out  some 
numerical  «ooi;putetioi|i  he  make®  a  further  approximation 
that  the  kinematic  hdundwry  condition  on  the  ’  pdy  -may  be 


. . . 
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surface.  Although  this  approximation  is  perfectly  consistent 
with  the  linearised  theory  in  certain  contexts,  as  will  be 
seen  in  section  21,  it  is  not  consistent  with  the  theory  as 
formulated  here  and  must  be  considered  to  be  a  further 
approximation  of  some  sort. 

Freely  floating  sphere.  Computation  of  the  motion  of 
a  freely  floating  sphere  with  its  center  at  th**  undisturbed 
water  level  can  be  carried  through  without  an  unreasonable 
amount  of  numerical  work.  The  procedure  for  the  heaving 
motion  has  been  carried  up  to  the  point  of  numerical  compu¬ 
tation  by  MacCamy  [1954].  Part  of  the  problem  has  already 
been  solved  tr.  section  19 (X  }  i.e.,  the  waves  resulting  from 
the  forced  motion . 

Since  the  phase  at  infinity  must  be  kept  arbitrary,  one 
must  replace  (19.48)  by 

JC  Goo  4>t Samj  &  t)  f  ^  Q*  (19.72) 

However,  the  solution  of  that  problem  may  be  taken  over 
with  practically  no  change,  for  the  velocity  potential 
in  the  notation  of  (19.66)  must  satisfy 

I or  x  4  o  ■  (19.73) 

Thus  we  need  only  set  I  in  (19.49)  and  later. 

In  fact,  from  formula  (19.52) 


(19.74) 
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Finding  eh*  diffracted  wave  requires  finding  an  outgoing 
e  satisfying 
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A  r  7  v  A  Cf>  C‘  -  (19.75) 
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» re  JL  ■  r  suw  $  ceo  c<  ,  ^  » r  &  A  «■  r  clw  $  i*u.  .  HacCamy 
U sws  that  J  °  can  be  found  as  a  series  in  functions  of 
form  (13.21)  with  Z  «  O  and  account  taken  of  certain 


series,  Let 
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(f  may  be  expressed  as  follows: 


y  |c*a|  2'H»  —iytx  I  if  »  If  2f  ■  •  •  f 
¥n  ■  I »  •  ■  ife  , 


#  ?P  *2.*  24+1  „***  ^  2-Ht-i  ttf5 


(19.77) 
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I  with  the  coefficients  .  x 


to  be  determined 
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by  (19.75).  MacCamy  proves  that  the  series  converge  if 
Va  <■  II  /  "I  and  analyzes  the  coefficients  further. 
Unfortunately,  no  numerical  data  seems  to  be  available. 

20.  Motions  which  may  be  treated  as  steady  flows. 

In  this  section  we  shall  consider  several  problems 
which  are  time -independent,  either  by  their  formulation  or 
by  introduction  of  moving  coordinates.  The  flow  associated 
with  a  constant  discharge  rat#  through  a  canal  Is  of  the 


first  type;  the  waves  associated  with  a  ship  which  has  moved 
with  constant  velocity  C  over  a  long  period  is  typical 
of  the  second. 


The  boundary  conditions  at  the  free  surface  have  been 
derived  in  sections  10  and  11.  For  three-dimensional  motion 
the  velocity  potential  must  satisfy  (see  eq.  11.3) 


(20.1) 


the  equation  of  the  free  surface  is 


(20.2) 


In  two -dimensional  motion,  if  the  complex  potential 
j  •  Lfr  i  Y  is  used,  then  the  boundary  condition  may 


be  written 


Z,ff(x*io)n^{'U^o)  j.o 


(20.3) 


If  the  potential  has  been  taken  in  the  form  F(^) *  -  + 

with  *  0  as  the  free-surface  streamline,  then 

U*i0)j  *<?;  <20-3> 

the  surface  is  given  by 


(20 • 4) 


On  obstructions,  which  are  now  all  fixed,  one  has  as  usual 

iff  P  *  \  (20*5) 


Far  ahead  of,  or  far  up»tr»am  of,  the  obstruction 


tii*  notion  must  approach  either  rest  ,  or  a  uniform  flow, 
respectively. 

The  general  theory  of  steady  free- surface  flow  about  a 
submerged  obstacle  in  infinitely  deep  fluid  has  been 
considered  by  Kochln  [ 1937)  for  both  two  and  three  dimensions, 
Hasklnd  [1945a,  b]  has  extended  Kochln' s  treatment  to  fluid 
c£  constant  finite  depth.  The  methods  used  for  waves  generated 
by  oscillating  bodies  carry  over  with  only  slight  change,  so 
that  we  shall  not  consider  here  the  general  aspects  of  the 
theory  but  consider  instead  several  special  problems. 

20  *  ngSLffltt£-^ajaiLeyen,  bottom . 

Let  ua  first  derive  the  proper  boundary  condition  on  the 
bottom ,  Me  shall  assume  that  the  bottom  may  be  represented 
In  the  form 

aj  *  "■  lx  +  £  (p  (X )  ( 20 .  6 , 

and  that  the  fluid  flows  from  the  right  with  discharge  rate 
Q  -  C  k  •  Aa  i-n  the  derivation  of  (10.19)  we  take 

^)"C2+£A>+eV<  ()>♦■■■.  (20.7) 

Then  the  condition  that  the  bottom  be  a  streamline  is 

-c (-k+t4l'\  ...  eV'U-fc. ♦  Ck.  (20.8) 


Expansion  in  the  manner  of  section  10  and  grouping  of 


3-  **■<& 


I 

coefficient.?  leads  to  the  boundary  condition  for  : 

Y'  •  C  4  (f\  (20.9) 

We  may  hereafter  write  Y  for  t  f  ^  and  4  -for 


We  note  that  the  choice  of  L  indicates  that  the  amplitude 
of  unevenness  of  the  bottom  must  be  Mall  compared  with  A. 
for  the  linearized  theory  to  be  applicable. 

Consider  now  a  bottom  of  the  form  (see  Lamb  [1932, 
p.  409],  Wien  [1900  p . 200 ] ) 


i|  *  —  lx,  +  . 


We  look  fo  •  a  solution  in  the  form 


{( J )  u  A  Oi&v  4c  J  -f  -U**,  4.  J 


(20.10) 


(20.11) 


where  A  and  B  are  complex.  Substitution  in  (20.9),  with 
t*'  4*  Ce ok.  x  ,  shows  that  A  must  be  oure 

imaginary,  say  iA  '  ,  and  B  real,  and  furt  «r  that 

A  4.4.  —  £>  4«4  «  C  ■  / 


(20.12) 


Substitution  in  (20.3),  i.e.,  Yy  (X,q)  -  c  Y(jc,  o) 


yields 


4  £  -  A:  A! 


(20.13) 


On®  then  finds  easily  that 
//_ >  V it%  t  (  it  C&3  it i 
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it  4  4  -  4-4.  *  I 


ct0,  V  *  -? 


(20.14) 
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An  Interesting  consequence  is  that  when  c  ■£,  jtj  £  1  f 
i.a.,  when  the  flow  ;t«  subcritical,  the  crests  and  trough** 

Just  oppose  those  of  the  bottom,  whereas,  if  , 

they  occur  together..  If  C  m  t  >  there  Is  no  steady 
flow.  Also,  when  cl^j^  is  close  to  Jf  ,  it  i»  clear  that 
the  assumption  of  small  perturbations  is  no  longer  satisfied. 

By  use  of  the  Fourier  Integral  Theorem  one  may  now 
construct  solutions  for  an  arbitrarily  shaped  bottom, 
within  the  limitations  of  the  '  neorem.  For  from 


4(Jt)  *  —  j  d&  f  4  (£)  Cm  djt 

t  -«* 


(20.15) 


one  may  derive 


ffj )  *  ~  d  4  Tm}) 4 d i  (20 . 16) 

i  ^  O  -Le>  4  6<>r»A.  4  4.  -  V  litJU,  U.  A.  / 
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An  examination  of  the  asymptotic  properties  of  these 
integrals  as  JC  -«*»+■  o©  shows  that  they  do  not  vanish  if 
v  lx  “  1 4-  jc L  >  t  *  Conditions  for  the  validity  of  the 
Fourier  Integral  Theorem,  say  that  4(Ol)  is  of  bounded 
variation  and  absolutely  integrable,  indicate  that  it  applies 
to  situations  in  which  the  bottom  unevenness  is  somewhat 
localised.  Hence,  it  is  reasonable  to  require  the  additional 


boundary  condition 

•&VM-  ^(x)»0  , 


«*  oO 


Thus  we  mu® t  asi 0-1*1  d  chc  solutions  (20 * 10)  If  j  jf 
by  adding,  respectively, 


»C _ _ 

C*®  l-i.  /fv *“  V^v  J 


«o 

^ •€(}) C4o^,0(^  -  j£  t  iL)dJ  j 


(20.17) 


. AgjbL .  j"  /6(  }  )  Uw  J  j 

6#S>-  'it  #  *“  V*  *•  g>o 

where  4to  is  the  real  solution  of  4*.Lmp 

We  note  that  the  other  boundary  conditions  are  not  spoiled, 


for  the  first  expression  in  (20.17)  satisfies  (20.3)  and  its 
imaginary  (stream- function)  part  vanishes  for  Aj  *•  -  it,  so 
that  (20.3)  Is  still  satisfied. 

Thus,  if  Cx7  tJ'b/  there  is  a  local  disturbance  of 
the  fluid  in  the  region  of  unevenness  which  eventually 
smooths  out.  If  <  q'h  there  is  also  a  local  disturb¬ 


ance  given  by  (20.16),  but  as  X  -*  —  oo  there  remains  a 
disturbance  given  by  twice  the  expressions  1  t  (20.17) 

We  remark  in  passing  that  we  might  have  obtained  this 
solution  by  distributing  along  the  bottom  dipoles  of  the  form 
(13.40)  with  oC*  O  and  with  moment  density  C^»fX) 
Various  special  cases  of  may  be  considered. 

Lamb  [1932,  p.  410]  replaces  the  unevenness  by  a  single 
dipole.  Wien  [1900,  p.  202]  takes  <>(*)  =  <2.X 

and  in  the  limit  lets  t  -*  «*©  in  order  to  f^d  the 

flow  over  a  small  step.  However,  Kochin  [1938]  has  treated 
this  problem  by  a  different  method  and  finds  that  Wisti  has 
made  an  error  by  a  factor  of  two  in  the  downstream  waves 
(he  had  not  satisfied  the  upstream  condition) ,  The  flow 
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about  *  vertical  plat#  in  the  bottom  may  be  treated  by 
distributing  vortices  (13,47)  along  the  plate  with  the 
intensity  to  be  determined  by  solving  an  integral  equation. 

One  will  find  an  attractive  discussion  of  the  subject 
in  four  papers  of  W.  Thomson  (Lord  Kelvin)  [1386,  1887], 

Italian  (1906]  has  applied  the  same  method  to  three  -dimensional 
flow.  First  he  finds  the  form  of  the  free  surface  over  a 
doubly  periodic  bottom,  then  applies  the  double  Fourier 
Integral  theorem  to  construct  flows  over  irregular  bottoms. 

He  analyses  the  asymptotic  behavior  of  the  surface  for  the 
case  of  an  isolated  dipole  on  the  bottom  and  presents  graphs 
showing  the  change  in  wave  amplitude  for  different;  radial 
sectloim.  The  method  of  analysis  may  also  be  extended  to 
superposed  fluids  of  different  densities  (see  Long  (1953,  §4]). 
20  .  Flow  about  submerged  obstacles^ 

y  neariza_tj,on.  The  procedure  for  linearizing  may  be 
carried  through  in  at  least  two  ways,  leading  to  somewhat 
different  boundary  condition®  for  the  body.  Consider  a 
body  moving  in  a  fluid.  For  the  time  being  in  order  to 
achieve  somewhat  greater  generality,  we  shall  not  restrict 
the  velocity  to  be  constant.  If  the  dimensions  of  the 
body  are  sufficiently  small  compared  with  the  depth  of 
submersion,  it  will  not  disturb  the  surface  appreciably, 
and  one  will  expect  to  be  able  to  use  the  infinitesimal- 
wave  approximation.  However,  the  same  end  la  obtained  if 
the  body  approximates  to  a  flat  disc  moving  in  Its  plane, 
a  line  segment  moving  along  its  line,  a  piece  of  a 


Mliii)iniiii!|iiu  MiiiiiSfe^  . . 


l  ro 
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cylindrical  surface  moving  along  the  cylinder,  etc.,  various 
combinations  being  easily  visualized.  We  consider  the  two 
situations  separately. 

let  m0  describe  the  surface  of  a 

bounded  body  at  time  t  ,  and  let  X  be  some  typical 

dimension  of  the  body,  say  its  maximum  diameter,  and  let  A. 
be  the  depth  of  submersion  measured  to  some  point  of  the 
bod’~.  Now,  consider  the  family  of  flows  associated  with  the 
family  of  surfaces 

FCi>U,y,},0  a  F(f  ’  -  4  »  J  ,i)  •O-  (20.18) 


where  l  -  a/f^  .  As  L-*  0  the  surface  F  ^  O 

contracts  to  a  point  and  the  fluid  approaches  a  state  of 


rest.  Hence,  as  in  section  10  or  ,  it  is  allowable,  to  expand 
<f>  and  ^  in  a  perturbation  series 


<!>  -  &  * £V  y ( t 


(20.19) 


The  boundary  condition  to  be  satisfied  on  the  body,  namely, 


(£) 


F  L  <f>  +■  F  £  '  -  0 


/ 


(20.20) 


becomes 


C^Vvd  F  *  €^tkd  f^+F{  f-  t  ^CUPf.  F*  (frajrf  f. 
and  must  satisfy 


m  O 


fynO'jpL  F  *  at  <f>  '  f  F  j.  »  O  a*.  F 


Thun,  one  finds  that,  in  this  method  of  Linearizing,  the 
boundary  condition  to  be  satisfied  on  the  body  ia  the 
exact  one 


cfrcckA  F  *  *  Ft  -  o  ow  FUf^f  j,f)  *  o  .  (20.21) 

The  boundary  condition  satisfied  by  <f>  on  the  free 
surface  is,  of  course,  the  linearized  one.  The  approximation 
to  the  exact:  solution  is  better,  the  deeper  the  relative 
submergence . 

The  second  method  of  linearization  will  be  illustrated 
with  the  so-called  thin -ship  approximation.  Let  the  equation 
of  a  ship  hull  be  given  in  the  form 

>  -  ±  F(*'j)  (20.22) 

in  coordinates  fixed  in  the  ship.  Let  us  write  this  in  the 
form 


-  r(0,-  -  \ 

J-  ±  *  F  (a,*]) 


(20.23) 


where  £  is,  say,  the  beam/ length.  Suppose  the  ship 
moves  in  direction  0*  with  velocity  C  (+)  and 
consider  the  family  of  flows  generated  by  the  motion  of 
such  bodies  for  different  £  .  Let  the  velocity  potential 

be  i  L)  .  Then,  since  es  £  •*  o  the 

hull  approaches  a  flat  disc  S«  ,  the  ship's  centerplane 
section,  the  motion  of  the  fluid  will  also  approach  a  state 
of  rest  and  ’it'  may  expand 


(20. 2 -t) 
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and  similarly  for  y  .  The  assumed  forms  for  <j>  and  y 
lead  immediately  as  In  eecClon  10  o<  m  the  linearised  free- 

t  (* ) 

surface  boundary  condition  for  f  '  ,  The  exact  condition 

on  the  hull  is 


F„  (Jt  -  fet r)  dl,  <f>x !x #  jtr( x.~  jcdVA..)^;  t  F,  <f.  -  ^  (20.25) 

-  c.  (t)  F*,  ■  o. 

After  substituting  (20.23)  and  (20.24)  in  (20.25),  one 
finds  that  <£)  ^  must  satisfy 

f}  }(*J  ±  0}i)»  To(t)  Fx  U-  fc(r)  dr.  )  ,  (20 . 26) 

r  must  satisfy 

^  (*’  y  Uj  *  4  9*  ^V/- J  }  (20 .27) 

and  a  relation  of  the  form 


{fay.*  *t)-±  s«">j, 


(20.28) 


where  is  afunctional  of  the  functions  in  bracts.  Me 
note  especially  that  it:  is  a  consequence  of  the  linearization 
that  the  boundary  condition  imposed  by  the  presence  of  the 
body  is  to  be  satisfied  on  the  centerplane  section  and  not 
on  the  actual  surface.  One  will  expect  this  linearized 
theory  to  be  more  accurate  the  smaller  &  is,  l.e. ,  the 
•mailer  the  beam-to-length  ratio. 

It  is  clear  that  one  may  proceed  similarly  in  the 


situations  mentioned  earlier.  We  record  the  renults  In 
several  cases  for  reference. 


ti'S 

F:lr*t  consider  the  thln-wir  approximation  for  two- 
dimensional  hydrofoil*.  In  *  coordinate  system  C 5  Jt  j 
fixed  in  the  hydrofoil  let  the  trailing  edge  of  the  hydro¬ 
foil  be  at  (-  Q-,  -  k ; ,  and  let  the  upper  and  lower  surfaces 
be  given  by 

i|«  -  tv  +  M.  (X)  <X»\A  y  a- -<'X  i'  X  £  CO/  (20.29) 
respectively.  Define 

r(X  )  "  i  fat*)  f  ^  (*}]  ,  5(i)  -  ^  [U(x)  -  4  (Z ;]  (  (20.30) 

so  that  now  the  top  and  bottom  are  given  by 

r(x)  +4(i>/'  and  Cj*-Ui-  f{x)  -  ^(x)l  ~o>4  X  4d>f  (20.31) 

respectively.  The  class  of  profiles  in  a  form  analogous  to 
(20.23)  Is  now  given  by 

y-k  t  e[r0)(xjs  *°U)] ,  -<x*Xi<x  •  (20.32) 

It  t®  clear  that  as  L  -*  0  ,  the  profiles  approach  the  line 
segment  -iT*.  -  lv  t  04  X  4  <x  t  ao  that  the  perturbation 
procedure  is  allowable.  The  analysis  leads  to  the  linearized 
boundary  condition 

tt(x,-b±v,t)  -  -  c(i)i-'(*-fc(r)dr)j  c(^s‘(x-/^d^20.33) 

1  ’  if 

-0.4  X-fcdt  4  «,  . 

Ih*  jillLpder-hody  approximation  is  also  consistent  with 
the  linearised  free- surface  condition.  Let  the  body  be 
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defined  by 

(<j  *  tv)  +  Y  -  r  X(ji)  *<?  fxi^CK/  k?°/  U0.34) 

In  a  coordinate  system  fixed  in  the  body,  if  one  considers 
the  class  of  bodies  defined  by  Lr{"(X. )  ,  ^hen  the.  appropriate 
condition  to  be  satisfied  by  <f>  O  is 

.  f  t 

f‘ (X.-jclVdr  -Ltcr(‘,cei9llr'*~6r)rl'\x-[cdr)  -  (20.35) 

We  note  that  the  same  problem  may  be  approached  by 
two  linearized  theories.  For  example,  in  approximating  the 
flow  about  a  hydrofoil,  one  may  either  consider  It  as  a 
relatively  deeply  submerged  body  and  satisfy  the  exact 
conditions  on  the  surface,  or  consider  it  as  a  thin  hydrofoil 
and  use  the  conditions  (20.33).  Each  method  will  have  its  own 
domain  of  excellence,  but  it  is  not  proper  in  the  present 
context  to  say  that  the  thin-hydrofoil  approximation  is  less 
exact  than  the  other  one,  even  though  this  is  true  in  an 
unbounded  fluid. 

The  H- functions.  Kochin ' a  H-function,  introduced  in 
section  19c*,  may  also  be  used  effectively  for  the  flows 
considered  in  the  present  section.  The  definition  for 
three  dimensions  is  identical  with  (19.14).  For  two 
dimensions  (19.21)  is  replaced  by 

{'(?)*$ . 

Cm 


(20.36) 


irr" 


However ,  the  foniiul.ee  for  the  force  on  the  body  axe  scaewhat 
alffarent.  For  three  distant  ions  they  ere 
hr  t 

X  "  -  ~f  I H  vtc’&f  &)l  Ate' &  , 


•O  jT 


f  }  V  “  jff%  f  0)  I  it.  A 0  i{te 


(20,37) 
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whexe  \/  is  th«  displaced  volirne  of  fluid  and 


v  -  ^  /c v 


nine  twe  ■■  dimensional  formulas  are 


A-  -  vf  /H(v)f 


7*4’  I  ,  (20,3£) 

Y-nA*rcr-&fmfa*%-  ^  < 

0  —  MS> 

Re  fi  H'(o)j  -f  RtfJ-futy  jfiZ)  dJC/ 

+v  H  (v)  H(W  +  ^  ry(  tiX'ii-i-v)  H(v-tey)  , 

vdiere  A  is  the  area  of  the  profile,  (  Xc  ,  ^  c  )  are 
the  coordinate#  of  it*  centroid,  and  F  is  the  circula¬ 
tion.  Hie  renarrk*  made  In  section  19  ex'  concerning  the  use 
of  the  H -function  «j|»p;l.y  also  here. 
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"NOTICE:  When  Government  or  other  drawings,  specifications  or 
other  data  are  used  for  any  purpose  other  than  in  connection  w 
a  definitely  related  Government  procurement  operation,  the  U.S. 
Government  thereby  incurs  no  responsibility,  nor  any  obligation 
whatsoever;  and  the  fact  that  the  Government  may  have  formula! 
furnished,  or  in  any  way  supplied  the  said  drawings,  specification: 
or  other  data  is  not  to  be  regarded  by  implication  or  otherwise 
in  any  manner  licensing  the  holder  or  any  other  person  or  corpc  1 
tion,  or  conveying  any  rights  or  permission  to  manufacture,  use  o 
sell  any  patented  invention  that  may  in  any  way  be  related  then 
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Submerged  circular  cylinder.  The  appropriate  lineari¬ 
zation  for  the  circular  cylinder  is  the  one  associated  with 
deep  submergence.  Hence,  one  must  try  to  satisfy  the  exact 
boundary  condition  on  the  cylinder. 

This  problem  has  been  treated  by  Lamb  [1913;  see  also 
1932,  $  247],  Havelock  [1927,  1929i;  1936] ,  Sretenskii  [1938], 
who  considers  also  finite  depth,  Kochin  [1937]  and 
Haskind  [1945a],  who  applies  Kochin's  methods  for  finite 
depth.  Coombs  [1950]  considers  the  flow  about  a  pair  of 
submerged  cylinders,  and,  as  a  preliminary,  also  about  a 
single  cylinder;  numerical  computations  are  carried  though 
for  two  cases,  one  with  the  centers  on  a  horizontal  line 
and  one  with  them  on  a  vertical  line,,  Coombs'  method  has 
wider  applicability  than  just  to  circular  cylinders.  In 
all  the  cited  papers,  with  the  exception  of  Havelock's 
and  Coombs',  the  problem  is  solved  by  placir^;  at  the  center 
of  the  circle  a  dipole  modified  to  satisfy  the  free-surface 
condition,  i.e.,  (13.45)  with  cx  »  <?  and  f"7  ■  2.T  c  a-  } 
where  CL  is  the  radius  and  c  the  velocity  (the  c 
of  (13.45)  is  taken  as  -  i  ,  k  the  depth  of  the  center). 
This  provides,  of  course,  only  an  approximate  solution,  for 
in  the  presence  of  a  free  surface  a  dipole  in  a  stream  no 
longer  generates  a  circle,  as  is  testified  to  by  the  fact 
that  the  contour  actually  generated  is  subject  to  a  moment. 
Havelock  [1927,  1929(1  gave  second  approximations  for  drag 
and  lift  an d  later  [1938 1)  a  c om plete  so.. u ti on . 


✓ 


The  problem  nay  b«  Created  by  a  combination  of 
Milne-Thomaon'a  Circle  Theorem  [1956,  p.  151)  and  a 
formula  of  lochia's.  The  former  states  that  If  f(x) 
la  the  complex  potential  for  a  flow  with  Ite  singularities 
all  at  a  distance  greater  than  <£^  from  the  origin  and 
with  no  aolld  boundaries,  then 


(20.39) 


la  the  complex  potential  for  a  flow  with  the  same  slngulari- 
zatlon  but  nov*  with  a  circle  of  radius  a .  and  canter  at 
the  origin  situated  in  the  fluid. 

Kochin  11937]  has  proved  that  if  f(i)  la  the  complex 
potential  for  a  bounded  contour  C  under  a  free  surface, 


then 


where  C,  la  any  contour  in  the  lower  half-plane  containing 
C.  The  formula  end  ite  proof  are  almost  Identical  with  that 
given  In  (17.15).  The  first  Integral  In  (20.40)  repreaente 
a  function  regular  avenrwhere  outeide  C,  the  aecond  integral 
e  function  regular  everywhere  In  the  lower  half-plane.  If 
one  starts  with  a  function  f(x)  whose  only  singularities 
are  contained  Inclde  0,  then  the  operation 


bX/fl'  (20,41) 


whl®'  XT1®* 


IS# 


(20.42) 
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'  0  W 


+  Zi. 


yields  a  complex  potential  function  satisfying  the  free- 
surface  condition  and  having  the  same  singularities  as  ^ 
in  the  lower  half-plane. 


On  the  other  hand,  if  one  starts  with  a  complex  potential 
f(z)  whose  only  singularities  are  in  the  upper  half-plane,  then 


u  nif-fi 


(20.43) 


where 


I  lj{  1  “  fi~~ |TX  f  1  c ^  ,  a<  L  ; 


(20,44) 

will  be  a  complex  potential  for  a  flow  about  a  circle  of 
radius  A.  and  center  —i  U,  and  with  the  same  singularities 
as  f  in  the  upper  half -plane,  the  singularities  of  Yfl{{} 
being  all  inside  the  circle. 

We  start  with  the  flow  f0(p  =  '~cJ  representing  a 
uniform  flow  from  the  right;  the  free-surface  condition  is 
ratisfisd  for  o  in  a  trivial  manner.  Now  formthe  sequence 

W.  -hi  [{*] ,  4*  h  If  J,  •  ■  L.r  Kff^lvo.w 


Then  fo  +  flf  t/r  ,  •  each  rapresent  flows  satis¬ 

fying  the  boundary  condition  on  the  circle;  hence,  also 


their  eum  If  the  series  converges.  On  the  other  hand, 
j0  /,+  fx  f  +  {q ,  eech  represent  flows  satisfying  the 
free  surface  condition,  and,  hence,  also  the  sun*  if  it  exists. 

Let  us  now  consider  the  two  operators  TO  ***d  M 
Tn.  i«  always  being  applied  to  a  function  regular  and 
bounded  in  the  lower  half-plane.  Since  a  ( ttL  maps 
the  exterior  of  the  circle  onto  the  the  interior,  the  maximum 
of  I m  [f  ]  /  for  J  in  the  lower  half-plane  does  not 
exceed  that  for  ff  |  within  or  on  the  circle.  We  write 
this  as 

Imffji*  Hf  H-  i-fi  •  <20-46) 


In  particular, 


II  £  II U II  •  (20-w) 

The  operator  T°(  is  always  applied  here  to  functions 
regular  everywhere  outside  and  on  the  boundary  of  the  circle. 
Hence  may  be  contracted  to  C  and  one  can  establish  the 

following  estimate  for  j  in  the  lower  half-plane: 

<  j  [■—.  +■  Zve^^lEi^l^lOl ] 

C  -*(+-*■)  -  (20.48) 

a  [j~-£  +  £■  Ji  v  5  '  a  £  i  /w*.  If  / , 


4  K  Ilf  II 

where  in  the  second  inequality  k  must  be  large  enough  thst 
>$A  •  For  fixed  values  of  va  one  may  aalact 
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kf£l'  large  enough  so  that  K  la  as  small  as  one 

wishes ,  in  any  case,  less  than  1  .  Thus,  in  particular, 

Hfz^hKlI  fL-.IKIh,.-,  II  (20.50) 

<K‘H  U III 

tfj£  K"inj 

From  this  It  follows  easily  that  the  series 

f°  +  {•  f  fi  +■  •  •  •  +  {l«,  +  {z^r,  +  •  ■  •  (20.51) 

with  terns  defined  by  (20.45)  converges  uniformly  in  the 
part  of  the  lower  half-plane  exterior  to  I'-j+c  U  I  <  a. 

One  nay  extend  the  method  to  flows  about  tuoxe  general 
cylinders  by  combining  the  operator  ft]  with  another 
defined  in  terms  of  conformal  mapping  of  the  given  profile 
into  a  circle.  The  procedure  carried  through  above  is  a 
natural  generalization  of  the  procedure  used  by  Havelock 
in  his  first  two  papers  [1926,  192*]  to  find  the  second 
approximation.  However,  in  his  later  paper  [ 1936b]  he  used 
a  different  procedure,  one  which  has  also  been  ured  by 
Ursell  in  analogous  problems.  This  consists  in  expressing 
the  potential  as  a  sum  of  multi  poles  situated  at  the  center 
and,  of  course,  already  modified  so  as  to  satisfy  the 
condition  on  the  free  surface  and  as  X  —  °o  .  This  leads 
to  an  infinite  set  of  equations  for  the  coefficients. 

The  method  is  quite  suitable  for  approximate  computation. 


li\ 


After  computation,  of  H((t)  ,  th«  foneilaa  (20.38) 
can  be  used  to  find  the  fore*  and  moment.  In  the  computa¬ 
tion,  of  H  only  the  terns  with  odd  indices  contribute. 

This  lead*  to  a  considerable  saving  in  effort.  For  example, 
1,£  one  had  approximated  the  flow  by  the  first  three  terms 
of  (20.51)  and  computed  the  force  by  integrating  the 
pressure  over  the  cylinder,  the  result  would  be  the  same 
as  that  obtained  by  using  the  H- function  evaluated  for  f, 
alone,  and  without,  the  need  of  finding  ft  .  ...aveiock  ha® 
frequently  mad#  use  of  this  device  without  specifically 
introducing  the  H- function.  Figure  20  from  Havelock  [19368] 
shown  Rb~X  *n4  Y  plotted  in  units  of  J ^  f  with 

abscissa  f  /  VvC  for  a/k»-\  .  The  curves  labelled 

ft ,  and  Yt  give  the  result  when  only  the  first 
approximation  la  used,  i.e., 

It 1* 


H { it,) m  2  ft  c  a1  It  t 

R- Xfj.1  t[f)  (**?■)  e 

v-  -  $  -  w  -  <m . 


(20.52) 
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Computation  of  M  gives,  on  this  approximation ,  the 
anomalous  result 

b 

n-kLi  ~~  fk  1 '  "R  • 

In  the  formula  for  Y  the  terms  resulting  from  buoyancy 
and  circulation  are  omitted.  Havelock  [1928b]  ha® 
investigated  the  form  of  tie  surface  over  a  moving  dipole 
i >  °ver  a  sphere  to  the  same  degree  of  approximation. 

Sgmg_Jjjbmerged  three-dimensional  bodies.  The  flow 
about  submerged  ellipsoids  and  bodies  of  revolution  in 
general  has  obvious  Interest  in  connection  with  the 
wave  resistance  of  submarines.  As  a  result  there  is  a 
considerable  amount  of  both  theoretical  and  experimental 
work  available,  and  even  some  tables  for  computation  of 
wave  resistance.  Most  of  the  theoretical  work  does  not 
go  beyond  the  approximation  in  which  one  represents  the 
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body  by  the  singularity  distribution  appropriate  to  an 
unbounded  fluid,  but  with  the  potential  function  for  the 
singularity  modified  to  satisfy  the  conditions  on  the 
free  surface  and  at  X  -  +  00  .  Thus ,  to  find  the  flow 
about  a  submerged  sphere  one  will  in  this  approximation 
use  a  modified  dipole  with  axis  In  the  direction  Ok  and 
moment  {  C  <x,h  .  One  should  realise,  however,  that 
the  boundary  condition  on  the  body  appropriate  to  deep 
submergence  has  not  been  satisfied.  The  necessary  refine¬ 
ments  could  be  carried  through  for  the  sphere  in  a  manner 
similar  to  that  used  for  the  circular  cylinder.  Since  the 
sphere  (and  even  more  the  circular  cylinder)  is  a  poor 
shape  to  which  to  apply  perfect-fluid  theory,  such  a 
computation  is  of  only  moderate  interest.  Both  Pond  (1951, 
1952]  and  Havelock  [1952]  have  considered  methods  for 
improving  the  accuracy  with  which  the  boundary  condition 
on  bodies  of  revolution  is  satisfied.  This  is  particularly 
important  in  estimating  the  moment  about  the  transverse 
horizontal  axis,  but,  as  Pond  shows,  of  less  importance 
for  the  wave  resistance. 

Havelock  [ 1931a]  treated  by  the  approximate  method 
the  wave  resistance  of  prolate  and  oblate  spheroids  moving 
both  along  and  at  right  angles  to  their  axes.  Latsr 
[1931b]  he  extended  the  results  to  general  ellipsoids 
moving  in  the  direction  of  the  longest  axis.  Velnblun 
[1936]  has  considered  bodies  of  revolution  using  the 
slender -body  approximation,  but  satisfying  it  only  in  the 
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4pp:rox,iin«t«  sense  described  above;  his  aim  was  to  find 
form®  of  minimum  wave  resistance.  Weinblum  [ 1931]  returned 


to  the  problem,  taking  up  in  particular  numerical  computation 
of  the  wav©  resistance  for  a  given  shape.  Tables  and 
graph®  are  given  to  facilitate  the  computation  for  certain 


classes  of  bodies.  Experiments  were  made  by  Wcinblum , 
Amtsberg  and  Bock  [1936]  on  three  forms  at  several  depths, 

f °,re  Meant  experiments  exist  whose  results 

are  not  publicly  available,  A  general  survey  of  the  rhenrv 
may  be  found  in  Besaho  [1957],  survey  or  tne  theory 


If  one  has  once  computed  the  function  H  (h  i&)  for 


1  M 


a  source  and  a  dipole,  it  is  usually  straightforward  to 


compute  it  for  bodies  generated  by  distributions  of  sources 
and  dipoles,  and  hence  to  compute  the  force.  Let  S,  be 
a  surface  containing  a  single  submerged  source  of  strength 
rw  at  the  point  (a,l,c)l  L>o  (i.e.,  (13.36)  multiplied 
by  -  Kv  ) ;  one  finds 


H  ffc ,  9)  -  4 J- K, S  c  4 U(a  *  *  c  Wl 9) 


(20.53) 


For  a  dipole  of  moment  M  in  the  direction  Ox 


one  finds 

(20.54) 
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These  may  now  be  superposed  as  necessary  for  either 
discrete  or  continuous  distributions.  Thus,  if  we  write 

for  the  function  (13.36)  with  ) 

replaced  by  (JFf*j,  J  ),  and  if  we  can  express  cp  for 
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soine  problem  by 

*l<])  di  ,  (20.55) 

s 

then  (cf.  (19.20)) 

U(^).--njj&u'iek(^OK9  +  ^ , 7)dn  (20.56) 

5 

Prolate  spheroid.  We  give  an  example  of  the  preceding 
remarks.  A  prolate  spheroid  of  major  semi-axis  a.  and 
minor  semi-axes  L  moving  with  velocity  C  in  the 
direction  of  its  major  axis  can  be  represented  in  an 
unbounded  fluid  by  a  distribution  of  dipoles  of  moment 
density 

,/*(£)  ^  Ac  ( oJ'e-  -  ^  )  f  /jl<ae,  (20.57) 


where 


A  1  -  -2  too 


»  e 


/-  A1 , 


placed  along  the  major  axis  between  the  two  foci.  Hence,  with 
the  center  at  depth  ^  on*  has  in  this  approximation 


H  (4t,  9)  -  43Fi  A  c  k  t  j  <£f 


-e.t 


(20.58) 


Substituting  in  the  first  formula  of  (20.37),  one  finds 
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Figure  21  from  Havelock  11931a)  shows  a  graphical  presence- 
tlon  .if  for  spheroids  with  various  ratios  of 

«•/*  and  for  k  - z  4  • 


In  comparing  the  different  curves  one  should  keep  In  mind 
the  selected  vertical  scale;  one  based  on  displaced  fluid, 
l.a.  R/flf  ,  would  give  the  comparison 

a  different  aspect. 

As  mentioned  earlier,  it  has  been  shown  by  both  Pond 
and  Havelock  that  this  approximate  treatment  of  the  boundary 
condition  on  the  body  is  inadequate  for  computation  of  the 
moment.  Figure  22  is  from  Pond  [1951]  and  shows  the  com¬ 
puted  moment  about  the  centei  for  a  Rankine  ovoid,  i.e,,  for 
the  body  generated  in  p-\  unbounded  fluid  by  a  source  and 
sink  of  equal  strengths  moving  in  the  direction  of  their 
axis.  The  dashed  curves  show  the  result  with  the 
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Approximate  computation;  th«  solid  curves  were  computed 
by  a  method  in  which  the  boundary  condition  on  the  body  ia 
more  closely  satisfied.  The  length  £  of  the  body  is 
10.5  times  the  maximum  diameter  d  -2.4*  A  positive  moment 
i  8  CI-CMMI  * 


Slender  bod* eg .  It  is  known  that,  for  bodies  of 
revolution  given  in  the  form  (20.34),  the  slender-body 
boundary  condition  (20.35)  can  be  satisfied  in  an  infinite 
fluid  by  a  distribution  of  sources  along  the  axis  of 
strength  density 

l( *)  -  -  c  r  (*)  ■ 


(20.60) 


If  one  assumes  that  this  same  distribution  of  the 
modified  sources  (13.36)  will  satisfy  approximately 
(20.35),  then 


and 

&) 


(20.61) 
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(20.62) 


From  this  one  finds  easily  from  (20.37) 
jr/t  T 

Rm-Xm  \GX%v' j d9uc  9t  "V 


~A  -CL- 


1C  I?v1/1 


(20.63) 


where 


A* 


“fly  -«/ 

As  mentioned  earlier,  Weinblum  [1951]  has  published  tables 
allowing  one  to  compute  R  for  ^  's  representable  as 
cert/ In  polynomials.  An  earlier  paper  [1936]  consider/  -he 
minimization  of  R  among  certain  class®*  .if  polynomial 
ft  's.  Pond  [1952]  treats  the  necessary  refinements  to 
(20.61)  In  order  to  compute  the  moment. 

Let  the  equation  of  the  hull  be  given  as 
in  (20.22)  by  J  »  ±  F (*>*])  l’*  a  coordinate  systura 

moving  with  the  ship,  which  we  take  to  «ove  with  constant 
velocity  c  in  the  di.octien  Ox  ,  If  we  assume  that  a 
steady  state  hail  been  reached,  the 


condition  for 


th«  bull  Appropriate  to  the  thin  ship  epp-oxination  is, 
from  (20 , 26)  with  f(x, j,  t)  - /(*-  ct,  y,  j  )  mad  a 

chang*  to  a  moving  coordinate  system, 

ft  (x.  y  i  o)a  +  C  Fk(  x,  y)  (20 . 64) 

for  in  S0  ,  th»  c«nt*rpl«n*  section  of  the  ship  at 

rest.  For  )  not  In  Sm  on*  ha*  fy  X,  y  ±  o)  ~  o 

from  «ynr  try  considerations.  f  mutt,  of  course,  alto 
satisfy  (20.1)  and  tbs  condition  of  vanishing  notion  m  £-+  <» 


The  boundary  conditions  nay  be  satisfied  immediately 
by  distributing  sources  (1j.26)  over  £•  .  If  we  again 
denote  the  potential  function  in  (13.36)  by  i  )> 

then,  for  Infinitely  deep  fluid,  tin*  solution  is 

¥('  *■  i  y  | )  *  jj  C(  x ,  y  | ;  J,  y  °)Fjl(J\\  )  d€.  ( 20 . 63) 

'  5# 

This  follows  easily  from  known  the or  sms  in  potential 
theory  [see,  e.g,,  Kellogg,  'Foundations  of  potential 
theory,  Springer,  Berlin,  1929,  pp.  160-166].  (The  part 
of  6  regular  in  y  <0  does  not  inter- fere  with  the 
satisfying  of  (20,64)  since  the  z- derivative  of  these 
tana  van! she*  for  «  *  0.) 

The  quantity  of  chief  Interest  is  the  resistance 
MWultlnf  Iron  th*  waves.  This  may  be  computed  by  using 
again  (20.36)  and  (MM?)  (and  remembaring  to  taka 
account  of  both  halve*  of  the  hull) ,  or  by  direct 


account  of  linearization,  I.e. 

f?»  2  ?C  [Jpx^^FjX'AjjcUdLy  ■  <20.6 6) 

If  the  latter  form  is  used,  only  the  single -integral 
term  In  G  gives  a  non-vanishing  contribution.  In 
either  case  one  finds  immediately,  again  for  infinitely 
deep  fluid, 

jr/x. 

R  .  f  u2 0  [ P'(0)+  0\  OfjdO, 

ft  Vii  i€c ^(2 

P-  c*>(i>a  StoQ)<d*  dy  I  (20.67) 

S# 

Q  -  °  a  ux.9)  cix  dy  . 

The  result  may  be,  and  has  been,  put  into  a  variety 
of  different  forms  by  change  of  variable  and  der 
of  integration.  We  give  one  of  them.  Let  A  ■  «-c  ®  . 
Then  one  may  verify  that 

R  =  ffd  Id^FJl(A,^)F(Jl>i)f‘i^,(x-/)/-f(^th)j  f 

rcVS«  So  1  (20.68) 

)m  J  . . .  £  ^  C&i  ^\X  (l  X  . 

This  expression  for  R  in  terms  of  the  hull  form  and 
velocity  waa  first  given  by  Michell  [1898],  but 
derivations  by  different  methods  have  since  been  given 

i 

by  many  other,  e.g.  Havelock  [1932,  1951],  Sretenskli 
[1937],  Kcc  iln  [1937],  Lunde  [1951],  end  Tinman  and 
Vossers  (195SJ .  It  is  ueuelly  called  "Michell' s 

Integral . ” 


"Ml"’  ■  1 4  MMWtIMM*  fhlllH# 
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Because  Mlchell's  integral  gives  R  directly  in 
terns  of  the  hull  geometry  it  has  been  intensively 
investigated  by  several  persona  in  order  to  throw  light 
upon  the  influence  of  variations  of  hull  form  upon 
wave  resistance.  Foremost  among  these  investigators 
has  been  Havelock,  who  in  a  series  of  notable  papers 
11923,  1925  a,b,  1926  a,  1932  a,  b,  1935]  studied  the 
effects  of  various  systematic  variations  described  by 
the  titles  of  the  papers.  Much  of  this  work  is  summarized 
in  Havelock  1 1  0 6 1 .  In  addition,  there  are  numerous 
papers  by  G  P.  Weinbluxn  and  W.  C.  S.  Wigley  devoted  to 
comparison  of  experiment  and  theory.  One  can  find 
surveys  of  much  of  this  and  related  work,  as  well  as 
further  bibliography,  in  Wigley  [1930,  1935,  1949],  Weln 
blum  [1S50],  Havelock  [1951],  Lunde  [1957],  end 
Wehauaen  [1957].  Lunde 's  1951  paper  contains  derivations 
of  practically  all  tha  ganaral  thaoratical  rasults, 
including  tha  affect  of  finite  depth,  walls,  and 
accelaratlon.  Takao  Inul  [1954]  has  given  an  axtanoive 
survey  of  Japanese  investigations  on  wava  resistance  and 
related  topics,  end  in  e  later  paper  [1957]  e  complete  survey. 

In  order  to  allow  better  exploitation  of  Mlcheli’® 
integral  much  attention  has  been  give  in  recent  years 
to  itf  numerical  computation.  One  can  find  a  general 
dinette eion  in  flirlhtf'tf  Kotik  •  [1954’  ,  and  various 


special  pfibpiliilj  in  lafeaeiiinsfeii  U947'}*  Reinov  [1951] 


i-n 


Gullloton  [1951]  and  Weinblura  [1955].  Tha  last  two  papers 

« 

both  contain  sat*  of  tables  to  be  used  In  evaluating 
Michel  1'  s  I.,  itegral. 

In  making  a  comparison  of  the  theoretically 
predicted  wave  resistance  with  measured  wave  resistance 
one  must  examine  critically  the  experimental  method 
ror  estimating  the  wave  resistance.  The  standard  method 
consists  in  measuring  the  total  resistance,  estimating 
the  part  of  the  total  resulting  from  the  effects  of 
visosity,  and  attributing  the  difference  to  wave  making. 
Thu*  the  accuracy  of  the  experimentally  estimated  wave 
resistance  depend*  upon  the  accuracy  of  the  estimated 
viscous  resistance  and  upon  the  validity  of  the  assumption 
that  the  two  may  be  added.  In  the  case  of  a  very  thin 
body  this  estimate  can  be  made  accurately  and,  in 
addition,  the  physical  assumption  in  the  thin-ship 
linearization  ia  well  realized.  Figure  23  from  a  report 
by  Weinblum,  Kendrick  and  Todd  [1952]  show*  a  comparison 
between  estlnaa  tad  and  computed  values  off?  */(  $  $  for 

a  towed  "friction  plane"  21  feet  long  with  parabolic 
ends  and  3  foot  draft.  These  experimental  data  present 
Michell's  integral  in  a  most  favorable  light =  For 
more  ship- like  forms  the  separation  of  viscous  from 
wavs -making  resistance  ia  mora  difficult  and  tha  compared 


Figure  23 

values  s«l don  show  such  striking  quantitative  agraanant, 
although  It  Is  still  fal“  In  many  cases.  Wa  call  attention 
to  the  fact  that  K‘ chall's  Integral  predicts  the  same 
wave-making  resistance  no  matter  In  which  direction  the 
ship  moves. 

So  far  we  have  discussed  vessels  moving  In  an 
infinitely  deep  fluid.  However,  If  for  our  function  G  we 
had  taken  (13.37)  instead  of  (13.36) ,  the  same  analysis 


would  have  lead  us  to  the  following  axpraasion,  first 


given  by  Sreteaskil  [1937] i 

OO  t 

R.maf( pV;*<?W 


oe«4 


Hirt  A  l  1(1  thi  non***°  solution  of  in  m  \>  fauU*  ,u  4, 
if  such  exists,  l.«.  Iff  c'y^k  >  t  ;  otherwise 
As  k-*oo  v  and  one  obtain*  on*  of  th#  form:.  of 

Mich*  11'*  integral.. 

An  expression  for  th*  wave  resistance  of  a  thin 
•hip  moving  down  th*  center  of  a  rectangular  canal  wait 
derived  independently  by  $r#tenskii  [1936,  1937]  and 
Keldysh  and  Sedov  [1937].  'The  result  may  be  found  in 
Lunde  [1951]. 

One  may  naturally  a*k  how  the  wave  pattern 
Illustrated  in  Figure  l  for'  a  moving  source  la  related, 
to  that  for  a  ship.  In  the  thin -ship  approximation, 
the  ship  in  replaced  by  a  distribution  of  sources  on 
th,«  eenterplene,  so  that  *.ach  infinitesimal  area  of 
the  eenterplene  contribute*  to  such  a  pattern 
according  to  Its  strength.  However,  in  many  large 
vees*  i  th*  middle  part  of  the  ship  is  cylindrical , 

»o  that  <xm0  in  this  region  and  only  the  bow  and 
stern  regions  contribute  a  nonvanishing  source  density. 
Thus,  if  one  replace*  the  ship  by  a  single  source  In 
the  bow  region  and  a  single  sink  in  the  stern  region, 
the  resulting  wave  pattern  will  approximate  to  that  of 
a  ship,  the  approximation  being  better  at  higher  values 
of  the  Fronde  number  cffyL  .  Depending  upon  th*  value 

of  c/ffL  ,  the  transverse  wave  systems  from  the  two 

| 

singularities  may  either  reinforce  or  partially  cancel 
one  another.  When  they  are  in  phase,  e  Larger  amount  of 
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Figure  2i* 


energy  it  being  left  behind  in  the  wave  ays ten  end  the  re* 
si stance  curve  ihowa  i>  maximum,  when  they  ere  Qttfc  of  pCIA'ffl1®  A 
minimum,  the  so-called  "humps  and  hollows"  of  the  resistance 
curve;  these  show  clearly  in  Figure  23.  Replacing  the  ship 
by  a  source  and  sink  is,  of  course,  a  gross  simplification. 
However,  It  serves  to  explain  qualitatively  certain  aspects 
of  a  ship's  wave  pattern  and  resistance  curve,  and,  in  fact, 
can  be  given  a  certain  amount  of  validity  as  an  approximate 
computation  of  Michel 1' a  integral  for  sufficiently  large 
c/fft  .  For  very  large  value*  of  »  the  wave  length 

of  the  transverse  waves  along  the  path,  If  c.1/^  ,  becomes  much 
larger  than  L  and  one  may  approximate  the  ship  by  a  dipole. 
Many  photographs  of  the  wave  pattern  made  by  a  fast  motor  boat 
fall  into  this  class.  The  photograph  reproduced  in  Figure  24 
shows  clearly  the  diverging  waves  from  the  bow  and  stem,  a 
third  set  possibly  originating  at  the  forward  shoulder,  and 
also  the  transverse  waves,  which  presumably  are  here  near-’ 
in  phase. 

The  angular  opening  of  the  wedge  containing  the  wave 
pattern  should  be,  according  to  (13.42)  and  Figure  1,  38°  56' 
in  deep  water.  This  is  usually  confirmed  approxltiately  in 
photographs,  although  difficulty  in  determining  boundaries 
makes  precise  voufirnuttion  difficult.  For  ships  moving  in 
water  of  finite  depth  k  the  angular  opening  changes  as 
shown  in  Figure  2,  and  for  supercritical  velocities ,  i »e . 

c^JU,  ,  there  are  .-o  transverse  waves, 

Jinnaka  [195/]  has  recently  published  a  brief  survey  of 
the  theory  of  ship  waves. 


hydrofoils.  Me  take  the  hydrofoil  as  described  In 
(20.29)  and  treat  the  problem  two-dimensional  1y.  Assuming 
constant  velocity  C  and  steady  motion  and  taking  our 
coordinate  system  moving  with  ttw  hydrofoil,  the  boundary 
condition  (20.33)  becomes 

%(x,-U.±o)--Cr(x)rcsU),  -***<*■•  (20.70) 

This  problem  has  been  treated  very  thoroughly  by 
Keldysh  and  Lavrent'ev  [1936].  They  follow  a  procedure 
quite  analogous  to  that  used  in  section  19c*  to  find  the 
waves  generated  by  a  vertical  oscillator  not  in  a  wall. 
Distribute  vortices  of  intensity  if (x)  and  sources  of  strength 
6(x)  along  the  line  0l£  >Abut  taking  them,  of  course, 

modified  as  in  (13.43)  in  order  to  satisfy  the  free-surface 
condition  and  the  conditions  at  infinity.  To  start  with, 
we  take  4>(x)  «- Z  C  5  *(x )  .  It  then  follows  from  the 

theorem  of  Flemelj ■ Sokhotskii  (cf.  17.18)  that 

£j-U+o)~  *fy(jc  L-<?)  --2c  s'(x),  (20.71) 

a  step  toward  satisfying  (20.70).  We  now  look  for  a 
complex  potential  in  the  form 

where  we  have  separated  the  source  and  vortex  potentials  in 
(13.43).  The  boundary  condition  (20.70)  now  yields  an 

integral  equation  for  Y(>*)  in  much  the  same  manner  that 
(17.18)  was  derived: 

(20.73) 

U  ([- 2c s'(7 )£(*- 

*  —  ^  £  Jt  £  • 


!!« . ,'L 
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Noting  from  the  third  expression  in  (13.43)  that  f,  and^f 
are  functions  of  the  difference  JC  -  J  ,  we  define 


X  x  - 2  i 


1  o  ;i^fV 

nd'2ue  /frzi 


2  nfv  (jc-  ^,-ilv)  «  £  + 


-  2.  vd 


-tVJt  (  t 


K(x)  +  i  l(x). 


Here,  for  example t  ^ 

to*)  -i  +  . 

The  integral  equation  (20.73)  can  now  be  written  :* n 


the  following  form: 

a 


)KU  -  $)dl--2JT  c  r(x)  +  \lcs(  ()H(X- 1 74) 


where  the  right-hand  side  Is  a  known  function. 

This  integral  equation  is  the  hydrofoil  analogue  of  the 


thin-wing  integral  equation  of  airfoil  theory: 


JlfiPjTJ  ~ -ZJTc  r’(jc)  ■ 


(20.75) 


In  the  latter  equation  the  kernel  is  simpler  and  in  addition 
only  the  function  describing  the  camber  and  the  angle 

1 

of  attack  occurs  on  the  right  side.  Since  the  wing  thickness 
does  not  enter  into  the  determination  of  %  in  (20.75)  it 
nay  be  neglected,  for  only  ^  is  needed  to  find  the  lift. 


The  situation  is  clearly  different  for  hydrofoils.  Even  a 
symmetric  wing  with  zero  angle  of  attack  may  nave  a  circulation, 
and  hence  lift.  This  is  a  consequence,  of  course,  of  the 
presence  of  the  free  surface  and  the  associated  wave  motion. 

Kochin  [1936]  has  also  considered  hydrofoils,  but  from 
a  somewhat  different  viewpoint.  He  has  essentially  used 
the  "deep -submersion"  linearization  described  first  in 
this  section.  Thus  he  must  satisfy  the  exact  boundary 
conditions  on  the  wing  as  well  as  the  Kutta- Joukowski 
condition.  However,  one  cannot  say  here,  as  one  could 
for  an  infinite  fluid,  that  his  method  is  more  exact  than 
that  of  Keldysh  and  Lavrent’ev.  Their  approximation  is 
more  accurate  the  thinner  the  wing,  for  a  given  submersion. 
Kochin' s  Is  more  accurate  the  deeper  the  submersion,  for 
a  given  wing. 

Equation  (20.74)  is  not  sufficient  to  determine  /(x) 
uniquely.  One  must  still  add  some  further  condition. 

We  shall  assume  a  finite  velocity  at  the  trailing  edge, 

*  -a 

Keldysh  and  Lavrent'ev  propose  solving  the  integral 
equation  (without  actually  proving  that  a  solution  exists) 
by  expanding  K,  H  and  Y  in  a  power  r  -it*  in  T*  and 

then  determining  recursively  the  coefficients.  Let 

KJzvLJT '(if 

w  t.  .it  " 


(20.77) 


wirw  >fi' 


. . 


**e 


'IIHflllW  mk#  1 1  iiiWH’  Hi  %  Wi  I  He-  '  m*m**wi 


/ 


l! 

■I' 


ZW 


The*»  (20,74)  gives  the  following  sequence  of  Integral 
equation*. 

/  sJf )£f’~2,cr'M  I 

-a,  } 


a 


■in.fzcs'(ruj, 

-  Ol 

m  0 


(20.78) 


to 


This  procedure  has  the  obvious  advantage  of  reducing  the 
solution  to  the  airfoil  integral  equation  for  which  an 
explicit  solution  satisfying  the  traillng-edge  condition 
is  known.  If  m  denote  temporarily  the  right-hand  sides 
of  the  equations  (20.78)  by  f^(^)  ,  respectively, 

then  the  general  solution  is  [see,  e.g.,  W.  Schmeidler, 

I 

Integralgleichungen  .  .  .  ,  Akademische  Verlagsgesellschaf t , 


SI 

lJX« 

n 

d" 


must  remain  finite  for  Jfr  .  We  assume  $  /-*') 

finite,  ^  may  possibly  have  a  singularity  of  the 

form  i  /faT^X 
Integral 


near  jd  e  -  a. 


However,  the 


i  r 
>  * 


./fe  &  d.\  -  /  JC»  ”*  ^ 

3L  ^ 

^c?r  />  0 


.Oo< 


is  a  polynomial  in  Jt  *  +  V  .  Thus  the  last 

two  summands  of  (20.80)  remain  finite  at  xm~o^. 
However,  the  first  summand  potentially  contributes 
terms  like  4, 

/; 


n 


-L  \l  x^-  o} 

Irt  order  to  avoid  this  singularity  at 


select  the  total  circulation 

cx* 


fl ~d> 


X  *  -  cu 
so  that 


we 


(20.81) 


Substituting  into  (20.79),  one  finds  finally 

fuvftT  -fh  ■  <20-82) 


- Os 


^(x'j  itself  is  given  by  the  sum  displayed  in  (20.77). 
Although  the  singularity  at  the  trailing  edge  has  been 
removed,  chere  is  still  one  at  the  leading  edge;  this 
occurs  also  in  thin-air foil  theory  and  corresponds  roughly 
to  the  fact  that  the  conditions  of  linearisation  (l.e.  of 
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small  dls turbine*)  are  not  satisfied  near  the  leading 
edge  stagnation  point. 

Keldysh  and  Lavrent’ev  compute  the  Integrals  which, 
will  be  necessary  If  rfa)  and  sfjc,)  are  given  as 
polynomials  and  apply  their  computational,  method  to  a 
flat-plate  and  circular-arc  airfoil  «t  a  small  angle  of 

m  f*  #"  dk  #*  lr 
f!E  L  Li  IEL>iE\  • 

In  order  to  find  the  force  and  moment  on  the  wing, 
it  is  convenient  to  fall  back  on  the  H- function.  'One 
finds  easily 

10. 

(20.83) 

K(t.)  f  4  cls'(  $  s'(jc)Jc*  U(  f~x)t 

+  l  cfy(j)6U)- (QOtf  f)]%i>AQ-x)jd^dx, 

Formulas  (’0.38)  allow  one  to  complete  the  calculation 
for  special  cases. 

The  theory  analogous  that  that  described  above 
for  fluid  of  finite  depth  4.0  has  been  carried 
through  by  Tikhonov  (1940],  He  has  applied  the 
method  to  calculate  the  lift  and  drag  coefficients  for 
a.  flat  plate  at  a  nail  angle  of  attack. 

Rather  than  reproduce  the  graphical  presentations 


of  Keldysh  and  Lavrent'ev  and  Tikhonov  for  the  flat  plate,  we 
shall  give  Instead  the  lift  and  drag  coefficients  for  a  sub¬ 
merged  vortex.  Here  one  may  give  relatively  simple  analytic 
expressions,  and  the  qualitative  behavior  of  the  curves  is 
similar  to  that  of  a  flat  plate.  The  formulas  for  lift  L  and 
drag  D  are  as  follows: 


Lo-trc-gr*1 


(20  *84) 
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where  r>  .f  is  the  real  root  of  "4^.  >  ^4.,.  For  finite  depth 


the  expression  for  D  stems  from  the  last  term  in  (13.4/ ) .  The 
dimensionless  coefficients  4- -  DC.  jfr 1  and  CLm  fi~  ?cT)  i^/frZ 
are  shown  in  Figure  25a  for  infinite  depth  as  functions  of 
and  in  Figures  25b  and  25c  as  fctnctions  of  c ’L/c}&.0  for  various 
values  of  A  ,  Foi  ^finite  depth  CL  starts  with  a 
value  I A lr  and  tends  asyt  totically  to  -iMtt  ,  crossing 

the  axis  at  c  /%lm2.41.  For  finite  depth  the  coefficients  have 
a  discontinuity  at  *  1  .  As  c.x/qI  r»o  C.-+a  t  and  as 

C.  //( .j  /  »  £p  -*  /*.  fb  ( I ~  ft>')  •  For  c/ ^.4,*  ^  1  ,  is 


always  negative  and  increasing  with  a  vertical  asymptote  at 


C  *  1  end  a  horizontal  on_  as  cv/<^^0  -*•  «o  «t 
,  these  curves  start  at  ^ 
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Further  development  of  hydrofoil  theory  ha*  taken 
piece  In  several  directions.  Haekind  11945a]  has  extended 
Kochtn ' s  "deep-submersion"  theory  to  water  of  finite  depth. 
However,  he  does  not  discuss  the  steps  necessary  for  fulfil¬ 
ment  of  the  Kutta-Joukowskl  condition,  ae  does  Koehin. 

The  lifting-line  theory  for  airfoils  of  finite  span  has 
been  extended  to  hydrofoils  by  Mu  [1954],  Sreslin  [1957], 
and  Hasklnd  [1956],  Parkin,  Perry  and  Mu  [1956]  and 
Lai tone  [1954,1955]  have  investigated  both  theoretically 
and  experimentally  the  effect  of  bringing  a  given 
hydrofoil  so  clot  to  the  surface  that  the  infinitesimal - 
wave  approximation  breaks  down  completely.  There  exists  also 
a  considerable  amount  of  work  on  flow  about  cavitating 
hydrofoils.  However,  since  the  effect  of  gravity  is 
neglected;  this  work  is  not  considered  in  the  present 
article.  Experimental  data  relevant  to  the  theoretical 
development  outlined  above  is  scanty.  Reports  by  Benson 
and  Land  [1942]  and  by  Land  [1943]  give  results  of  an 
experimental  investigation  of  the  effect  of  depth  of  sub¬ 
mersion.  However,  the  investigations  were  not  designed  to 
test  the  validity  of  the  theory  and  do  not,  for  example, 
include  the  region  of  maximum  C p  •  Ausman  ] 1953] , 
in  connection  with  an  experimental  investigation  of  the 
pressure  distribution  on  the  upper  surface  of  a  hydrofoil, 
measured  the  lift  coefficient  and  compared  It  with 
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that  predicted  by  the  thin -hydrofoil  theory.  The  theory  felled 
»h«n  too  —U  b.c.u.0  the  ...oct.ted  free  eurf.ee 

over  the  hydrofoil  no  longer  approximated  i  u  if  1.  n  !&  t  s  i-Tn  I,  'w av  6  b  , 
or,  in  other  words,  the  thin  hydrofoil  was  not  thin  enough  for 
these  values  of  for  the  theory  to  be  applicable.  It 

should  also  be  emphasized  that  for  small  value®'  of  ^  l*,  /c 
the  occurrence  of  cavitation  on  the  upper  surface  must  be 
taken  into  account  for  a  complete  theory. 

A  comprehensive  survey  of  hydrofoil  theory  Is  given  In 
a  recent  paper  by  Nlshiyama  [1957}, 

20  Y  •  Planing  surfaces. 

The  following  discussion  Is  limited  to  two-dimensional 
motion,  for  the  theory  of  three-dimensional  planing  surfaces 
for  flows  with  gravity  do^s  not  appear  to  have  been  developed. 

For  the  linearized  problem  it  is  natural  to  consider  the 
planing  surface  or  glider  as  an  approximation  to  a  flat  plate 
moving  along  the  surface  of  the  undisturbed  fluid,  i.e.  the 
curvature,  angle  of  attack  and  vertical  displacement  are  all 
assumed  small.  In  order  to  formalize  the  perturbation  pro¬ 
cedure,  let  the  planing  surface  be  represented  by 


y  =  k  -+  F(x)  ,  )*|  $  j  F  (-a.)  =  Oy 


(20.85) 


. . . 
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in  coordinates  fixed  in  space,  and  let  the  fluid  have 
velocity  —  c  at  «£  *•  -f  •  Thua,  we  are  going  to 
consider  the  flow  to  be  a  perturbation  of  a  uniform 
flow. 

First  let  us  consider  briefly  in  a  qualitative  fashion 
the  exact  solution.  There  will  be  a  stagnation  point  A 
somewhere  behind  the  loading  edge  and  a  jet  td.ll  be 
thrown  out  ahead  of  the  glider.  We  take  it  to  be  of 
thickness  b  and  to  make  an  angle  with  OX.  If  $  -  -  CX  +  T(^>  J ) 

and  *¥  *  **  u  y  f  Y  (x,  are  potential  and  stream 
function,  respectively,  we  shall  take  the  free  surface 
ahead  of  the  glider  to  be  ;iven  by  f  «  -be  and 
behind  the  glider  by  f  -  0  (see  Figure  26). 


^=0 


2fL./ 


Figure  26 

Then  b/a,  and  Ail  a  will  all  be  functipns  of  v 

and  .  It  will  be  assumed  as  one  of  the  boundary 

conditions  of  the  problem,  in  analogy  with  the  Kutta- 

! 

Joukowski  condition,  that  the  velocity  is  continuous 
at  the  trailing  edge.  It  is  obvious  that  the  flow 
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near  the  leading  edge  cannot;  conform  to  the  requirement 
that  it  be  a  small  perturbation  of  a  uniform  flow. 

However,  we  shall  give  arguments  below  to  indicate  that, 
except  In  the  neighborhood  of  the  leading  edge,  this 
effect  is  of  the  second  order. 

In  order  to  get  some  idea  of  the  relative  size 
of  the  jet,  consider  the  simpler  problem  of  a  flat 
plate  of  length  J  and  angle  of  attack  PC  gliding  on  a 
weightless  fluid.  This  problem  can  be  solved  exactly, 

(see,  e.g. , M llne-Thomson  [1956,  $  12.26];  A.E.  Green 
[1935,  1936]).  The  asymptotic  expression,  for  small 

jT  ^  ^ 

<X  ,  of  both  the  ratios  b/i  and  AL/f  s  ^  |  +•  * 

i.e.  they  are  both  of  the  second  o^der.  We  shall  suppose 
that  this  relation  continues  co  nold  when  gravity  is 
acting. 

We  now  carry  through  the  perturbation  procedure  of 
section  10  (V  (see  especially  (10.16)),  writing 

j>  --cjc  +  £  f  +  ■  •  •  ,f  -  -  C-  r  6  •  •  •,  /  i~"  ' 

U^-F(oc)  -  t  F°  U)  i  £  k(l>+-  •  •  •  ,  (20.86) 

Substitution  in  the  exact  boundary  conditions  then  yields 


the  following  linearized  conditions: 


1 

flHWff! 


*\b  0  h  . 
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{j£  |  7  ^ 


(20.87) 
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The  free  surface  Is  given  by 

Y^W)-  %  <SVu,°),  U  Da. 

3 


(20.88) 


Ws  require  at  usual  that  the  disturbance  vanish  as  X“*a&  , 
One  will  explet  that  the  behavior  near  the  leading  edge 
will  reflect  in  some  manner  the  inconsistency  of  the 
exact  solution  with  the  notion  of  a  small  perturbation. 

It  will  turn  out  that  it  will  be  necessary  to  allow1  a 
singularity  at  the  leading  edge.  (Almost  the  same 
situation  exists  in  the  thln-hydrofoil  or  thin -wing 
theory  since  the  stagnation  point  near  the  leading  edge 
also  prevents  the  flow  in  that  region  from  being  a  small 
perturbation  of  a  uniform  flow.) 

A  singularity  at  the  trailing  edge,  although 
mathmatically  possible,  has  been  specifically  proscribed. 
The  strength  of  the  singularity  and  the  elevation  h  of 
the  trailing  edge  will  be  determined  as  functions  of 
f  c  t  th®  course  of  solving  the  problem, 

The  problem  as  formulated  above  has  been  considered, 
for  infinite  depth,  by  Sretenskli  [1933,  1940],  Sedov 
[1937],  Kochln  [1938],  and  Maruo  [1951].  Haskind  has 


extended  Sad.ov'a  analysis  to  finite  depth  [1943a],  and 
later  [1955]  has  treated  a  glider  moving  on  a  wavy 
surface.  Yu.  S.  Chaplygin  [1940]  has  apparently  carried 
through  a  fairly  comprehensive  numerical  analysis  for 
a  flat  plate  waking  use  of  Sedov's  method  of  analysis 
I:*##  p. '83]).  Sretenskli 'b  paper® 
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ara  expounded  In  tens  of  a  flat  plate,  but  the  method  clearly 
has  wider  applicability,  as  he  remarks  In  his  first  paper. 
Sretenskil 1 s  1940  paper  gives  the  results  of  rather  extensive 
calculations  for  flat  plates.  Mario's  paper  is  conceptually 
very  similar  to  those  of  Sretenskil,  but  his  method  is  not 
quite  as  efficient  for  computation.  However,  he  also  gives 
computational  results  and  Includes  a  correction  to  take  ac¬ 
count  of  the  failure  of  the  linearised  theory  near  the  lead¬ 
ing  edge . 

Both  Sedov  and  Kochln  introduce  the  complex  potential 
■f())  *  f  tlY  and  thereafter  the  function  fW' 

Although  the  two  methods  are  not  by  any  means  the  same,  they 
have  much  In  caramon  with  the  treatment  of  hydrofoils  given 
above.  Consequently,  we  shall  outline  below  the  method  fol¬ 
lowed  by  Sretenskil 

As  e  preliminary  we  need  a  result  from  section  21  below. 
Suppose  that  a  pressure  distribution  p>(Jt  )  ,  which  we  take 
to  be  absolutely  integrable,  is  given  on  the  free  surface. 

Then  the  complex  velocity  potential  must  satisfy 


and  the  free  surface  Is  given  by 

-j  »  f*.}  *  £  Y  (X,0)  - 


,.jj/  . 11  + 

.  . ■  •ttHpRil1;  K 

i  ij^nnfj  r,  i  ■!  I*"  'i|. 

_ 


(20.89) 


(20.90) 
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The  function  which  satisfies  (20.89)  and  'which 


vanishes  n»  £-*«©  cun  b*  written  In  several  forms, 

of  which  we  select  the  following  (see  eq. (21.3*1 )): 

<?  *®  ~i>( J-f)  7°  -tf/'l-/) 

«i)  -fTf-  j  <*  hi  V  Pi ~\-y  fl  jf(  3  f( 

-oO  o  ~oO 

The  fret  surface  Is  given  by 


£(20.91) 


«"  -hMm  UZZ 

?  c  —oO 

the  reason  for  leaving  y,  explicitly  in  the  formulas  will 
appear  below. 

When  a  gilder  is  moving  on  a  free  surface,  the 
streamline  ^sC'^(* ,o)  will  consist  partly  of 
free  surface,  where  p(x)  «•  0,  and  partly  of  the 
wetted  surface  of  the  glider,  where  f(x)  is  some 
unknown,  function. 

Equation  (20.92)  may  then  be  written  as  the  following 
Integral  equation  for  this  unknown  function  (x) : 

L  +•  F(jl)  -Xvw  f  f  ®  €h  d\  f- 

Y~°  °  (20.93) 

>  W<*- 

\  V1  *  fi* 

Once  f(x)  has  been  determined,  one  may  substitute 
back  into  (20.92)  in  order  to  find  the  fona  of  the  free 
surface  for  |Jt  I  >  , 


Zjf(l)  «>  <20 

0  V  +  -1— _  fp(f)ii~)>(x 


J  'll'1 
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It  la  possible  to  work  directly  with  (20.93),  and 
this  Is  the  procedure  followed  oy  Maruo.  However, 
Sretcnskii  differentiates  twice  with  respect  to  x  and 
adds  V1,  times  (20.93).  This  yields 


(X,  oil  j. 

f4(*>+VlF(x)*v1t-  fd^(f)j(X^)c0»)i(^t)t  JdX 


Although  this  last  equation  is  a  necessary  condition 
for  p(x) ,  It  obviously  cannot  determine  it  uniquely, 
for  the  last  term  of  (20.93),  assuring  vanishing  of 


the  disturbance  far  ahead  of  the  glider,  was  lost  in 
the  formation  of  (30.94),  Thus  one  still  has  need  for 


(20.93).  Equation  (20.94)  is  essentially  the  aquation 


derived  by  Sretenskil. 

Let  us  now  -Integrate  (20.94)  with  respect  to  x 


f.  m  x 


to  x,  and  denote 


*  57 T_1 


Then  equation  (20.94)  becomes 


(20.95) 


IA> 

F  (*)-  F(-<^ f vv/  F( (jd\  r  . -v  at 


ffl 


II*"  . 


in* 


Where  an  additive  constant  has  bean  discarded  since  h 
Use  If  is  an  undetermined  constant.  Equation  <20. 95) 

Is  Just  Frandtl ‘ •  Integro-dlfferential  equation  for 

the  circulation  about:  an  airfoil  of  finite  span  [see,  e.g., 

N. I .  Muskhellrhvili  »  Singular  integral  equations, 

Noordhaff,  Groningen,  1953,  ch.  171.  Thus  known  method® 
for  solving  the  airfoil  equation  can  be  carried  over 
to  the  study  of  this  equation.  However,  the  solutions 
themselves  -snnot  be  taken  over  directly,  for  different 
boundary  conditions  are  imposed:  in  the  airfoil  equation 
the  unknown  function  is  the  circulation  P^jt)  and  it  is 
usually  assumed  that  P(-jt)  •  V {X)  and  f  (-&)  •  I  (&■)"  &  \ 
in  the  present  problem  Pi*)  is  not  necessarily 

symmetric  and  F  (-a)  *  P'  (-a)  *  0;  but  I  (a)  is  not 
restricted  except  to  be  finite.  The  theory  of  the 
Frandtl  lntegro-differential  equation  without  the  customary 
additional  requirement®  aseociated  with  airplane  wings  has 
bean  developed  by  L.G.  Magnaradze  l Soobehch ,  Akad.  Nauk 
Gruzin,  SSI  J(1942),  503-508]. 

The  equation  can  be  solved  by  an  extension  of  Glauert's 
tacthod  [t-g-j  H.  Glauert,  'The  element*  of  airfoil  and 
airscrew  theory,  Cambridge,  1943,  ch  XI].  This  is  the 
method  which  has  been  used  by  both  Maruo  and  Sretenskii. 
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results  a  equations  of  the  form 

vtw  -  Q,  CL  at  R,vtv+  S,  , 


V  k  -  y,  0*>  t  K ,  V  W  +  0 ,  ,  (20,99  > 

F,(-a)-Q1a.t  R*vL  +  St  , 

where  Qi  ,  R  i  ,  Si  are  function#  of  VO/  ;  these 
equations  may  be  used  to  determine  viv  and  a,  as  function# of 
VOw  *  As  long  as  Q  there  will  be  a  singularity  at 
the  leading  edge. 

Once  p(x)  has  been  determined  approximately,  one  can 
compute  the  lift,  drag  and  moment  about,  say,  the  center. 

To  the  order  of  approximation  appropriate  to  the  linear 


theory  they  are 


!..  ■  J  p  t  “t') « t  x 


i 

P  m  j  p(.*)  (U)  (kX) 


(20.100) 


ou 

p'l  m  I"  p(x)X 


For  the  flat -plate  glider  it  ia  possible  to  give  the  following 
asymptotic  expressions  for  these  quantities  when  VCX  -»  O . 


.101) 


L‘  R-*L 

M .  4  r  a*f  cV[(  ■ -  *  #] + 1 9  a< h va')  ■ 


me-e  were  first  given  by  Sedov,  but  are  also  derived  in 
the  papers  by  Koeliin  and  Sretenskll. 


Figure  2?  reproduces  several  of  Sretenskll ' s 
computed  pressure  distribution!?  far  a  flat  plate, 
lb*  predictions  of  the  linearized  theory  cannot,,  of  course , 
be  expected  to  be  Accurate  near  the  leading  edge.  Maruo 
has  corrected  his  computed  points  in  this  region  by  using 
the  exact:  theory  for  *  weightless  f  uid.  Both  Maruo  and 
Sretentkii  give  further  computational  results  which  we,  do  not 
reproduce . 
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21.  Wiwmw  resultiii^^  pressure  distributions 

In  the  situations  considered  up  to  no*  in  this  chapter  the 
pressure  at  the  free  surface  ha*  been  taken  as  constant,  We  now 
consider  the  result  of  allowing  the  pressure  over  the  free  sur¬ 
face  to  he  a  given  function  of  both  position  and  u.  ?.  Other¬ 
wise  the  fluid  is  taken  to  be  infinite  in  horizontal  extent  and 
to  be  either  infinitely  deep  or  of  uniform  depth  lu  ,  The  time 
variation  in  pressure  will  be  limited  to  two  cases.  In  section 
21  oC  a  periodically  varying  pressure  i®  considered;  in  section 
21  /S  the  pressure  is  taken  to  move  with  uniform  velocity; 
section  21  give®  some  references  to  a  combination  of  these  two. 

In  section  22  waves  from  pressure  distributions  will  be  considered 
again  in  connection  with  Initial-value  problems.  Since  the  meth¬ 
ods  for  finding  the  velocity  potential  are  similar  in  most  res¬ 
pects  to  those  used  in  finding  the  velocity  potential  for  a 
source,  we  shall,  with  one  exception,  give  the  results  without 
proof. 

Just  as  in  the  cases  of  the  stationary  source  of  periodic 
strength  and  the  moving  source  of  constant  strength  treated  in 
section  13  <jf  ,  we  must  In  the  present  situation  Impose  boundary 
conditions  at  infinity  in  order  to  ensure  a  unique  solution. 

Tbs  ii&poised  conditions,  ssMly  the  radiation  condition  and  tbs 
vanishing  of  the  fluid  motion  far  ahead,  respectively,  are 
•eleptad  as  being  physically  reasonable.  However,  one  may  pro- 
M»d  dl If erectly,  derive  formulas  analogous  to  (13.50)  and  (13,51) 
and  tb«i  fled  the  Mlit  »:#  t  #'  «•©  .  The  resulting  velocity 
ootemtiai®  eutoaatically  satisfy  the  correct  b*u*d*rv  conditions 


at  Infinity.  Thin  method  has  been  used,  for  example,  by  G.  Green 
[1043]  In  the  two-dimensional  problem?  considered  In  the  follow¬ 
ing  two  sections,  and  also  by  Stoker  [1053,  1954]. 


The  theory  of  wave  generation  by  pressure  distributions  has 
an  obvious  application  in  oceanographic  problems.  However,  the 
theory  was  apparently  first  developed  In  an  attempt  to  explain 
the  wave  pattern  produced  by  a  ship.  We  shall  not  attempt  to 
disentangle  the  history  of  the  subject.  For  the  material  covered 
in  section  21  £  we  call  attention  to  a  survey  by  J.K.  Lunde 
[1951|>  ]  which  also  contains  a  useful  bibliography. 

21  •  Py«««ure  distributions  periodic  in  time 

^ree^dl^ensions^.  The  boundary  conditions  have  already 
been  given  in  section  11.  If  and  o  Are  represented  by 


?(■*, pi* 


••  u  e  T 


(21.1) 


thon  the  condition  on  the  free  surface  may  be  written 


^(*,0,))-  f  pr(x'p' 


(21.2) 


and  the  form  of  the  surface  Is  given  by 


l<*,  J-0  --  p(x,})J 
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In  addition,  a  radiation  condition  i»  assumed  at  Infinity  (sen 
(13,9),  «q.  5)  and  a  condition  appropriate  to  Hit  depth  of  fluid. 


We  shall  also  assiuM  p(j£f^)  to  1m  absolutely  integrable. 

The  velocity  potential  can  be  expressed  as  follows; 

•".ao  ^ 

tfif  .  T I  ,  "  t  ^  I  f  t  r  T  \  ,  T  ,  ,  n,/  Itl  1  {}  rf  t‘>7  “  ^  V  6 


-<*' O 


(21.3) 


+ ifr  if fH(v  ^  -  f >l*(}  -tf)  << 1  *  t,  -a*, 

•  -OO  -«o  # 

and  in  cylindrical  coordinates  ,  j  *  in  the  fonn 

/ft,*  J,  t  fptRyiR'dZ'A*' 

«  0 


+  <*  v  e*  ^  / j  pfR.'* 7  J0^  jfcVR' -  28-R‘,c*a(«'-«)  J  &  dfc 
"L  f  ^  J  j 


Ifl  crO 


(21.4) 
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The  addition  theorem  for 


[see 


Watson,  A  treatise  on  the 


theory  of  Bessel  function®,  Cambridge,  1944,  p.  333]  allows  on© 


to  write 


...  ..  «• 

4^;  M  O 

tm  «  I  ,  "  2  /  ^  ;  •  (12.5) 
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If  p  Is  independent  of  <K  ,  one  may  derive  easily 

OO  ^ 

fi'K  y) «  ft* (l R‘)R' d  K 1  rv M^)U^  RJ) d  k 

*  a 

+  S±Zp-  J,  (v  <?.)  fp{tZ')J.(:<R')R’ci/l'  . 


(21. 6) 


The  asymptotic  form  for  large  2  of  (21,6)  is  a  relatively  simple 
expression: 


~  lih  ^ * "J>  ft&'MWM1 . 


(21.7) 


We  note  in  passing  that  the  potential  function  (21.3)  or  (21.4) 
can  also  be  obtained  as  a  distribution  of  sources  on  the  surface 
(see  Hudimac  [1953,  p.  78]).  This  may  be  easily  verified  as 
follows.  In  (13.17")  let  ^-  0.  Then,  using  (13.12),  one 
obtains  (substituting  f, $  for  o,c  ) 

2.  nf°^-e^  d(i  \l(x-f)+(} -  S)1-  at i-i  2/v/1 


A  distribution  of  these  sources  over  the  plane  ^-0  of  strength 

yields  (21.3)  (we  recall  that  a  source  of 
strength  m  behaves  l.-.ke  -m/r  near  the  singularity). 

The  rate  at  which  the  pressure  distribution  does  work  upon 
the  fluid  can  be  calculated  directly  or  by  us lug  equation  (8.2). 
Consider  the  volume  of  fluid  contained  in  a  large  cylinder  of 
radius  R0  .  Then,  from  (8.2),  after  appropriate  linear  iant  1  OO' 'i 
the  rate  of  incre  tse  of  energy  of  the  fluid  is  given  by 
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dt 
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+ 1 
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Now  substitute  (21.1)  and  take  the  average  over  a  period,  which 
will  clearly  be  aero.  The  result  may  ee  written: 
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2  / /  f  fj  R 
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8) 


The  first  integral  gives  the  average  rate  W**  at  whioj  the 
pressure  distribution  Is  working  upon  the  fluid.  It  must  equal 
mluus  the  second  integral.  If  p(J2,o<  j  «  *>(%-)  •  then  we  may 

apply  (21.7)  to  obtain  a  relatively  simple  expression  for  the 
average  rate  over  the  whole  fluid; 
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(21.0) 


To  carry  through  the  computation  when  p  la  not  circularly  sym¬ 
metric  Is  more  complicated  arithmetically,  but  can  be  carried 
through  by  use  of  (21.5). 

One  can  find  an  investigation  of  the  waves  resulting  from 
a  doubly  modulated  pressure  distribution  over  a  rectangular 
domain, 
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In  a  paper  of  Sretenskli  [1956]. 

If  the  fluid  is  of  uniform  depth  /v  »  the  expression  for  the 
velocity  potential  is 


|  m, mi  1 1  in  a  mu  i -T  -in-  -nniinL  mr—  \ 


(21,10) 
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where,  as  usual,  nv0  is  the  real  solution  of 

/\V,  0  k.  -V  -  o  . 


Other  forms  of  this  expression  similar  to  (21.4),  (21.6)  and 

(21.7)  can  be  found  with  no  difficulty.  We  give  only  the  analogue 
of  (21.9) : 
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(21. 11) 


The  identities  iollowing  (13.18)  may  be  used  to  put  both  (21.10) 
and  (21.11)  into  other  forms. 

Two  dimensions.  The  derivation  of  the  velocity  potential 
will  be  carried  through,  at  least  in  part,  since  it  illustrates 
a  nice  application  of  tl.e  Plemelj-Sokhotakii  formulas.  Two  com¬ 
plex  units  will  be  introduced,  as  described  at  the  end  of  section 
11.  That  is,  w©  shall  write 
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(j>'x i y  t/m  ft (*> ^ )  o*i 6t  t  fi  ( Jt, y )  sUm.  6i  «  12  e.j  f 

.1 

p,f x  t)  =  p((x)oii6<  f /z>z(>x)u~6t  *  pe  ^ 


(21.12) 


and  also  Introduce  a  stream  function  T m  ^  ^2.  and  a 


p  1,  ©  x  po  i  ©  n  t»  1  fli  1 


(21.13) 


Then  the  boundary  condition  on  the  free  surface  may  bj 

formulated  as  follows 

Ibx.j  lO)i  Lvf  (Ji  ~10)J  *  “  “  £  ' 


(21.14) 


The  definition  of  ^  *  f  V  tvf  may  be  extended  to  the  whole 
plex  j  -plane  by  reflection,  i.e. 

(?  I  X  +  <■  W  )  *  0  (vt  **  l,  u  }  >  I  2 
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Then  the  condition  (21.14)  may  be  written  in  the  form 


I^[t'(x±  oo)  +iv{(x±  io)  j-±  ^  ^p(Jt) 


( 2  X e  15) 
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I 

We  shall  suppose  to  be  absolutely'  intsgrabls  on  the 

■A 

Infinite  interval.  In  addition,  we  shall  suppose  that  either 
pfj t)  satisfies  a  Holder  condition  (and  is  hence  continuous)  on 
the  whole  infinite  interval,  or  else  that  there  are  a  finite 

number  of  segment*  (  -  «»,  4 )  "d  2,^2)  *  •  *  *  *  ,  °° )  * C,  ^  ^  <  1  » 

such  that  pfx)  satisfies  a  Holder  condition  on  any  closed  seg¬ 
ment  interior  to  one  of  the  above  segmen  ts^  and  at  an  end-point 
may  be  expressed  in  the  form 

fW-,%  I  ,  C-^orl,  , 

1  (x  -  o) 

where  ^(Ot)  satisfies  a  Holder  condition  at  the  end.  Here  a 
Holder  condition  moans  that  for  any  pair  of  points  JC ,  ,  tXj»  ,  pfx) 

satisfies  ,  . 

lp(Jc‘)-r('x'-)l *  A  ,  a>0 


In  the  first  case  -f  will  be  assumed  to  have  no  singularities 
In  the  whole  lower  half-plane.  In  the  second  case  the  behavior 
of  •f  (’j)  near  an  end-point  C  will  b®  restricted  so  that  it 
must  satisfy 

As  usual  It  will  be  assumed  that  1/7  is  bounded  as  j -»  00  and 
that  only  outgoing  waves  are  generated. 

The  solution  of  this  boundary- value  problem  for  the  function 


{ V  i i*  determined,  up  to  are  Additive  real,  constant 
which  may  be  discarded  here,  by  [see,  e<g*,  Muskhelishvili , 
Singular  integral  equations,  Noordhoff,  Groningen,  1953,  §  £  43 

78] 


After  integrating  the  differential  equation  and  selecting  the 
solution  i  is  to  represent  outgoing  waves  at  ,  one  ob 
tains  finally  [the  derivation  is  similar  to  that  of  (13.28)] 
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(21.17) 


where  the  path  of  Integration  for  r  is  ta.ien  in  the  lower  half¬ 
plane.  The  asymptotic  form  of  the  time -dependent  velocity  poten¬ 
tial  is  given  by 


i.  e-(v)  T a)j as 

—  oO 


(21. 18) 


arid  the  asymptotic  form  of  the  free  surface  by 


yjf,0  ~  I«jj| 
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(21. 19) 
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From  thlf  last  expression  one  can  easily  derive  the  average 
rate  at  which  the  pressure  ays ten  is  transferring  energy  to  the 


fluid: 
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The  expression  for  •({%)  can  be  put  into  several  different 
forms  by  changing  variables  and  deforming  the  path  of  integra¬ 
tion  appropriately.  Thus,  if  one  introduces  a  new  variable  A 
hy  k  (  if  )  =  "  A  ( ~6)  and  deforms  the  resulting  path  to  the 


X  -axis,  one  obtains 
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(21.21) 


A  different  deformation  of  the  path  leads  to 
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For  fluid  of  depth  ^  an  expression  for  the  complex  velo¬ 
city  potential  analogous  to  (21.10)  and  (21.21)  is 
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On©  will  find  both  the  two- and  the  three-d  linens  Iona  i  case 
of  a  periodic  pressure  distribution  over  infinitely  deep  water 
discussed  in  Lamb  [1004,  pp.  387-303].  Stoker  [j.957,  ch.  4] 
discusses  in  considerable  detail  the  two-dimensional  problem 
of  waves  generated  by  a  periodic  uniform  pressure  applied  over 
a  finite  interval. 

21/S  .  Moving  pressure  distributions. 

In  this  section  we  shall  suppose  that  a  fixed  pressure 
distribution  is  moving  with  a  constant  velocity  c-  Thus  the 
motion  may  be  treated  as  time-independent  in  a  coordinate  system 
moving  with  the  pressure  distribution.  The  boundary  condition 
at  the  free  surface  is  given  by  [see  (11.3)] 


and  the  form  of  the  free  surface  by 

7^3)  “  ^  ^  f°  (x>) )  ■ 


(21.24) 


(21,25) 
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In  addition,  w©  shall  assume  vanishing  of  the  fluid  motion  far 
anead,  i.e.  as  ,  and  the  usual  conditions  appropriate 

to  infinite  or  finite  depth,  p  will  be  assumed  to  b©  absolute¬ 
ly  integrable  and  to  vanish  for  sufficiently  large  values  of 
X  %  ;  however,  the  latter  condition  can  be  weakened. 


Results  will  be  given  without  proof  since  their  derivations 

are  similar  to  those  in  13  £  and  21  «  .  The  results  for  two 

and  three  dimensions  will  be  separated. 

Three  dimensions.  The  expression  for  the  velocity  poten¬ 


tial  for  infinite  depth  of  fluid  can  be  given  as  follows: 
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(21.26) 
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The  rate  at  which  the  pressure  distribution  is  transferring 
energy  to  the  fluid  is  given  by 


DO 


-  oO 


(21.271 


This  may  be  conjputed  directly  from  (21,26).  The  first  term 
gives  no  contribution  since  it  is  an  odd  function  of  x  -£  (of. 
the  evaluation  of  20.66).  The  final  result  may  be  expressed  as 
follows : 
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If  the  pressure  distribution  Is  given  in  cylindrical  coordinates, 
p~r(R,<x)  ,Jt=E^cx  ,  >  then  one  may  express,  sa’  , 

p(J-)  in  the  form 

OO  2rl  f  A  -J 

P(>)  mfdZf  dot  Z  ft  %->  *)  Cm e*  [v  Z  -  ^)J 


o  o 


(21.29) 


and  a  similar  formula  for  Q  ('J')  .  If  f  depends  only  upon 
then  Q  (-J- )  —  0  and 


PfJ^)  or  pfiflfiiJl  (y  K’  maHA)  i  1? . 


(21.30) 


If  the  fluid  is  of  depth  ^  ,  the  velocity  potential  is  given 


by 
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is  the  positive  real  root  of 
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The  rat©  of  transfer  of  energy  may  again  be  computed  from  f21.27) 
and  again  only  the  second  integral  gives  a  non vanishing  contribu¬ 
tion.  The  result  may  oe  expressed  in  several  forms  analogous  to 
(21.26) -(21. 30) : 
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In  cylindrical  coordinate®  formulas  (21,29)  and  (21.30)  carry  over 
to  the  present  situation  with  V  replaced  by  iz<>  • 

The  asymptotic  form  of  tue  free  surface  for  either  infinite 
or  f  <  nite  depth  is  much  more  complicated  to  analyze  than  for  the 
stationary  periodically  oscillating  pressure  distribution  of  the 
preceding  section.  Although  it  is  not  strictly  necessary  to  do 
so,  it  has  been  customary  in  this  type  of  analysis  to  consider 
the  special  case  of  a  'concentrated  pressure  point".  To  derive 
the  velocity  potential  for  the  pressure  point  consider  the  pres¬ 
sure  distribution  defined  by 
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Substitute  in  (21.26)  or  (21.31)  and  take  the  limit  as  k.  „ 
Then  (21.26)  becomes 
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The  velocity  potential  for  the  pressure  point  in  fluid  of 

* 

depth  is  derived  similarly .  The  equation  representing  the 

surface  may  now  be  obtained  immediately  from  (21. 25) : 


finite 

free 


(21.35) 


The  velocity  potential  (21*34)  is  very  similar  to  that  of  a 
submerged  source  in  steady  motion  (see  1.3,36)  and  the  method 
sketched  in  section  13  for  the  derivation  of  the  asymptotic  ex¬ 
pression  (13.42)  can  be  carried  over  directly  to  the  moving 
pressure  point. 

The  result,  expressed  in  cylindrical  coordinates,  is  as 

#  1  1  AUJ'JCfi  " 

•X  Tiitir  x  x  Hjf  w  wp  * 


for 


0  ^  &<  <  Jf  ~  Arc  ^  *  oC 


at  Q[ 


o(wrl) 


} 


for 


o(  m  oK  f 

v 


,  DC  q  J  * 


) 


}o.VV  'V£  ^ 

■I  -  9  3  I 

V?  % 


A 


(21.36) 


3  F~(%) 


for  <X  c 


£  £>r 


<  J 


f  ??  )<X  ) 


ISC 


p 

if 


V 


1 


IJT 


I 


(A 


;z> 


.  £ 


1 7S  *  f“c  v, *. 


S/2 


J. 


i 

s* 


i 


,.*•*  .  uwwymmtwtiioni 


. ~  . . 


<IY 


IliUlulMlIMlIli 


,*jiuiuLlh;J  U imimj  y  lull  till 

Iflttln  W  pi!  vmi  H  Iff  » 


4i. .  ♦  fttflWHSHIHWHH*  >■**»»<■« M*  <  ‘ '  *  '  *"  * 


for  0<  a »  JT 


yf  R,  r)  -  - 


0 


Ff<j 


umrnAmmmmm mmmm 


tJT 


V  f*. 


f  V  rv*  *f* 


v  ^  Jl 

r'"** 


The  variables  t  /<-,  »  /■‘t  ar«  th0  8aine  on«s  defined 


following  (13.42),  where  certain  properties  are  also  given. 

Figure  1  is  not  quite  accurate  as  a  description  of  the  wave  crests 


in 


region  |T-oc  j  <£  oc  c  for  the  pressure  point  since  the 


phase  in  (21,36;  nas  been  shifted  by  *  JT  ;  the  wave  crests  in 

C 

Figure  1  should  be  moved  back  a  distance  ~  T  • 

The  wave  pattern  resulting  from  a  moving  pressure  distribu¬ 
tion  has  been  the  subject  of  many  investigations,  starting  ap¬ 
parently  with  Kelvin  [1906].  Hie  aim  was  to  explain  the  typical 
wave  pattern  found  behind  a  ship.  The  procedure  is  quite  rea¬ 
sonable  as  a  method  for  obtaining  a  qualitative  prediction  of  a 
ship’s  wav©  pattern,  since  a  moving  »ihip  has  associated  with  it 


a  pressure  distribution  around  the  wetted  hull. 


obvious 


disadvantage  of  the  method  is  that  it  gives  no  connection  be¬ 
tween  the  geometry  of  the  hull  and  the  wave  pattern.  For  this 
the  ’thin  ship’  approximation  of  section  20  A  is  better  within 
its  rang©  of  applicability,  The  single  pressure  point  can  be 
taken  to  rep.  ©sent  approximately  a  ship  moving  at  high  speed 
(more  accurately,  at  high  Froude  number  C/fla  •  where  L  is 
the  length),  say  a  fast  motor  boat. 

For  detailed  invostigation  of  the  asymptotic  expression  on© 
should  refer  to  Hogner  [1923],  Peters  [1949],  Bartels  and 


Downing  [  1955]  (who  do  not  restrict  themselves  to  a  pressure 
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point)  and  Stoker  [1957,  ch.  8].  The  necesuary  eodiilc.  ions 
for  finite  depth  have  been  made  by  Havelock  [  11(03)  aud  Te-  >ro 
Inul  [1936]  and  are  describeu  qualitatively  in  the  discussion 
following  (13.42) .  One  can  find  an  exposition  of  the  theory  of 
naves  generated  by  moving  prticuff  distributions  io  a  report  by 
Lund©  [ 1951b] •  Several  papers  by  Havelock  [ 1909,  1914b,  1919, 
1922]  take  up  the  wave  resistance  ( -  W/c)  of  n  pressure  distri 
button.  Hogue r  (1928]  has  also  considered  the  wave  reallowance 


and  gives  essentially  (21.2®). 

In  the  preceding  considerations 


have  assumed  that,  as 


/ 


X ,  ,  x 

jt  4-  ’  -*►  t»o  * 


1 


pfjt,'})  approached  aero  sufficiently  quickly  so 


that  It  might  toe  represented  as  a  Fourier  integral.  It  is  also 

L  ‘  x 

possible  to  proceed  somewhat  differently,  assume  periodic 

in  one  or  both  variables  and  use  a  Fourier  series  representation. 


This  has  been  done,  for  example,  by  Volt  [1957a],  who  has  con¬ 
sidered  for  both  infinite  and  finite  depth  a  moving  pressure 
distribution  of  the  following  form: 


4P& 

m  (  P\j)  £  i  /w  fe.  X  /  I  V  ^  *"  u 

i  3  i  /J  I  >  . 

The  waves  resulting  from  a  pressure  point  moving  parallel  to 
beachas  forming  angles  of  30°  and  45°  have  been  treated  by  Hanson 
[ 1926] ;  in  the  same  paper  he  also  treats  the  waves  formed  by  a 
pressure  point  moving  over  a  two— 1 ayt red  fluid.  A  detailed  inves 
tigatiou  of  this  last  topic  is  gl  on  in  a  paper  of  Sretenskii 's 

[1934]. 
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Two  dimensions.  By  introducing  a  stream  functi 

mmmumz. 


>n  T  ^  ) 

and  a  complex  potential  {  *  V*  ^  k  T  »  the  free  surface  boundary 
condition  can  toe  put  into  a  form  analogous  to  (21.14),  namely, 


- 1 


Ri  |{  (Jc~  L0)+ tv 
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(21.37) 
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In  addition,  we  assume 


-*■  and  also 


I  /  |  bounded  as  J 
|f'|o.  We  shall  assume  f(  *  )  subject  to  the  same  limits 

>1  «o 

tions  as  in  section  21  oc  . 

One  may  apply  the  same  method  of  analysis  to  derive  the 
following  forms  for  the  complex  velocity  potential: 
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where 


taken  in  the  lower  half-plane.  The  rate  at  which  the  pressure 
distribution  transfers  energy  to  the  fluid  is  easily  found  from 
formula  (21.27)  and  the  second  expression  for  -fQ)  to  be 
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If  the  fluid  is  of  depth  fu  »  the  complex  velocity  potential 


may  be  given  in  a  form  analogous  to  the  second  formula  of  (21.37) 
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where 
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is  the  real  positive  root  of 
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the  last  ter 


and  exl<  Im  only  if  v4,  *  J 4  jc, 1  2  /  ;  if  *v 4  *>  f  , 
in  (21. 40)  must  be  deleted.  The  rate  at  which  the  pressure  i 
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doing  work  upon  the  fluid  to  given  by 
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(21.41) 


The  absence  of  the  second  term  in  (21.40)  and  the  vanishing 
of  W  when  vL  <  1  correspond  to  the  absence  of  an  Infinite  train 
of  trailing  waves.  A  similar  Phenomenon  occurs  behind  a  moving 
singularity  in  two  dimensions  (cf .  (13. 46) -(13. 48)  and  the  follow 
ing  remarks]. 

For  either  (21.37)  or  (21.39)  the  form  of  the  free  surface 
can  be  written  down  immediately  from  the  formula 


(21.41) 


We  shall  not  carry  out  the  details.  The  asymptotic  form  of  the 
surface  behind  a  two-dimensional  "pressure  point",  or  also  a 
distributed  pressure,  is  much  easier  in  two  than  in  three  dimen¬ 
sions  and  we  again  omit  a  detailed  statement.  However,  the 
problem  has  been  treated  by  itt  Ivin  [1905]  and  is  discussed  in 
Lamb  [1932,  242-5],  both  for  infinite  and  finite  depth.  It 

is  also  discussed  in  the  paper  of  Peters  [1949]  already  cited  In 
connection  with  the  three-dimensional  problem. 

Derivations  of  the  complex  velocity  potential  may  be  found 
in  the  papers  of  Sretenskii  [1934,  1940],  Bedov  [1936],  Kochin 
[1939]  and  Hasklnd  [1943a]  already  cited  in  connection  with 
plating  surfaces.  We  refer  also  to  papers  of  Dean  [1947]  and 

T Iranian  and  Vossers  [1955]. 

20  Y  .  Moving  periodic  pressure  distributions . 

mmmmmmmmmmmmmmmmmmhmimrni  . . 

It  is  clearly  possible  to  combine  the  cases  considered  in 
sections  20  <%  and  20  A  and  consider  the  wavt.  resulting  from  a 
pressure  distribution  expressible  in  the  form 


f(x  ,  t)  •  p^-ct  ^  T  p%(&-  ct^)  t  J 


where  the  coordinates  are  fixed  in  space.  The  resulting  velocity 

potential  will  be  analogous  to  (13.52)  for  the  three-dimensional 

case,  if  one  is  dealing  with  a  "pressure  point". 

We  shall  not  reproduce  the  formulas  here.  However,  tue 

analogue  of  (13.49)  for  pressure  distributions  may  be  found  in 

the 

the  cited  report  of  Lunde  [1951b],  and  f r on  thisArequired  velo¬ 
city  potential  may  be  found.  For  two-dimensional  motion  the  de¬ 
tails,  carried  out  by  this  procedure,  may  be  found  in  papers  by 
Kaplan  [1957]  and  Wu  [1957]. 
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